{'fraction_non_alphanumeric': 0.03999905829174122, 'fraction_numerical': 0.03048780487804878, 'mean_word_length': 4.1406268909597, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Exploring GLIMPSE Bubble N107: Multiwavelength Observations and Simulations

Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of numerical simulations, we aim to estimate the bubble age and other parameters which cannot be derived directly from observations. Methods. From the observations of the HI (I-GALFA) and 13CO (GRS) lines we derive the bubble's kinematical distance and masses of the atomic and molecular components. With the algorithm DENDROFIND, we decompose molecular material into individual clumps. From the continuum observations at 1420 MHz (VGPS) and 327 MHz (WSRT), we derive the radio flux density and the spectral index. With the numerical code ring, we simulate the evolution of stellar-blown bubbles similar to N107. Results. The total HI mass associated with N107 is 5.4E3 Msun. The total mass of the molecular component (a mixture of cold gasses of H2, CO, He and heavier elements) is 1.3E5 Msun, from which 4.0E4 Msun is found along the bubble border. We identified 49 molecular clumps distributed along the bubble border, with the slope of the clump mass function of -1.1. The spectral index of -0.30 of a strong radio source located apparently within the bubble indicates nonthermal emission, hence part of the flux likely originates in a supernova remnant, not yet catalogued. The numerical simulations suggest N107 is likely less than 2.25 Myr old. Since first supernovae explode only after 3 Myr or later, no supernova remnant should be present within the bubble. It may be explained if there is a supernova remnant in the direction towards the bubble, however not associated with it.


Introduction
Bubble N107 (SIMBAD: CPA2006 N107) is one of the largest bubbles in the catalogue by Churchwell et al. (2006). This catalogue is based on the infrared observations made with the Spitzer Space Telescope, (GLIMPSE survey, Benjamin et al. 2003) and contains more than 300 bubbles found in the interstellar medium (ISM) near the Galactic Plane. The catalogue was later expanded by Simpson et al. (2012) to more than 5000 bubbles, which were identified by a community of over 35 000 volunteers. These bubbles, with typical radii less than 10 pc, are examples of shell-like structures -features found in the ISM all over our Galaxy. Other examples are H i shells in catalogues by Heiles (1979Heiles ( , 1984; Ehlerová & Palouš (2005), with sizes up to 3500 pc; or dusty "loops" discovered in the IRAS full-sky survey by Kiss et al. (2004) and Könyves et al. (2007).
The varying sizes and masses of these structures suggest different mechanisms responsible for their creation. These mechanisms, which do not have to act solely, include: These different models of stellar feedback were also discussed by Tenorio-Tagle & Bodenheimer (1988); Ehlerová & Palouš (1996) and Efremov et al. (1999).
The shell-like structures can trigger star formation, since they consist of a layer of cold and dense material. In the socalled collect-and-collapse scenario (Elmegreen & Lada 1977;Elmegreen 1994;Whitworth et al. 1994), the shell fragments and creates a new generation of stars. Radiation driven implosion (Duvert et al. 1990;Lefloch & Lazareff 1994) considers the compression of preexisting condensations (globules) by the pressure of the ionised gas. A brief review of these and other triggering mechanisms is given in a study of GLIMPSE bubbles by Deharveng et al. (2010).
In this paper, we present a multiwavelength study of one of the largest GLIMPSE bubbles: N107, which lies in the central plane of the Galactic disc (l ≈ 51 • .0, b ≈ 0 • .1). A multiwavelength, colour picture of N107 is shown in fig. 1.
The paper is organised in the following way: In section 2, we complement the Spitzer infrared observations with H i line observations from the I-GALFA survey (Koo et al. 2010), 13 CO (J = 1−0) line observations from the GRS survey (Jackson et al. 2006) and radio continuum observations at 1420 and 327 MHz from the VGPS (Stil et al. 2006) and WSRT (Taylor et al. 1996) surveys. From the line observations, we derive the bubble's LSR 1 radial velocity (RV), kinematical distance and masses of the molecular and atomic components. From the radio continuum observations, we derive the flux densities and corresponding spectral indices of two radio sources lying apparently inside the bubble. In section 3, we use a recently published algorithm DENDROFIND (Wünsch et al. 2012) to decompose the molecular gas associated with the bubble into individual clumps and derive the slope of the clump mass function. In section 4, we use the numerical code ring (Palouš 1990;Ehlerová & Palouš 1996) to simulate the bubble's evolution in order to estimate its age, energy input and the size and mass of the molecular cloud from which the bubble evolved. Finally, a discussion and conclusions are presented in sections 5 and 6.

Observation of Bubble N107
2.1. Data Sources GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire, Benjamin et al. 2003) and MIPSGAL (Carey et al. 2009) are two complementary infrared surveys mapping the Galactic disc with the Spitzer Space Telescope. GLIMPSE used instrument IRAC (Infrared Array Camera) which has four channels centred around 3.6, 4.5, 5.8 and 8.0 µm, for which the point response function FWHM is between 1 . 7 and 2 . 0. MIPSGAL used the instrument MIPS (Multiband Imaging Photometer for Spitzer) with a broadband channel centred around 24 µm, for which the point spread function FWHM is 6 .
I-GALFA (Inner-Galaxy ALFA Low-Latitude H i Survey, Koo et al. 2010) is part of a large project called GALFA-H i (Peek et al. 2011), which plans to map the entire sky observable by the Arecibo Observatory in the H i line. The FWHM beam size of the Arecibo telescope at 1420 MHz is 3 .35 with the GALFA-H i radial velocity channel width of 0.184 km/s. The data are provided in FITS cubes of the brightness temperature T b . The single-channel standard deviation of the gaussian noise σ noise is about 0.25 K. 1 LSR = local standard of rest The GRS (Galactic Ring Survey, Jackson et al. 2006) is a 13 CO (J = 1−0) survey focused on the Galactic molecular ring. The survey used the FCRAO 14-m telescope, with a FWHM beam size of 46 and a radial velocity channel width of 0.21 km/s. The data are provided in FITS cubes of the antenna temperature T * A , with a typical noise σ * noise ≈ 0.13 K. We converted the data to the main-beam brightness temperature: T b = T * A /η MB and σ noise = σ * noise /η MB , where η MB = 0.48 is the main beam efficiency (value suggested by Jackson et al. 2006).
VGPS (VLA Galactic Plane Survey, Stil et al. 2006) is a 1420 MHz continuum and H i line survey of the Galactic disc based on the interferometric observations done with the VLA (Karl G. Jansky Very Large Array). We use only the continuum observations since for the H i, data from I-GALFA are available with significantly lower noise. The FWHM beam size is 1 and the gaussian noise is 0.3 K, for the continuum observations.
WSRT (Westerbork Synthesis Radio Telescope survey, Taylor et al. 1996) is a 327 MHz continuum survey of the Galactic plane utilising the WSRT interferometer in the Netherlands. The FWHM resolution is 1 × (1 / sin δ) in the right ascension and declination (δ) direction, respectively. This yields for N107 (δ ≈ 16 • ) the FWHM resolution of 1 ×3 .6. The median gaussian noise is 2.5 mJy/beam.
The UKIDSS (UKIRT Infrared Deep Sky Survey) project is defined in Lawrence et al. (2007). UKIDSS uses the UKIRT (United Kingdom Infrared Telescope) Wide Field Camera (WF-CAM; Casali et al. 2007) and a photometric system described in Hewett et al. (2006). The pipeline processing and science archive are described in Irwin et al. (2004) and Hambly et al. (2008). We used data from the 7th data release.

Dust Component
The bubble's outline is most prominent in the 8 µm maps from the GLIMPSE survey, where it was discovered by Churchwell et al. (2006), who derived its basic properties: mean position of l ≈ 51 • .0, b ≈ 0 • .1; mean radius of 11 .4 and a mean thickness of 2 .3. Most emission at 8 µm is coming from the bubble edges, forming a ring-like structure. The 8 µm channel is dominated by the emission of polycyclic aromatic hydrocarbons (PAHs) (Draine 2003;Reach et al. 2006, fig. 1;Pavlyuchenkov et al. 2013). PAHs, which are destroyed in the ionised regions, are tracing the photodissociation regions (PDRs), edges of molecular clouds illuminated by the UV radiation from a massive star. The UV radiation pervades the PDRs and excites the PAHs herein.
Counterpart emission at 24 µm was observed in MIPSGAL (Carey et al. 2009). Contrary to the 8 µm continuum, emission at 24 µm is filling part of the bubble's interior. The 24 µm continuum is dominated by the emission of very small grainsdust grains larger than PAHs -heated by a nearby massive star (Pavlyuchenkov et al. 2013).
A dark hole is present in the southwestern part of the bubble's interior. Almost no emission from dust, and also no emission of 13 CO or H i (see below), is observed there. Fig. 1 shows a multiwavelength, colour picture of N107 composed of the 8 µm and 24 µm continuum emission and the 13 CO line emission integrated over the radial velocities of 38.5 to 47.6 km/s.

Atomic Component
We searched the H i data cubes by eye for features possibly associated with bubble N107. An H i shell, morphologically similar Fig. 1. Bubble N107 -false colour multiwavelength image. Red corresponds to 24 µm continuum, green to 8 µm continuum and blue to 13 CO (J = 1−0) line integrated over the radial velocities of 38.5 to 47.6 km/s. Note that the 24 µm continuum (red) is dominated by the emission of very small dust grains heated by a nearby massive star. The 8 µm continuum (green) is dominated by the emission of PAHs -large molecules found within molecular clouds. The emission of PAHs is tracing the photodissociation regions (PDRs) -edges of these clouds illuminated by UV radiation from a massive star. The CO line (blue) traces the cold gas of molecular clouds. The distinct red-green structure seen in the lower-left in the direction of the bubble's opening is likely not related to the bubble, since the radial velocity of the associated molecular material (≈ 60 km/s) differs significantly from the radial velocity of the N107 complex (≈ 43 km/s). to N107, is located around the central LSR radial velocity of 43 km/s. Note that an incomplete ring of CO clumps is present at similar radial velocities (see sec. 2.4). The H i emission is located along the bubble edges, protruding outside farther than the CO emission, forming an atomic envelope of the whole structure. Maps of the H i brightness temperature in several velocity channels are shown in fig. 2.   Fig. 3 shows a map of the galactic longitude versus radial velocity (l-v diagram). The front wall is very faint, located at a radial velocity of ≈ 34 km/s. The back wall is more distinct, located at a radial velocity of ≈ 50 km/s. Between these radial velocities, a hole is apparent inside the bubble. The relative expansion velocity of the front versus back wall is ≈ 16 km/s.  In order to measure the H i mass associated with the bubble, we define the following 3D aperture (aperture 1, see also fig. 4, left panel): radius of two times that of N107, r ap1 = 2 · r N107 = 22 .8, and the radial velocity channels from 32.7 to 52.0 km/s inclusive. All the pixels lying within this aperture contribute to the measured H i mass.
To derive the mass of the neutral hydrogen, we use the optically thin limit of H i line radiative transfer. Then, the H i column density N H i can be computed from the formula (Rohlfs & Wilson 1996): where T b is the observed brightness temperature of the H i gas and, in our case of an I-GALFA data cube, the integral takes form of a sum over a set of pixels with dv = 0.184 km/s. The Galactic H i emission has a strong background component. To estimate it, we assume that the bubble's interior is evacuated, so the H i emission observed apparently near the bubble centre is due to this H i background, not associated with the bubble. We defined another 3D aperture with the radius half that of the N107 and the radial velocity channels from 39.1 to 43.7 km/s (see fig. 4). This aperture encloses a volume in the central part of the bubble. The mean T b within this region is 72 K, which we adopt as the H i background component and subtract it before computing the mass.
Furthermore, in the process of mass derivation, we need to know the distance to N107. Considering also the 13 CO emission (see below), we adopt the central radial velocity of 42.8 km/s. The Galactic rotation model of Brand & Blitz (1993) gives for that radial velocity two kinematical distances: 3.6 and 7.1 kpc. We adopt the near one of 3.6 kpc, which yields the bubble mean radius of 11.9 pc (11.4 ). The near distance is favoured by the estimated virial masses of the molecular clumps found along the bubble borders. For the far distance of 7.1 kpc, we got 22 of 49 clumps with supervirial masses, while for the near distance of 3.6 kpc, we got only 2 clumps with virial masses above 1 (clumps 7 and 21 in tab. A.1). Note that Churchwell et al. (2006) suggest distances of 4.7/6.0 kpc. They use the same rotation model as we do, however with a less reliable determination of the radial velocity obtained from the observation of coincident H ii regions. They advocate also the near kinematical distance.
To estimate the uncertainty of the measured mass, we use the standard Taylor formula: where f is the measured value (in our case the H i mass), σ f is its uncertainty, which is derived from uncertainties σ x i of variables x i on which f depends (in our case e.g. the distance to N107).
The resulting H i mass associated with N107, i.e. the total H i mass within aperture 1 after subtracting the background, is 5.4 × 10 3 M . The major sources of uncertainty are the distance to N107 and the background emission. Assuming the distance uncertainty of 10 % (0.36 kpc) and the background uncertainty of 1 K, the total H i mass uncertainty is about 20 %, i.e. 1 × 10 3 M .

Molecular Component
Around the radial velocity of 43 km/s, CO emission is present along the bubble borders. Maps of the CO emission ( fig. 5) show that CO gas forms clumps, which are distributed in a ring-like structure along the edges of the bubble, protruding outside. Almost no emission is coming from the interior. Around the radial velocity of 48 km/s, the CO emission merges with another complex of molecular clouds, while in the channels below 38 km/s, CO emission vanishes.
To measure the associated molecular mass, we use two 3D apertures (see fig. 4, right panel). Aperture 1 is same as that for measuring the H i mass. Aperture 2 is smaller, with radius r ap2 = 1.36 · r N107 = 15 .5, radial velocity channels from 38.8 to 46.9 km/s inclusive and covers only the immediate vicinity of the bubble, i.e. only the mass accumulated in the expanding shell. Note that for comparing the observation and simulations, we use the molecular mass found within aperture 2. We treat the radiative transfer in the line of 13 CO (J = 1−0), described in Rohlfs & Wilson (1996). First, we derive the optical depth: where T ex is the excitation temperature, which we assume to be 20 K. T b is the brightness temperature of the CO line; T CMB = 2.7 K is the temperature of the cosmic microwave background and where h is the Planck constant; k is the Boltzmann constant and f = 110 GHz is the 13 CO (J = 1−0) line frequency. Assuming local thermodynamical equilibrium, characterised by a single excitation temperature T ex , we can compute the 13 CO column density from the formula (Rohlfs & Wilson 1996): where, in our case of the GRS data cube, dv = 0.21 km/s and the integral becomes a sum over a desired range of pixels. From the 13 CO column density, the H 2 column density is derived assuming the abundance (Binney & Merrifield 1998) N H 2 = 6.5 × 10 5 N13 CO . Finally, the total mass of the molecular matter m mol is computed using the standard solar hydrogen abundance X: and m H 2 is the mass in the form of the molecular hydrogen. Note that we use term molecular matter or molecular mass to refer to the matter found in molecular clouds, i.e. a mixture of H 2 , He, Article number, page 5 of 15 A&A proofs: manuscript no. sidorin_etal_2013_n107 CO and trace amounts of other molecules, elements and dust. The distance to N107 (3.6 kpc) is discussed above in section 2.3. The total molecular mass within aperture 1 is 1.3 × 10 5 M , from which 9.8 × 10 4 M in in the form of H 2 , rest in helium and heavier elements. The total molecular mass within aperture 2 (in the expanding shell only) is 4.0 × 10 4 M , from which 2.9×10 4 M is in the form of H 2 . The major source of uncertainty is the 13 CO abundance. We used N H 2 = 6.5×10 5 N13 CO (Binney & Merrifield 1998). Observations and models show, however, that the abundance varies for different clouds and depends also on the physical conditions inside individual clumps (Sonnentrucker et al. 2007;Visser et al. 2009;Wilson 1999). Another source of the uncertainty we consider is the distance to N107. Using eq. 2 and assuming the 13 CO abundance uncertainty of 20 % and the distance uncertainty of 10 % (0.36 kpc), the uncertainty of the molecular masses given above is about 30 %.

Radio Continuum and Spectral Index
Since an H ii region was identified in the direction of the bubble (Paladini et al. 2003, region 576) and also the 24 µm emission is present inside the bubble, we expected to observe its radio continuum counterpart.
Both the VGPS (1420 MHz) and WSRT (327 MHz) surveys show a distinct diffuse radio source located apparently within the bubble, in its eastern (left) part, slightly below the area, where the 24 µm emission is strongest (see fig. 6). Beside this source, the bubble interior is filled with much fainter radio emission. Another diffuse radio source is present further east of the bubble (l ≈ 51 • .35, b ≈ 0 • ). This source is associated with a complex of molecular material found at radial velocities of about 60 km/s, so it is likely not associated with the bubble. Fig. 7 shows a comparison of the VGPS radio continuum with the emission at 8 µm and 24 µm. Inside the bubble, the radio continuum follows the bubble outline defined by the PAHs emission at 8 µm. In the northern (top) and western (right) part, the radio continuum is mixed with the emission of very small grains at 24 µm, while in the eastern part (left), the 24 µm emission vanishes and only the radio continuum is observed.
We measured the radio flux densities in four regions: two located apparently within the bubble -defined by apertures A and B -and two located outside the bubble -defined by apertures C and D ( fig. 6). Aperture A covers the distinct source in the eastern part of the bubble (source "A"), while aperture B covers the faint emission in its western (right) part (source "B"). Aperture C covers part of the fainter emission outside the bubble (source "C") and aperture D covers the distinct radio source located eastern of the bubble (source "D"). We subtracted a background derived from flux within an aperture covering a relatively homogeneous region free of visible radio sources ( fig. 6). Then we used the two fluxes at 1420 and 327 MHz to derive the spectral index α for both sources A and B: where S 1 , S 2 are the fluxes at frequencies f 1 , f 2 . Note that the theoretical spectral index for a classical H ii region in the optically thin limit is −0.1 (Condon & Ransom 2010), although the usually observed values are around 0 or more (Rohlfs & Wilson 1996;Kellermann & Verschuur 1988). In distinction to H ii regions, the mean observed spectral index of supernova remnants is −0.5 (average value in catalogue of Green 2009). The measured flux densities and the derived spectral indices are presented in table 1. The spectral index of source A of −0.3 suggests a nonthermal (synchrotron) contribution to the received radio flux: a supernova No known supernova remnant is associated, yet, with that radio source. For source B, the spectral index of ≈ 0 suggests a thermal origin of the flux: a classical H ii region. The spectral indices of sources C (−0.2) and D (−0.1) suggest that part of their emission is also nonthermal.

Analysis of Molecular Clumps
We decomposed the molecular ring associated with N107 into individual molecular clumps with the program DENDROFIND, described in the appendix of Wünsch et al. (2012). DEN-DROFIND processes data cubes -typically position-positionvelocity data cubes of the brightness temperature T b -and outputs a list of molecular clumps organised in a hierarchical structure, which can be represented by a dendrogram (Rosolowsky et al. 2008). In our case, the data cube is a cutout of the 13 CO GRS data and has two spatial dimensions (galactic longitude, galactic latitude) and one radial velocity dimension. Each clump consists of a dense core surrounded by a less dense molecular envelope. The core is observed as a local T b peak. DEN-DROFIND is controlled by four parameters: Nlevels sets the number of steps by which it goes through the data cube when searching for clump peaks; Npxmin sets the minimal number of pixels which can form a clump; Tcutoff sets the minimal T b which is yet considered a signal, pixels with lower T b are discarded; dTleaf sets the minimal height of a local peak to be considered a clump's core; a higher dTleaf will merge more clumps. Parameters Tcutoff and dTleaf are designed to be set to a value of ≈ 3 × the noise level. We ran DENDROFIND with the following parameters: Nlevels = 1000, Npixmin = 5, Tcutoff = dTleaf = 1.0 ≈ 3σ noise , where σ noise ≈ 0.34 K is the estimated standard deviation of the gaussian noise found within our N107 cutout of the 13 CO GRS data. Note that this value is slightly greater then the 'typical' value of the entire GRS of 0.27 K given in the paper of Jackson et al. (2006). In order to select only the clumps associated with bubble N107, we took only those clumps of which cores lie within the aperture 2 used for measuring the total molecular mass (sec. 2.4). Table 2 lists the identified molecular clumps. Note that only the first 10 clumps with the highest core T b are shown there; a complete list is available in an online appendix (tab. A.1). The meaning of the columns is as follows: N is the number identifying a clump in our list; l peak , b peak , v peak are the galactic longitude, galactic latitude and LSR radial velocity, respectively, of the brightest pixel of the clump. T peak is the brightness temperature of the brightest pixel. R is the clump equivalent radius, defined as where A cl is the clump projected area in the sky. FWHM is the clump full width at half maximum computed as where x ∈ {l, b, v}. The second term under the big square root represents the pixel resolution term and ensures that the FWHM is not smaller than the pixel size. For the GRS data cube, the pixel size ∆x is 22 . 1 × 22 . 1 × 0.21 km/s, which corresponds at N107's distance of 3600 pc to 0.39 pc × 0.39 pc × 0.21 km/s. If ∆x goes to zero, the FWHM approaches the gaussian FWHM. σ 2 x is the variance weighted by the brightness temperature T b : where with the summation taken over all pixels belonging to the clump. M is the clump mass, which is the total mass of pixels belonging to the clump, derived from equations 3-6 (sec. 2.4). The error of the clump mass, assuming gaussian noise and neglecting pixel correlation within the beam, is where N px is the number of pixels of the clump and ∆M px is the σ noise equivalent mass. In our case of the GRS 13 CO data cubes, σ noise = 0.34 K ⇒ ∆M px = 0.21 M . M vir is the clump virial mass, estimated as (MacLaren et al. 1988): The formula above assumes a spherical clump with a density profile ∝ 1/r. ∆M vir is the virial mass error derived from the first order Taylor expansion of eq. 13, assuming the clump's area and FWHM v errors of one pixel. Fig. 8 shows a histogram of clumps masses. Overlaid is a line representing a fit by a function dN/dM = C(M/M ) α , where α is the slope of the clump mass function and C is some constant. We used least-square fitting on a logarithmic scale, i.e. we fit the log of the bin heights log(dN/dM) by a function log(C(M/M ) α ). The best fitting slope is α = −1.1. The bins in the histogram have variable widths, set so that each contains 7 clumps. This ensures that all bins have the same significance.

Numerical Simulations of Bubbles and Comparison with N107
We simulated the evolution of bubbles formed in the ISM around massive stars in order to estimate parameters of N107, which cannot be derived directly from observation (esp. its age). We used a version of the FORTRAN code ring (Palouš 1990;Ehlerová & Palouš 1996), which was originally written for simulations of expanding H i shells in the Galactic ISM.   Inside the bubble, the radio continuum is well correlated with the 24 µm emission, except for the eastern part, where the 24 µm emission is absent and only the radio continuum is observed.

Code Description
The simulation starts with a small spherical shell, which is set according to the wind-blown bubble model of Weaver et al. (1977) in the pressure-driven snowplough phase. The infinitesimally thin shell (thin-shell approximation) is divided into nl × np elements. We solve the equations of motion including the pressure difference between the shell interior and exterior and the mass accumulation of the ambient medium. We do not include gravity or any other force acting on the shell. The momentum equation for an element is  where m is the element mass, v is the element velocity, S is the element surface and P i and P e are the interior and exterior pressures. The total number of elements is fixed (we used 30 × 60), so their surface and mass grow with time. The mass growth is given by the equation where ρ is the local density of the ambient medium. For the sake of simplicity, an energy source is represented as a continuous and constant energy inputĖ. For each timestep a corresponding amount of energy is added to the bubble interior and a new value of the interior pressure is computed as where E i is the interior thermal energy and V is the bubble volume. Eq. 16 assumes that the interior is filled with monoatomic gas expanding adiabatically. In the case of a shell expanding into a uniform medium, this type of energy source gives for the shell radius, as a function of time, a power-law with an index 3/5. This is in between the power-law index 2/5, given by a single (supernova) explosion and the power-law index 4/7, valid for an expanding H ii region. The motion of the bubble elements is integrated using the fourth-order Runge-Kutta method with adaptive timestep.

Simulation Setup
We assume a bubble expanding in a molecular cloud with a 3D gaussian density profile (see fig. 9): where n is the particle density, n 0 is the particle density in the centre of the cloud; σ x , σ y , σ z are the gaussian thicknesses (standard deviation) in the x, y, z direction; x 0 , y 0 , z 0 specify the location of the cloud centre. We also assume the centre of expansion is shifted/dislocated from the centre of the cloud by ∆r 0 towards direction (−x, −y). In the code, the expansion centre (bubble centre) is located at x = y = z = 0, so we shift the cloud centre by ∆r 0 in the xy plane: We assume a hydrogen abundance of X = 0.735, helium abundance of Y = 0.245, metallicity of Z = 0.02, and a mean relative particle weight of metals of 16.65. We also assume that all hydrogen in the cloud is in the form of H 2 . This yields a mean relative particle weight of the molecular matter of µ = 2.34. The molecular material of the cloud is isothermal with the temperature set at 20 K. We also consider the atomic gas component, which is represented by a density background with the constant value of 4 hydrogen atoms per cm 3 and a temperature of 6000 K. The total initial energy of the bubble is set to 1 × 10 49 erg and each time step δt the energy is increased byĖ δt.

Comparison with Observation
We ran a series of simulations in order to find the parameters which produce a bubble with similar mass distribution and expansion velocity as N107. The varying input parameters for ring were the energy input (Ė), the shift/dislocation of the expansion centre from the cloud centre (∆r 0 ), the gaussian thickness of the cloud in the direction of our line of sight (σ z ) and the total mass Article number, page 9 of 15  Fig. 9. Simulation setup for a stellar-blown bubble. The bubble is formed around a massive star(s) inside a molecular cloud, with a 3D gaussian density profile. The gaussian thickness σ x and σ y , in the sky plane, is set to 12 pc, while σ z , which runs parallel to out LOS, is a parameter we vary. The centre of expansion is dislocated from the centre of the cloud by ∆r 0 , which is another parameter we vary. The dislocation is towards the (−x, −y) direction, perpendicular to the LOS. The molecular cloud is overlaid over a background of atomic gas with the particle density of 4 cm −3 . of the molecular cloud (M 0 ). The gaussian thickness in the sky plane σ xy = σ x = σ y was kept constant at 12 pc, which is comparable to the radius of N107. For each set of these parameters, the code evolves a bubble in time and outputs a chain of its snapshots in given time intervals.
Because the bubble is expanding in an inhomogeneous, nonspherical density field, the elements facing a lower density expand faster. This causes the bubble's geometrical centre to move progressively towards the direction of a lower density. To compensate for this drift, we centred each snapshot at the bubble's geometrical centre, which is computed as the mean (x, y, z) coordinates of the bubble elements weighted by their surfaces. We projected these centred snapshots onto FITS data cubes with two spatial (XY) and one radial velocity (RV) dimensions with the same pixel size as in the GRS data cube.
The radiation from OB stars illuminating molecular clouds creates a layer of dissociated matter. This layer then serves as a shield against the radiation and prevents it from penetrating further inside. The dissociated matter will not contribute to the 13 CO line brightness temperature. To account for that, before computing the 13 CO line brightness temperature, we decreased the surface density of each bubble element by 4.35 M pc −2 . This value is based on the suggestion for shielding by Krumholz et al. (2009). We assume that the 13 CO abundances are the same as we used for measuring the observed molecular mass in sec. 2.4. Then we derived for each pixel the corresponding brightness temperature of the 13 CO (J = 1−0) line. Finally, we applied convolution to simulate the FCRAO telescope resolution used in the GRS. The resulting FITS thus represents an image, which would be seen by FCRAO if it observed the simulated bubble at a given evolution time.
We define a goodness-of-fit parameter δ, to compare the synthetic data cubes to the observations: . Ranges of the input parametersĖ, r 0 , σ z , M 0 and the constraint of the evolution time (t) we used for the numerical simulations of expanding bubbles.

Range of values Steps
First series ( Lower values of δ mean a better fit to (better agreement with) the observed properties of N107. The meaning of the terms is as follows: represents the square of the difference of the bubble expansion velocity between the observed value of N107 (V exp, o ) and simulated value (V exp, s ). We adopt the observed bubble expansion velocity of V exp, o = 8 km/s, which is derived from the H i observations and represents half of the relative velocity of the shell's front and back wall.
is a sum taken over four circular sectors of the aperture 2 used for measuring the bubble molecular mass (see sec. 2.4). This term compares the angular mass distribution between N107 and the simulated bubbles. A visualisation of the sectors is show in fig. 10. m i, s and m i, o are the simulated ("s") and observed ("o") molecular masses in circular sector i. r i, s and r i, o give the simulated and observed mean radii in sector i, weighted by the brightness temperature: where the sums are taken over the pixels in sector i. First, we ran a series of simulations with a wider span of the input parametersĖ, r 0 , σ z , M 0 with steps on a logarithmic scale and the evolution time up to t = 100 Myr with snapshots taken every 0.5 Myr. Then we located a subset of the input parameters with best fits to the observation (lowest δ) and ran a second series of simulations with finer parameter resolutions. In the second series, we used linear steps for parameters r 0 and σ z and logarithmic steps forĖ and M 0 . The snapshots were taken every 0.25 Myr up to t = 8 Myr. Note that the internal time step of ring Table 4. First ten sets of parameters best fitting the observation of N107. m 1,2,3,4 are the masses of four sectors of the bubble. m 1 /m 3 is the ratio of masses from the first and third sector, which quantify the mass contrast between the most and least massive sector. δ is a goodness-of-fit parameter (lower is better). The simulations are divided into two groups A and B, discussed at the end of sec. 4.3. Further explanation of the table columns is also given in sec. 4.3. The asterisk (*) marks two simulations -best of their group -which are visualised in fig. 10 To conclude, we simulated the evolution of bubbles formed around massive stars in molecular clouds and compared their properties to N107. We got the best agreement for two groups of parameters: 1. Group A represents bubbles produced with a lower energy input (Ė < 0.4 × 10 50 erg/Myr), greater centre-of-expansion dislocation (r 0 ≈ 16 . . . 22 pc), formed within an oblate molecular cloud with the oblateness of (σ z /σ xy = 6/12), and with a greater evolution time (t ≈ 1.75 . . . 2.25 Myr). 2. Group B represents bubbles produced with a higher energy input (Ė = 2.7 × 10 50 erg/Myr), smaller centre-of-expansion dislocation (r 0 ≈ 10 . . . 12 pc), formed within a prolate molecular cloud (σ z /σ xy ≈ (20 . . . 28)/12), with a smaller evolution time (t = 1 Myr).
Fits show little sensitivity to the parameter M 0 -the initial mass of the surrounding molecular cloud (M 0 ≈ 5.8 . . . 15.5×10 5 M ). Both groups of the best-fittings-simulations suggest that bubble N107 is relatively young (t 2.25 Myr) and was formed at the very edge of its parental molecular cloud. Fig. 10 shows snapshots from the best fitting simulations from both groups A and B and the GRS observation of the CO line.

Numerical Simulation Setup, Assumptions and Validity of Approximations
An important assumption for the simulation setup is the shape and density profile of the parental molecular cloud. In our simulations, we chose the gaussian density profile and tested different levels of the oblateness (ratio of σ z to σ xy ). The oblate shape of the parental molecular clouds which host bubbles similar to N107 is supported by a study of Beaumont & Williams (2010), who concluded that the Spitzer bubbles are formed in oblate, sheet-like, molecular clouds with a thickness of a few parsecs. Our simulations, however, favour neither oblate nor prolate clouds. The numerical code ring assumes the thin-shell approximation, which is applicable so long as the bubble expansion is supersonic, i.e. expanding with a velocity higher than the sound speed of the ambient cold gas: ≈ 0.3 km/s (for T = 20 K, γ = 7/5, µ = 2.34). In our runs, almost all the bubble's elements move faster through the ambient gas during the whole bubble evolution, so the thin-shell approximation is reasonable. Moreover, if the speed of an element decreases below the sound speed, the code stops the mass accumulation for that element.

Source of Radio Continuum
The radio continuum emission coming from the direction of the bubble interior has both thermal and nonthermal components. While the strong radio source A is dominated by nonthermal radiation, fainter source B is emitting thermally (fig. 6). The nonthermal source A is brighter than the thermal source B by almost one order of magnitude, so we can expect that the thermal emission is present in the whole bubble volume, only outshone in aperture A by the nonthermal source. Mixing of thermal and nonthermal radiation in aperture A is also suggested by the visual comparison of the radio continuum and the 24 µm emission ( fig. 7). The northern part of aperture A features both the radio continuum and the 24 µm emission, while in the southern part only the radio continuum is observed. Further, merging of thermal and nonthermal radiation in aperture A explains why the measured spectral index (−0.3) is closer to 0 than the typical value for supernova remnants (−0.5, the average value in the SNRs catalogue of Green 2009).
The negative spectral indices of sources C (−0.2) and D (−0.1) suggest that part of their radio continuum is also emitted by a nonthermal source. The nonthermal contribution to the flux is stronger for source C than for source D. This is expected, since aperture D covers a distinct 24 µm (thermal) source, while there is almost no 24 µm emission in aperture C ( fig. 7).
Given that first supernovae in an OB association explode only after 3 Myr or later, and considering the bubble evolution time estimated from simulations (< 2.25 Myr), no supernova remnant should be present inside the bubble, yet. In other words,  our analysis suggests that a supernova remnant, not associated with the bubble, may be present in the direction towards the bubble.

Missing OB Association
We searched the SIMBAD database and a catalogue of Galactic OB stars by Reed (2003) for sources which could be responsible for the bubble creation. However, no OB star, supernova remnant or pulsar is known to lie in the direction of the bubble interior or in its vicinity.
We also performed a search for massive stars possibly related to N107 using the UKIDSS catalogue 7th data release (Lawrence et al. 2007), stellar evolutionary tracks by Bressan et al. (2012) and the standard model of extinction (Rieke & Lebofsky 1985). We derived the extinction for the massive star candidates in the field of N107 and compared it to the values of extinction given for that direction by Stead & Hoare (2009). The search for massive stars involves many uncertainties, but the comparison of the derived extinctions with independent values is a good test of the credibility of our results. We ended up with six massive star candidates, which lie at distances comparable with the distance of N107 (3.6 kpc). None of these six stars have measured radial velocities. Three of them lie at the rim of the bubble, so they can hardly be its progenitors, but the others are possible energy sources. The selected massive stars are shown in fig. 11 and listed in table 5.

Molecular Clump Mass Spectrum
The best fitting slope for the mass spectrum of molecular clumps found along the borders of N107 is − 1.1 (fig. 8). This is a significantly shallower slope than that of the classical stellar initial mass function (Salpeter: −2.35) and also shallower than the clump mass function slope of −1.6 to −1.8 derived by Kramer et al. (1998). Note that Kramer et al. (1998) used gaussian decomposition for clump identification, which is a fundamentally different method than we used and, moreover, studied different types of clouds -not shell-like as is our case of N107.
The least massive clumps we found have masses about 3.5 M . This limit is given by the angular resolution and sensitivity of the observations of the 13 CO line. With better resolution and sensitivity, we expect that some clumps, especially the most massive ones, would be identified as several less massive clumps. This would alter the shape of the clump mass spectrum, and likely redistribute the mass of the largest clumps into less massive ones, so that the slope would be steeper. Table 5. Possible stellar progenitors of bubble N107 found in the UKIDSS catalogue. First column gives the source label as marked in fig. 11. Second column gives the ID from the UKIDSS catalogue. Third and fourth columns give the galactic longitude and latitude. Fifth column gives the estimated mass of the source and sixth column gives the estimated distance. label in fig. 11

Conclusion
We studied dust bubble N107, discovered by Churchwell et al. (2006) with infrared observations by the Spitzer Space Telescope. Radio observations of the H i and 13 CO (J = 1−0) lines revealed its atomic and molecular components, which allowed us to determine the bubble radial velocity, kinematical distance and mass of these components. Using the code DENDROFIND (Wünsch et al. 2012) we decomposed the molecular ring associated with the bubble into 49 individual clumps and found that the clump mass spectrum slope is about −1.1. Radio continuum observations at 1420 and 327 MHz suggest the radio flux coming from the bubble direction has both thermal and nonthermal components; i.e., besides a classical H ii region, a supernova remnant contributes to the emission. We simulated the evolution of stellar-blown bubbles within molecular clouds with the numerical code ring (Palouš 1990;Ehlerová & Palouš 1996). We were able to produce two groups of bubbles with a similar angular mass distribution and expansion velocity as N107: (A) Bubbles formed within an oblate molecular cloud (σ z /σ xy = 6/12), with the energy inputĖ < 0.4 × 10 50 erg/Myr and the evolution time t ≈ 1.75 . . . 2.25 Myr. (B) Bubbles formed within a prolate molecular cloud (σ z /σ xy = (20 . . . 28)/12), with the energy in-putĖ ≈ 2.7 × 10 50 erg/Myr and the evolution time t = 1 Myr. Considering the estimated bubble evolution time (< 2.25 Myr), no supernova remnant should be present inside the bubble, yet. This may be explained by a supernova remnant present in the direction towards the bubble, however not associated with it. A summary of the parameters we derived for bubble N107, together with the parameters given by Churchwell et al. (2006), is presented in table 6.