Did galaxies reionize the universe?

The epoch of reionization remains one of the last uncharted eras of cosmic history, yet this time is of crucial importance, encompassing the formation of both the first galaxies and the first metals in the universe. In this thesis, I present four related projects that both characterize the abundance and properties of these first galaxies and uses follow-up observations of these galaxies to achieve one of the first observations of the neutral fraction of the intergalactic medium during the heart of the reionization era.

First, we present the results of a spectroscopic survey using the Keck telescopes targeting 6.3 < z < 8.8 star-forming galaxies. We secured observations of 19 candidates, initially selected by applying the Lyman break technique to infrared imaging data from the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST). This survey builds upon earlier work from Stark et al. (2010, 2011), which showed that star-forming galaxies at 3 < z < 6, when the universe was highly ionized, displayed a significant increase in strong Lyman alpha emission with redshift. Our work uses the LRIS and NIRSPEC instruments to search for Lyman alpha emission in candidates at a greater redshift in the observed near-infrared, in order to discern if this evolution continues, or is quenched by an increase in the neutral fraction of the intergalactic medium. Our spectroscopic observations typically reach a 5-sigma limiting sensitivity of < 50 AA. Despite expecting to detect Lyman alpha at 5-sigma in 7-8 galaxies based on our Monte Carlo simulations, we only achieve secure detections in two of 19 sources. Combining these results with a similar sample of 7 galaxies from Fontana et al. (2010), we determine that these few detections would only occur in < 1% of simulations if the intrinsic distribution was the same as that at z ~ 6. We consider other explanations for this decline, but find the most convincing explanation to be an increase in the neutral fraction of the intergalactic medium. Using theoretical models, we infer a neutral fraction of X_HI ~ 0.44 at z = 7.

Second, we characterize the abundance of star-forming galaxies at z > 6.5 again using WFC3 onboard the HST. This project conducted a detailed search for candidates both in the Hubble Ultra Deep Field as well as a number of additional wider Hubble Space Telescope surveys to construct luminosity functions at both z ~ 7 and 8, reaching 0.65 and 0.25 mag fainter than any previous surveys, respectively. With this increased depth, we achieve some of the most robust constraints on the Schechter function faint end slopes at these redshifts, finding very steep values of alpha_{z~7} = -1.87 +/- 0.18 and alpha_{z~8} = -1.94 +/- 0.23. We discuss these results in the context of cosmic reionization, and show that given reasonable assumptions about the ionizing spectra and escape fraction of ionizing photons, only half the photons needed to maintain reionization are provided by currently observable galaxies at z ~ 7-8. We show that an extension of the luminosity function down to M_{UV} = -13.0, coupled with a low level of star-formation out to higher redshift, can fit all available constraints on the ionization history of the universe.

Third, we investigate the strength of nebular emission in 3 < z < 5 star-forming galaxies. We begin by using the Infrared Array Camera (IRAC) onboard the Spitzer Space Telescope to investigate the strength of H alpha emission in a sample of 3.8 < z < 5.0 spectroscopically confirmed galaxies. We then conduct near-infrared observations of star-forming galaxies at 3 < z < 3.8 to investigate the strength of the [OIII] 4959/5007 and H beta emission lines from the ground using MOSFIRE. In both cases, we uncover near-ubiquitous strong nebular emission, and find excellent agreement between the fluxes derived using the separate methods. For a subset of 9 objects in our MOSFIRE sample that have secure Spitzer IRAC detections, we compare the emission line flux derived from the excess in the K_s band photometry to that derived from direct spectroscopy and find 7 to agree within a factor of 1.6, with only one catastrophic outlier. Finally, for a different subset for which we also have DEIMOS rest-UV spectroscopy, we compare the relative velocities of Lyman alpha and the rest-optical nebular lines which should trace the cites of star-formation. We find a median velocity offset of only v_{Ly alpha} = 149 km/s, significantly less than the 400 km/s observed for star-forming galaxies with weaker Lyman alpha emission at z = 2-3 (Steidel et al. 2010), and show that this decrease can be explained by a decrease in the neutral hydrogen column density covering the galaxy. We discuss how this will imply a lower neutral fraction for a given observed extinction of Lyman alpha when its visibility is used to probe the ionization state of the intergalactic medium.

Finally, we utilize the recent CANDELS wide-field, infra-red photometry over the GOODS-N and S fields to re-analyze the use of Lyman alpha emission to evaluate the neutrality of the intergalactic medium. With this new data, we derive accurate ultraviolet spectral slopes for a sample of 468 3 < z < 6 star-forming galaxies, already observed in the rest-UV with the Keck spectroscopic survey (Stark et al. 2010). We use a Bayesian fitting method which accurately accounts for contamination and obscuration by skylines to derive a relationship between the UV-slope of a galaxy and its intrinsic Lyman alpha equivalent width probability distribution. We then apply this data to spectroscopic surveys during the reionization era, including our own, to accurately interpret the drop in observed Lyman alpha emission. From our most recent such MOSFIRE survey, we also present evidence for the most distant galaxy confirmed through emission line spectroscopy at z = 7.62, as well as a first detection of the CIII]1907/1909 doublet at z > 7.

We conclude the thesis by exploring future prospects and summarizing the results of Robertson et al. (2013). This work synthesizes many of the measurements in this thesis, along with external constraints, to create a model of reionization that fits nearly all available constraints.

Quantum interference effects of two coherent population trapping states on the atomic spectral lines of a F-e=0 -> F-g=1 transition

We investigate the fluorescence spectrum in a nearly degenerate atomic system of a F-e = 0 -> F-g = 1 transition by analytically solving Schrodinger equations. An ultranarrow fluorescence spectral line in between the two coherent population trapping windows has been found. Our analytic solutions clearly show the origin of the ultranarrow spectral line. Due to quantum interference effects between two coherent population trapping states, the width and intensity of the central spectral line can be controlled by an external magnetic field. Such an effect may be used to detect a magnetic field.

Manipulating spectral anomalies of focused pulses in a medium with electromagnetically induced transparency

We study the possibility of manipulating the focusing properties of a medium with electromagnetically induced transparency. In the focal region of focused ultraslow light pulses, the spectral anomalous behaviors can be actively modified by varying the control field intensity. Unlike the case in free space, we find in slow light focusing that the spectrum bandwidth of the incident field needed to produce observable spectral changes can be reduced by several orders. Numerical simulations with accessible parameters clearly show that spectral anomalies of focused mu s pulses are observable.

Sensitivity map of laser tweezers Raman spectroscopy for single-cell analysis of colorectal cancer

Raman spectroscopy on single, living epithelial cells captured in a laser trap is shown to have diagnostic power over colorectal cancer. This new single-cell technology comprises three major components: primary culture processing of human tissue samples to produce single-cell suspensions, Raman detection on singly trapped cells, and diagnoses of the cells by artificial neural network classifications. it is compared with DNA flow cytometry for similarities and differences. Its advantages over tissue Raman spectroscopy are also discussed. In the actual construction of a diagnostic model for colorectal cancer, real patient data were taken to generate a training set of 320 Raman spectra and, a test set of 80. By incorporating outlier corrections to a conventional binary neural classifier, our network accomplished significantly better predictions than logistic regressions, with sensitivity improved from 77.5% to 86.3% and specificity improved from 81.3% to 86.3% for the training set and moderate improvements for the test set. Most important, the network approach enables a sensitivity map analysis to quantitate the relevance of each Raman band to the normal-to-cancer transform at the cell level. Our technique has direct clinic applications for diagnosing cancers and basic science potential in the study of cell dynamics of carcinogenesis. (C) 2007 Society of Photo-Optical Instrumentation Engineers.

Cross spectral purity and its influence on ghost imaging experiments

Pseudo-thermal light has been widely used in ghost imaging experiments. In order to understand the differences between the pseudo-thermal source and thermal source, we propose a method to investigate whether a light source has cross spectral purity (CSP), and experimentally measure the cross spectral properties of the pseudo-thermal light source in near-field and far-field zones. Moreover we present a theoretical analysis of the cross spectral influence on ghost imaging. (c) 2006 Elsevier B.V. All rights reserved.

Ultranarrow absorptive spectral line induced by microwave field

This paper investigates the absorptive spectral lines of four-level atomic system driven by a coupling, probe and microwave fields. Due to the perturbation of the microwave field, the original electromagnetically induced transparency is changed to electromagnetically induced absorption and the absorptive spectral line can be very narrow. This ultranarrow spectral line has potential applications to the microwave atomic frequency standard and the measurement of very weak magnetic field.

A variational framework for spectral discretization of the density matrix in Kohn Sham density functional theory

Kohn-Sham density functional theory (KSDFT) is currently the main work-horse of quantum mechanical calculations in physics, chemistry, and materials science. From a mechanical engineering perspective, we are interested in studying the role of defects in the mechanical properties in materials. In real materials, defects are typically found at very small concentrations e.g., vacancies occur at parts per million, dislocation density in metals ranges from $10^{10} m^{-2}$ to $10^{15} m^{-2}$, and grain sizes vary from nanometers to micrometers in polycrystalline materials, etc. In order to model materials at realistic defect concentrations using DFT, we would need to work with system sizes beyond millions of atoms. Due to the cubic-scaling computational cost with respect to the number of atoms in conventional DFT implementations, such system sizes are unreachable. Since the early 1990s, there has been a huge interest in developing DFT implementations that have linear-scaling computational cost. A promising approach to achieving linear-scaling cost is to approximate the density matrix in KSDFT. The focus of this thesis is to provide a firm mathematical framework to study the convergence of these approximations. We reformulate the Kohn-Sham density functional theory as a nested variational problem in the density matrix, the electrostatic potential, and a field dual to the electron density. The corresponding functional is linear in the density matrix and thus amenable to spectral representation. Based on this reformulation, we introduce a new approximation scheme, called spectral binning, which does not require smoothing of the occupancy function and thus applies at arbitrarily low temperatures. We proof convergence of the approximate solutions with respect to spectral binning and with respect to an additional spatial discretization of the domain. For a standard one-dimensional benchmark problem, we present numerical experiments for which spectral binning exhibits excellent convergence characteristics and outperforms other linear-scaling methods.

Gauge Invariant Spectral Cauchy Characteristic Extraction of Gravitational Waves in Computational General Relativity

We present a complete system for Spectral Cauchy characteristic extraction (Spectral CCE). Implemented in C++ within the Spectral Einstein Code (SpEC), the method employs numerous innovative algorithms to efficiently calculate the Bondi strain, news, and flux.

Spectral CCE was envisioned to ensure physically accurate gravitational wave-forms computed for the Laser Interferometer Gravitational wave Observatory (LIGO) and similar experiments, while working toward a template bank with more than a thousand waveforms to span the binary black hole (BBH) problem’s seven-dimensional parameter space.

The Bondi strain, news, and flux are physical quantities central to efforts to understand and detect astrophysical gravitational wave sources within the Simulations of eXtreme Spacetime (SXS) collaboration, with the ultimate aim of providing the first strong field probe of the Einstein field equation.

In a series of included papers, we demonstrate stability, convergence, and gauge invariance. We also demonstrate agreement between Spectral CCE and the legacy Pitt null code, while achieving a factor of 200 improvement in computational efficiency.

Spectral CCE represents a significant computational advance. It is the foundation upon which further capability will be built, specifically enabling the complete calculation of junk-free, gauge-free, and physically valid waveform data on the fly within SpEC.

Lock Acquisition and Sensitivity Analysis of Advanced LIGO Interferometers

Laser interferometer gravitational wave observatory (LIGO) consists of two complex large-scale laser interferometers designed for direct detection of gravitational waves from distant astrophysical sources in the frequency range 10Hz - 5kHz. Direct detection of space-time ripples will support Einstein's general theory of relativity and provide invaluable information and new insight into physics of the Universe.

Initial phase of LIGO started in 2002, and since then data was collected during six science runs. Instrument sensitivity was improving from run to run due to the effort of commissioning team. Initial LIGO has reached designed sensitivity during the last science run, which ended in October 2010.

In parallel with commissioning and data analysis with the initial detector, LIGO group worked on research and development of the next generation detectors. Major instrument upgrade from initial to advanced LIGO started in 2010 and lasted till 2014.

This thesis describes results of commissioning work done at LIGO Livingston site from 2013 until 2015 in parallel with and after the installation of the instrument. This thesis also discusses new techniques and tools developed at the 40m prototype including adaptive filtering, estimation of quantization noise in digital filters and design of isolation kits for ground seismometers.

The first part of this thesis is devoted to the description of methods for bringing interferometer to the linear regime when collection of data becomes possible. States of longitudinal and angular controls of interferometer degrees of freedom during lock acquisition process and in low noise configuration are discussed in details.

Once interferometer is locked and transitioned to low noise regime, instrument produces astrophysics data that should be calibrated to units of meters or strain. The second part of this thesis describes online calibration technique set up in both observatories to monitor the quality of the collected data in real time. Sensitivity analysis was done to understand and eliminate noise sources of the instrument.

Coupling of noise sources to gravitational wave channel can be reduced if robust feedforward and optimal feedback control loops are implemented. The last part of this thesis describes static and adaptive feedforward noise cancellation techniques applied to Advanced LIGO interferometers and tested at the 40m prototype. Applications of optimal time domain feedback control techniques and estimators to aLIGO control loops are also discussed.

Commissioning work is still ongoing at the sites. First science run of advanced LIGO is planned for September 2015 and will last for 3-4 months. This run will be followed by a set of small instrument upgrades that will be installed on a time scale of few months. Second science run will start in spring 2016 and last for about 6 months. Since current sensitivity of advanced LIGO is already more than factor of 3 higher compared to initial detectors and keeps improving on a monthly basis, upcoming science runs have a good chance for the first direct detection of gravitational waves.

基于3×3耦合器的马赫-曾德尔干涉仪的光纤光栅波长解调技术

介绍了一种基于3×3和2×2光纤耦合器构成的非平衡马赫-曾德尔干涉仪的波长解调方案。理论分析和数据对比表明，相对于由两个2×2光纤耦合器构成的马赫-曾德尔干涉仪，本干涉仪具有宽谱的灵敏度、能跟踪波长的变化方向和相位展开的优点，实验方案用于测量固定在悬臂梁上的传感光纤布拉格光栅（FBG）的峰值波长变化，获得了±1pm的静态波长解调精度，在10Hz处的动态分辨率为27nε/√Hz。相位展开算法使得应变测量范围达到了2014με，对应的相位变化为3．22π。

Sensitivity of sturgeons to environmental hypoxia: a review of physiological and ecological evidence

In this essay, three lines of evidence are developed that sturgeons in the Chesapeake Bay and elsewhere are unusually sensitive to hypoxic conditions: 1. In comparison to other fishes,sturgeons have a limited behavioral and physiological capacity to respond to hypoxia. Basal metabolism, growth, feeding rate, and survival are sensitive to changes in oxygen level, which may indicate a relatively poor ability of sturgeons to oxyregulate. 2. During summertime, temperatures >20°C amplify the effect of hypoxia on sturgeons and other fishes due to a temperature oxygen "squeeze" (Coutant 1987). In bottom waters, this interaction results in substantial reduction of habitat; in dry years, sturgeon nursery habitats in the Chesapeake Bay may be particularly reduced or even eliminated. 3. While evidence for population level effects due to hypoxia is circumstantial, there are corresponding trends between the absence of Atlantic sturgeon reproduction in estuaries like the Chesapeake Bay where summertime hypoxia predominates on a system-wide scale. Also, the recent and dramatic recovery of shortnose sturgeon in the Hudson River (4-bid increase in abundance from 1980 to1995) may have been stimulated by improvement of a large portion of the nursery habitat that was restored from hypoxia to normoxia during the period 1973-1978.

Quantifying Earthquake Collapse Risk of Tall Steel Moment Frame Buildings Using Rupture-to-Rafters Simulations

There is a sparse number of credible source models available from large-magnitude past earthquakes. A stochastic source model generation algorithm thus becomes necessary for robust risk quantification using scenario earthquakes. We present an algorithm that combines the physics of fault ruptures as imaged in laboratory earthquakes with stress estimates on the fault constrained by field observations to generate stochastic source models for large-magnitude (Mw 6.0-8.0) strike-slip earthquakes. The algorithm is validated through a statistical comparison of synthetic ground motion histories from a stochastically generated source model for a magnitude 7.90 earthquake and a kinematic finite-source inversion of an equivalent magnitude past earthquake on a geometrically similar fault. The synthetic dataset comprises of three-component ground motion waveforms, computed at 636 sites in southern California, for ten hypothetical rupture scenarios (five hypocenters, each with two rupture directions) on the southern San Andreas fault. A similar validation exercise is conducted for a magnitude 6.0 earthquake, the lower magnitude limit for the algorithm. Additionally, ground motions from the Mw7.9 earthquake simulations are compared against predictions by the Campbell-Bozorgnia NGA relation as well as the ShakeOut scenario earthquake. The algorithm is then applied to generate fifty source models for a hypothetical magnitude 7.9 earthquake originating at Parkfield, with rupture propagating from north to south (towards Wrightwood), similar to the 1857 Fort Tejon earthquake. Using the spectral element method, three-component ground motion waveforms are computed in the Los Angeles basin for each scenario earthquake and the sensitivity of ground shaking intensity to seismic source parameters (such as the percentage of asperity area relative to the fault area, rupture speed, and risetime) is studied.

Under plausible San Andreas fault earthquakes in the next 30 years, modeled using the stochastic source algorithm, the performance of two 18-story steel moment frame buildings (UBC 1982 and 1997 designs) in southern California is quantified. The approach integrates rupture-to-rafters simulations into the PEER performance based earthquake engineering (PBEE) framework. Using stochastic sources and computational seismic wave propagation, three-component ground motion histories at 636 sites in southern California are generated for sixty scenario earthquakes on the San Andreas fault. The ruptures, with moment magnitudes in the range of 6.0-8.0, are assumed to occur at five locations on the southern section of the fault. Two unilateral rupture propagation directions are considered. The 30-year probabilities of all plausible ruptures in this magnitude range and in that section of the fault, as forecast by the United States Geological Survey, are distributed among these 60 earthquakes based on proximity and moment release. The response of the two 18-story buildings hypothetically located at each of the 636 sites under 3-component shaking from all 60 events is computed using 3-D nonlinear time-history analysis. Using these results, the probability of the structural response exceeding Immediate Occupancy (IO), Life-Safety (LS), and Collapse Prevention (CP) performance levels under San Andreas fault earthquakes over the next thirty years is evaluated.

Furthermore, the conditional and marginal probability distributions of peak ground velocity (PGV) and displacement (PGD) in Los Angeles and surrounding basins due to earthquakes occurring primarily on the mid-section of southern San Andreas fault are determined using Bayesian model class identification. Simulated ground motions at sites within 55-75km from the source from a suite of 60 earthquakes (Mw 6.0 − 8.0) primarily rupturing mid-section of San Andreas fault are considered for PGV and PGD data.

Dynamic full-range Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We propose a technique for dynamic full-range Fourier-domain optical coherence tomography by using sinusoidal phase-modulating interferometry, where both the full-range structural information and depth-resolved dynamic information are obtained. A novel frequency-domain filtering algorithm is proposed to reconstruct a time-dependent complex spectral interferogram from the sinusoidally phase-modulated interferogram detected with a high-rate CCD camera. By taking the amplitude and phase of the inverse Fourier transform of the complex spectral interferogram, a time-dependent full-range cross-sectional image and depth-resolved displacement are obtained. Displacement of a sinusoidally vibrating glass cover slip behind a fixed glass cover slip is measured with subwavelength sensitivity to demonstrate the depth-resolved dynamic imaging capability of our system. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Improvement of sensitivity and refractive-index changes in LiNbO3 : Ce : Cu crystals with UV light

A nonvolatile recording scheme is proposed using LiNbO3:Ce:Cu crystals and modulated UV light to record gratings simultaneously in two centres and using red light to bleach the grating in the shallow centre to realize persistent photorefractive holographic storage. Compared with the normal UV-sensitized nonvolatile holographic system, the amplitude of refractive-index changes is greatly increased and the recording sensitivity is significantly enhanced by recording with UV light in the LiNbO3:Ce:Cu crystals. Based on jointly solving the two-centre material equations and the coupled-wave equations, temporal evolutions of the photorefractive grating and the diffraction effciency are effectively described and numerically analysed. Roles of doping levels and recording-beam intensity are discussed in detail. Theoretical results confirm and predict experimental results.

A new method of measuring the waveguide dispersion factor and the thermo-optic coefficient of long-period fiber gratings

Based on the Mach-Zehnder effect between the core mode and the cladding modes, the interference fringes are formed by a pair of cascaded long-period fiber gratings (CLPFGs). Theoretical analyses show that the spectral spacing and the wavelength of these fringes are functions of the waveguide dispersion factor gamma, which is a characterizing parameter to LPFG and with theoretical and applicational significance. By measuring the characteristics of the transmission spectra of CLPFGs, the absolute value of gamma can be obtained. At the same time, the thermo-optic coefficient of effective refractive index difference between core and cladding modes, p, can also be obtained by measured the temperature sensitivity of these fringes. In the experiments, \gamma\ and mu were measured by this method to be 0.874 and 4.08 x 10(-5) degreesC(-1), respectively, for LPFGs with period of 450 mum and with a HE14 resonant peak at 1554 nm. (C) 2004 Elsevier B.V. All rights reserved.