Starbursts triggered by cloud compression in interacting galaxies

We propose a physical mechanism for the triggering of starbursts in interacting spiral galaxies by shock compression of the pre-existing disk giant molecular clouds (GMCs). We show that as a disk GMC tumbles into the central region of a galaxy following a galactic tidal encounter, it undergoes a radiative shock compression by the pre-existing high pressure of the central molecular intercloud medium. The shocked outer shell of a GMC becomes gravitationally unstable, which results in a burst of star formation in the initially stable GMC. In the case of colliding galaxies with physical overlap such as Arp 244, the cloud compression is shown to occur due to the hot, high-pressure remnant gas resulting from the collisions of atomic hydrogen gas clouds from the two galaxies. The resulting values of infrared luminosity agree with observations. The main mode of triggered star formation is via clusters of stars, thus we can naturally explain the formation of young, luminous star clusters observed in starburst galaxies.

Thermal instability and the feedback regulation of hot haloes in clusters, groups and galaxies

We present global multidimensional numerical simulations of the plasma that pervades the dark matter haloes of clusters, groups and massive galaxies (the intracluster medium; ICM). Observations of clusters and groups imply that such haloes are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally motivated heating prescriptions, we show that local thermal instability (TI) can produce a multiphase medium with similar to 104 K cold filaments condensing out of the hot ICM only when the ratio of the TI time-scale in the hot plasma (tTI) to the free-fall time-scale (tff) satisfies tTI/tff? 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of similar to 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling flow is similar to 10-3 for clusters and decreases for lower mass haloes; supernova heating may be energetically sufficient to balance cooling in galactic haloes. We further argue that the ICM self-adjusts so that tTI/tff? 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives tTI/tff? 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminositytemperature () relation. The high-velocity clouds observed in the Galactic halo can be due to local TI producing multiphase gas close to the virial radius if the density of the hot plasma in the Galactic halo is >rsim 10-5 cm-3 at large radii.

The intergalactic and circumgalactic medium surrounding star-forming galaxies at redshifts 2 < z < 3

We present measurements of the spatial distribution, kinematics, and physical properties of gas in the circumgalactic medium (CGM) of 2.0<z<2.8 UV color-selected galaxies as well as within the 2<z<3 intergalactic medium (IGM). These measurements are derived from Voigt profile decomposition of the full Lyα and Lyβ forest in 15 high-resolution, high signal-to-noise ratio QSO spectra resulting in a catalog of ∼6000 HI absorbers.

Chapter 2 of this thesis focuses on HI surrounding high-z star-forming galaxies drawn from the Keck Baryonic Structure Survey (KBSS). The KBSS is a unique spectroscopic survey of the distant universe designed to explore the details of the connection between galaxies and intergalactic baryons within the same survey volumes. The KBSS combines high-quality background QSO spectroscopy with large densely-sampled galaxy redshift surveys to probe the CGM at scales of ∼50 kpc to a few Mpc. Based on these data, Chapter 2 presents the first quantitative measurements of the distribution, column density, kinematics, and absorber line widths of neutral hydrogen surrounding high-z star-forming galaxies.

Chapter 3 focuses on the thermal properties of the diffuse IGM. This analysis relies on measurements of the ∼6000 absorber line widths to constrain the thermal and turbulent velocities of absorbing "clouds." A positive correlation between the column density of HI and the minimum line width is recovered and implies a temperature-density relation within the low-density IGM for which higher-density regions are hotter, as is predicted by simple theoretical arguments.

Chapter 4 presents new measurements of the opacity of the IGM and CGM to hydrogen-ionizing photons. The chapter begins with a revised measurement of the HI column density distribution based on this new absorption line catalog that, due to the inclusion of high-order Lyman lines, provides the first statistically robust measurement of the frequency of absorbers with HI column densities 14 ≲ log(N_HI/cm^-2) ≲ 17.2. Also presented are the first measurements of the column density distribution of HI within the CGM (50 <d < 300 pkpc) of high-z galaxies. These distributions are used to calculate the total opacity of the IGM and IGM+CGM and to revise previous measurements of the mean free path of hydrogen-ionizing photons within the IGM. This chapter also considers the effect of the surrounding CGM on the transmission of ionizing photons out of the sites of active star-formation and into the IGM.

This thesis concludes with a brief discussion of work in progress focused on understanding the distribution of metals within the CGM of KBSS galaxies. Appendix B discusses my contributions to the MOSFIRE instrumentation project.

Aerosol-cloud-precipitation interactions in marine stratocumulus clouds

Marine stratocumulus clouds are generally optically thick and shallow, exerting a net cooling influence on climate. Changes in atmospheric aerosol levels alter cloud microphysics (e.g., droplet size) and cloud macrophysics (e.g., liquid water path, cloud thickness), thereby affecting cloud albedo and Earth’s radiative balance. To understand the aerosol-cloud-precipitation interactions and to explore the dynamical effects, three-dimensional large-eddy simulations (LES) with detailed bin-resolved microphysics are performed to explore the diurnal variation of marine stratocumulus clouds under different aerosol levels and environmental conditions. It is shown that the marine stratocumulus cloud albedo is sensitive to aerosol perturbation under clean background conditions, and to environmental conditions such as large-scale divergence rate and free tropospheric humidity.

Based on the in-situ Eastern Pacific Emitted Aerosol Cloud Experiment (E-PEACE) during Jul. and Aug. 2011, and A-Train satellite observation of 589 individual ship tracks during Jun. 2006-Dec. 2009, an analysis of cloud albedo responses in ship tracks is presented. It is found that the albedo response in ship tracks depends on the mesoscale cloud structure, the free tropospheric humidity, and cloud top height. Under closed cell structure (i.e., cloud cells ringed by a perimeter of clear air), with sufficiently dry air above cloud tops and/or higher cloud top heights, the cloud albedo can become lower in ship tracks. Based on the satellite data, nearly 25% of ship tracks exhibited a decreased albedo. The cloud macrophysical responses are crucial in determining both the strength and the sign of the cloud albedo response to aerosols.

To understand the aerosol indirect effects on global marine warm clouds, multisensory satellite observations, including CloudSat, MODIS, CALIPSO, AMSR-E, ECMWF, CERES, and NCEP, have been applied to study the sensitivity of cloud properties to aerosol levels and to large scale environmental conditions. With an estimate of anthropogenic aerosol fraction, the global aerosol indirect radiative forcing has been assessed.

As the coupling among aerosol, cloud, precipitation, and meteorological conditions in the marine boundary layer is complex, the integration of LES modeling, in-situ aircraft measurements, and global multisensory satellite data analyses improves our understanding of this complex system.

Deep near-infrared spectroscopy of high-redshift galaxies: the physical growth of passive systems

The assembly history of massive galaxies is one of the most important aspects of galaxy formation and evolution. Although we have a broad idea of what physical processes govern the early phases of galaxy evolution, there are still many open questions. In this thesis I demonstrate the crucial role that spectroscopy can play in a physical understanding of galaxy evolution. I present deep near-infrared spectroscopy for a sample of high-redshift galaxies, from which I derive important physical properties and their evolution with cosmic time. I take advantage of the recent arrival of efficient near-infrared detectors to target the rest-frame optical spectra of z > 1 galaxies, from which many physical quantities can be derived. After illustrating the applications of near-infrared deep spectroscopy with a study of star-forming galaxies, I focus on the evolution of massive quiescent systems.

Most of this thesis is based on two samples collected at the W. M. Keck Observatory that represent a significant step forward in the spectroscopic study of z > 1 quiescent galaxies. All previous spectroscopic samples at this redshift were either limited to a few objects, or much shallower in terms of depth. Our first sample is composed of 56 quiescent galaxies at 1 < z < 1.6 collected using the upgraded red arm of the Low Resolution Imaging Spectrometer (LRIS). The second consists of 24 deep spectra of 1.5 < z < 2.5 quiescent objects observed with the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE). Together, these spectra span the critical epoch 1 < z < 2.5, where most of the red sequence is formed, and where the sizes of quiescent systems are observed to increase significantly.

We measure stellar velocity dispersions and dynamical masses for the largest number of z > 1 quiescent galaxies to date. By assuming that the velocity dispersion of a massive galaxy does not change throughout its lifetime, as suggested by theoretical studies, we match galaxies in the local universe with their high-redshift progenitors. This allows us to derive the physical growth in mass and size experienced by individual systems, which represents a substantial advance over photometric inferences based on the overall galaxy population. We find a significant physical growth among quiescent galaxies over 0 < z < 2.5 and, by comparing the slope of growth in the mass-size plane dlogR_e/dlogM_∗ with the results of numerical simulations, we can constrain the physical process responsible for the evolution. Our results show that the slope of growth becomes steeper at higher redshifts, yet is broadly consistent with minor mergers being the main process by which individual objects evolve in mass and size.

By fitting stellar population models to the observed spectroscopy and photometry we derive reliable ages and other stellar population properties. We show that the addition of the spectroscopic data helps break the degeneracy between age and dust extinction, and yields significantly more robust results compared to fitting models to the photometry alone. We detect a clear relation between size and age, where larger galaxies are younger. Therefore, over time the average size of the quiescent population will increase because of the contribution of large galaxies recently arrived to the red sequence. This effect, called progenitor bias, is different from the physical size growth discussed above, but represents another contribution to the observed difference between the typical sizes of low- and high-redshift quiescent galaxies. By reconstructing the evolution of the red sequence starting at z ∼ 1.25 and using our stellar population histories to infer the past behavior to z ∼ 2, we demonstrate that progenitor bias accounts for only half of the observed growth of the population. The remaining size evolution must be due to physical growth of individual systems, in agreement with our dynamical study.

Finally, we use the stellar population properties to explore the earliest periods which led to the formation of massive quiescent galaxies. We find tentative evidence for two channels of star formation quenching, which suggests the existence of two independent physical mechanisms. We also detect a mass downsizing, where more massive galaxies form at higher redshift, and then evolve passively. By analyzing in depth the star formation history of the brightest object at z > 2 in our sample, we are able to put constraints on the quenching timescale and on the properties of its progenitor.

A consistent picture emerges from our analyses: massive galaxies form at very early epochs, are quenched on short timescales, and then evolve passively. The evolution is passive in the sense that no new stars are formed, but significant mass and size growth is achieved by accreting smaller, gas-poor systems. At the same time the population of quiescent galaxies grows in number due to the quenching of larger star-forming galaxies. This picture is in agreement with other observational studies, such as measurements of the merger rate and analyses of galaxy evolution at fixed number density.

Galactic Edge Clouds I. Molecular Line Observations and Chemical Modelling of Edge Cloud 2

Edge Cloud 2 (EC2) is a molecular cloud, about 35 pc in size, with one of the largest galactocentric distances known to exist in the Milky Way. We present observations of a peak CO emission region in the cloud and use these to determine its physical characteristics. We calculate a gas temperature of 20 K and a density of n(H2)~10^4 cm-3. Based on our CO maps, we estimate the mass of EC2 at around 10^4 Msolar and continuum observations suggest a dust-to-gas mass ratio as low as 0.001. Chemical models have been developed to reproduce the abundances in EC2, and they indicate that heavy element abundances may be reduced by a factor of 5 relative to the solar neighborhood (similar to dwarf irregular galaxies and damped Lya systems), very low extinction (A_V <4 mag) due to a very low dust-to-gas mass ratio, an enhanced cosmic-ray ionization rate, and a higher UV field compared to local interstellar values. The reduced abundances may be attributed to the low level of star formation in this region and are probably also related to the continuing infall of primordial (or low-metallicity) halo gas since the Milky Way formed. Finally, we note that shocks from the old supernova remnant GSH 138-01-94 may have determined the morphology and dynamics of EC2.

Observations towards early-type stars in the ESO-POP Survey - II. Searches for intermediate- and high-velocity clouds

We present Ca it K and Ti it optical spectra of early-type stars taken mainly from the ultraviolet and visual echelle spectrograph (LIVES) Paranal Observatory Project, plus H 1 21-cm spectra, from the Vila-Elisa and Leiden-Dwingeloo Surveys, which are employed to obtain distances to intermediate- and high-velocity clouds (IHVCs). H I emission at a velocity of -117 km s(-1) towards the sightline HD 30677 (l, b = 190 degrees.2, -22 degrees.2) with column density -1.7 x 10(19) cm(-2) has no corresponding Ca Pi K absorption in the LIVES spectrum, which has a signal-to-noise ratio (S/N) of 610 per resolution element. The star has a spectroscopically determined distance of 2.7 kpc, and hence sets this as a firm lower distance limit towards Anti-Centre cloud ACII. Towards another sightline (HD 46185 with 1, b = 222 0, -10 degrees.1), H1 at a velocity of +122 km s(-1) and column density of 1.2 x 10(19) cm(-2) is seen. The corresponding Ca Pi K spectrum has a S/N of 780, although no absorption is observed at the cloud velocity. This similarly places a firm lower distance limit of 2.9 kpc towards this parcel of gas that may be an intermediate-velocity (IV) cloud. The lack of IV Ca it absorption towards HD 196426 (1, b = 45 degrees.8, -23 degrees.3) at a S/N of 500 reinforces a lower distance limit of -700 pc towards this part of complex gp, where the H I column density is 1.1 x 1019 cm(-2) and velocity is +78 km s(-1). Additionally, no IV Cart is seen in absorption in the spectrum of HD 19445, which is strong in H I with a column density of 8 x 10(19) cm(-2) at a velocity of - -42 km s(-1), placing a firm although uninteresting lower distance limit of 39 pc to this part of IV South. Finally, no high-velocity Call K absorption is seen towards HD 115363 (l, b = 306.0,-1.0) at a S/N of 410, placing a lower distance of -3.2 kpc towards the HVC gas at velocity of - +224 km s(-1) and WE column density of 5.2 x 10(19) cm(-2). This gas is in the same region of the sky as complex WE (Wakker 2001), but at higher velocities. The non-detection of Ca it K absorption sets a lower distance of -3.2 kpc towards the HVC, which is unsurprising if this feature is indeed related to the Magellanic System.

DENSE GAS IN NEARBY GALAXIES .9. A SURVEY FOR HNC

Ten detections and five tentative detections of hydrogen isocyanide (HNC) J=1-0 emission are reported from a survey including sixteen galaxies. Full maps are presented for the nuclear regions of NGC 253 and IC 342, partial maps for Maffei 2, M 82, and M 83. Toward IC 342, the HNC and HCO+ distributions differ from those observed in 12CO, 13CO, HCN, CS, and NH3. This is likely a consequence of the density structure. Relative HNC abundances are with 10(-10)-10(-9) much smaller than those measured in nearby dark clouds and appear to be slightly smaller than those in regions of massive star formation of the Galactic disk. This is consistent with the presence of dense warm gas or a frequent occurrence of shocks in the nuclear regions of the galaxies observed. As in prominent Galactic star forming regions, 3 mm HNC line emission tends to be weaker than the corresponding emission from HCN and HCO+. Toward Arp 220, however, the 3 mm HNC/HCN line intensity ratio is > 1. HNC/HCO+, HNC/CO, and HNC to 20 cm radio continuum luminosity ratios are also particularly large. A possible interpretation is the presence of cool quiescent gas outside the central region which contains the starburst. In the other ultraluminous galaxy observed, NGC 6240, X(HNC) 10 smaller than in Arp 220, demonstrating that the molecular composition in ultraluminous galaxies is far from being uniform.

Distance limits to intermediate- and high-velocity clouds

We present optical spectra of 403 stars and quasi-stellar objects in order to obtain distance limits towards intermediate- and high-velocity clouds (IHVCs), including new Fibre-fed Extended Range Optical Spectrograph (FEROS) observations plus archival ELODIE, FEROS, High Resolution Echelle Spectrometer (HIRES) and Ultraviolet and Visual Echelle Spectrograph (UVES) data. The non-detection of Ca II K interstellar (IS) absorption at a velocity of −130 to −60 km s−1 towards HDE 248894 (d ∼ 3 kpc) and HDE 256725 (d ∼ 8 kpc) in data at signal-to-noise ratio (S/N) > 450 provides a new firm lower distance limit of 8 kpc for the anti-centre shell HVC. Similarly, the non-detection of Ca II K IS absorption towards HD 86248 at S/N ∼ 500 places a lower distance limit of 7.6 kpc for Complex EP, unsurprising since this feature is probably related to the Magellanic System. The lack of detection of Na I D at S/N = 35 towards Mrk 595 puts an improved upper limit for the Na I column density of log (NNaD <) 10.95 cm−2 towards this part of the Cohen Stream where Ca II was detected by Wakker et al. Absorption at ∼ −40 km s−1 is detected in Na I D towards the Galactic star PG 0039+049 at S/N = 75, placing a firm upper distance limit of 1 kpc for the intermediate-velocity cloud south (IVS), where a tentative detection had previously been obtained by Centurion et al. Ca ´ II K and Na I D absorption is detected at −53 km s−1 towards HD 93521, which confirms the upper distance limit of 2.4 kpc for part of the IV arch complex obtained using the International Ultraviolet Explorer (IUE) data by Danly. Towards HD 216411 in Complex H a non-detection in Na D towards gas with log(NH I) = 20.69 cm−2 puts a lower distance limit of 6.6 kpc towards this HVC complex. Additionally, Na I D absorption is detected at −43.7 km s−1 in the star HD 218915 at a distance of 5.0 kpc in gas in the same region of the sky as Complex H. Finally, the Na I/Ca II and Ca II/H I ratios of the current sample are found to lie in the range observed for previous studies of IHVCs.

Blue luminous stars in nearby galaxies: Quantitative spectral analysis of M33 B-type supergiant stars

We present the detailed spectral analysis of a sample of M33 B-type supergiant stars, aimed at the determination of their fundamental parameters and chemical composition. The analysis is based on a grid of non-LTE metal line-blanketed model atmospheres including the effects of stellar winds and spherical extension computed with the code FASTWIND. Surface abundance ratios of C, N, and O are used to discuss the chemical evolutionary status of each individual star. The comparison of observed stellar properties with theoretical predictions of massive star evolutionary models shows good agreement within the uncertainties of the analysis. The spatial distribution of the sample allows us to investigate the existence of radial abundance gradients in the disk of M33. The comparison of stellar and H II region O abundances ( based on direct determinations of the electron temperature of the nebulae) shows good agreement. Using a simple linear radial representation, the stellar oxygen abundances result in a gradient of -0.0145 +/- 0.005 dex arcmin(-1) (or -0.06 +/- 0.02 dex kpc(-1)) up to a distance equal to similar to 1.1 times the isophotal radius of the galaxy. A more complex representation cannot be completely discarded by our stellar sample. The stellar Mg and Si abundances follow the trend displayed by O abundances, although with shallower gradients. These differences in gradient slope cannot be explained at this point. The derived abundances of the three alpha-elements yield solar metallicity in the central regions of the disk of M33. A comparison with recent planetary nebula data from Magrini and coworkers indicates that the disk of M33 has not suffered from a significant O enrichment in the last 3 Gyr.

A Stubbornly Large Mass of Cold Dust in the Ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 mu m data and improved imaging quality at 100 and 160 mu m compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 mu m [O-I] line flux, eliminating the possibility that line contaminations distort the previously estimated dustmass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 mu m flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5 +/- 0.1M(circle dot) of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3M(circle dot) of amorphous carbon and 0.5M(circle dot) of silicates, totalling 0.8M(circle dot) of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

Expanded iron uta spectra - Probing the thermal stability limits in AGN clouds

The Fe unresolved transition arrays (UTAs) produce prominent features in the 15-17 Å wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in AGNs. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration. © 2013. The American Astronomical Society. All rights reserved.

Selecting superluminous supernovae in faint galaxies from the first year of the Pan-STARRS1 Medium Deep Survey

The Pan-STARRS1 (PS1) survey has obtained imaging in five bands (grizyP1) over 10 Medium Deep Survey (MDS) fields covering a total of 70 square degrees. This paper describes the search for apparently hostless supernovae (SNe) within the first year of PS1 MDS data with an aim of discovering superluminous supernovae (SLSNe). A total of 249 hostless transients were discovered down to a limiting magnitude of MAB ∼ 23.5, of which 76 were classified as Type Ia supernovae (SNe Ia). There were 57 SNe with complete light curves that are likely core-collapse SNe (CCSNe) or type Ic SLSNe and 12 of these have had spectra taken. Of these 12 hostless, non-Type Ia SNe, 7 were SLSNe of type Ic at redshifts between 0.5 and 1.4. This illustrates that the discovery rate of type Ic SLSNe can be maximized by concentrating on hostless transients and removing normal SNe Ia. We present data for two possible SLSNe; PS1-10pm (z = 1.206) and PS1-10ahf (z = 1.1), and estimate the rate of type Ic SLSNe to be between 3+3−2×10−53+3−2×10−53+3−2×10−5 and 8+2−1×10−58+2−1×10−58+2−1×10−5 that of the CCSN rate within 0.3 ≤ z ≤ 1.4 by applying a Monte Carlo technique. The rate of slowly evolving, type Ic SLSNe (such as SN2007bi) is estimated as a factor of 10 lower than this range.

Ultraluminous Supernovae as a New Probe of the Interstellar Medium in Distant Galaxies

We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultraluminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z = 1.566 with a peak brightness of M UV ≈ -22.3 mag. PS1-11bam is one of the highest redshift spectroscopically confirmed SNe known to date. The spectrum exhibits broad absorption features typical of previous ULSNe (e.g., C II, Si III), and strong and narrow Mg II and Fe II absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [O II]λ3727 emission line at the same redshift. The equivalent widths of the Fe II λ2600 and Mg II λ2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of τ* ≈ 15-45 Myr and a stellar mass of M * ≈ (1.1-2.6) × 109 M ⊙ (for Z = 0.05-1 Z ⊙). The star formation rate inferred from the UV continuum and [O II]λ3727 emission line is ≈10 M ⊙ yr-1, higher than in previous ULSN hosts. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the ISM in distant galaxies. The depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z ~ 1-2); the future combination of LSST and 30 m class telescopes promises to extend this technique to z ~ 4.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.