Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “small sample size”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

A protected annealing strategy to enhanced light emission and photostability of YAG:Ce nanoparticle-based films

A significant obstacle in the development of YAG:Ce nanoparticles as light converters in white LEDs and as biological labels is associated with the difficulty of finding preparative conditions that allow simultaneous control of structure, particle size and size distribution, while maintaining the optical properties of bulk samples. Preparation conditions frequently involve high-temperature treatments of precursors (up to 1400 C), which result in increased particle size and aggregation, and lead to oxidation of Ce(III) to Ce(IV). We report here a process that we term protected annealing, that allows the thermal treatment of preformed precursor particles at temperatures up to 1000 C while preserving their small size and state of dispersion. In a first step, pristine nanoparticles are prepared by a glycothermal reaction, leading to a mixture of YAG and boehmite crystalline phases. The preformed nanoparticles are then dispersed in a porous silica. Annealing of the composite material at 1000 C is followed by dissolution of the amorphous silica by hydrofluoric acid to recover the annealed particles as a colloidal dispersion. This simple process allows completion of YAG crystallization while preserving their small size. The redox state of Ce ions can be controlled through the annealing atmosphere. The obtained particles of YAG:Ce (60 {+-} 10 nm in size) can be dispersed as nearly transparent aqueous suspensions, with a luminescence quantum yield of 60%. Transparent YAG:Ce nanoparticle-based films of micron thickness can be deposited on glass substrates using aerosol spraying. Films formed from particles prepared by the protected annealing strategy display significantly improved photostability over particles that have not been subject to such annealing.

36 MATERIALS SCIENCE↗

Focusing adaptive-optics for neutron spectroscopy at extreme conditions

Neutron Spectroscopy employing extreme-conditions sample environments is nowadays a crucial tool for the understanding of fundamental scientific questions as well as for the investigation of materials and chemical-physical properties. For all these kinds of studies, an increased neutron flux over a small sample area is needed. The prototype of a focusing neutron guide component, developed and produced completely at the neutron source FRM II in Garching (Germany), has been installed at the time-of-flight (TOF) disc-chopper neutron spectrometer TOFTOF and came into routine-operation. The design is based on the compressed Archimedes' mirror concept for finite-size divergent sources. It represents a unique device combining the supermirror technology with Adaptive Optics, suitable for broad-bandwidth thermal-cold TOF neutron spectroscopy (here optimized for 1.4–10 Å). It is able to squeeze the beam cross section down to a square centimeter, with a more than doubled signal-to-background ratio, increased efficiency at high scattering angles, and improved symmetry of the elastic resolution function. We present a comparison between the simulated and measured beam cross sections, as well as the performance of the instrument within real experiments. This work intends to show the unprecedented opportunities achievable at already existing instruments, along with useful guidelines for the design and construction of next-generation neutron spectrometers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prototypes of self-powered radiation detectors employing intrinsic high-energy current

Purpose: The authors experimentally investigate the effect of direct energy conversion of x-rays via selfpowered Auger- and photocurrent, potentially suitable to practical radiation detection and dosimetry in medical applications. Experimental results are compared to computational predictions. The detector the authors consider is a thin-film multilayer device, composed of alternating disparate electrically conductive and insulating layers. This paper focuses on the experiments while a companion paper introduces the fundamental concepts of high-energy current (HEC) detectors. Methods: The energy of ionizing radiation is directly converted to detector signal via electric current induced by high-energy secondary electrons generated in the detector material by the incident primary radiation. The HEC electrons also ionize the dielectric and the resultant charge carriers are selfcollected due to the contact potential of the disparate electrodes. Thus, an electric current is induced in the conductors in two different ways without the need for externally applied bias voltage or amplification. Thus, generated signal in turn is digitized by a data acquisition system. To determine the fundamental properties of the HEC detector and to demonstrate its feasibility for medical applications, the authors used a planar geometry composed of multilayer microstructures. Various detectors with up to seven conducting layers with different combinations of materials (250 μm Al, 35 μm Cu, 100 μm Pb) and air gaps (100 μm) were exposed to nearly plane-parallel 60–120 kVp x-ray beams. For the experimental design and verification, the authors performed coupled electron–photon radiation transport computations. The detector signal was measured using a commercial data acquisition system with 24 bits dynamic range, 0.4 fC sensitivity, and 0.9 ms sampling time. Results: Measured signals for the prototype detector varied depending on the number of layers, material type, and incident photon energy, and it was in the range of 30–150 nA/cm{sup 2} for unit air kerma (1 Gy), which is viable for practical applications. The experiments had an excellent agreement with the computations. Within the examined range of 60–120 kVp, the energy dependence of the HEC (normalized to the x-ray tube output) was relatively small. Conclusions: Based on the experimental results for 100 ms sampling time, it would be possible to measure the time dependence of x-ray beams for x-ray tube current of 0.1 mA or higher. Significant advantages of the HEC device are that generation of its signal does not require external power supply, it can be made in any size and shape, including flexible curvilinear forms, and it is inexpensive. It remains to be determined, which of the potential applications in medical dosimetry (both in vivo and external), or radiation protection would benefit from such selfpowered detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Clarifying remotely-retrieved precipitation of shallow marine clouds from the NSF/NCAR Gulfstream V

Precipitation is a key process within the shallow cloud lifecycle. The Cloud System Evolution in the Trades (CSET) campaign included the first deployment of a 94 GHz Doppler radar and 532 nm lidar. Despite a larger sampling volume, initial mean radar/lidar retrieved rain rates (Schwartz et al. 2019) based on the upward-pointing remote sensor datasets are systematically less than those measured by in-situ precipitation probes in the cumulus regime. Subsequent retrieval improvements produce rainrates that compare better to in-situ values, but still underestimate. Retrieved shallow cumulus drop sizes can remain too small and too few, with an overestimated shape parameter narrowing the raindrop size distribution too much. Three potential causes for the discrepancy are explored: the gamma functional fit to the dropsize distribution, attenuation by rain and cloud water, and an underaccounting of Mie dampening of the reflectivity. A truncated exponential fit may represent the dropsizes below a showering cumulus cloud more realistically, although further work would be needed to fully evaluate the impact of a different dropsize representation upon the retrieval. The rain attenuation is within the measurement uncertainty of the radar. Mie dampening of the reflectivity is shown to be significant, in contrast to previous stratocumulus campaigns with lighter rain rates, and may be difficult to constrain well with the remote measurements. An alternative approach combines an a priori determination of the dropsize distribution width based on the in-situ data with the mean radar Doppler velocity and reflectivity. This can produce realistic retrievals, although a more comprehensive assessment is needed to better characterize the retrieval errors.

54 ENVIRONMENTAL SCIENCES↗

Surfactant-Free Synthesis of Three-Dimensional Perovskite Titania-Based Micron-Scale Motifs Used as Catalytic Supports for the Methanol Oxidation Reaction

We synthesized and subsequently rationalized the formation of a series of 3D hierarchical metal oxide spherical motifs. Specifically, we varied the chemical composition within a family of ATiO 3 (wherein “A” = Ca, Sr, and Ba) perovskites, using a two-step, surfactant-free synthesis procedure to generate structures with average diameters of ~3 microns. In terms of demonstrating the practicality of these perovskite materials, we have explored their use as supports for the methanol oxidation reaction (MOR) as a function of their size, morphology, and chemical composition. The MOR activity of our target systems was found to increase with decreasing ionic radius of the “A” site cation, in order of Pt/CaTiO 3 (CTO) > Pt/SrTiO 3 (STO) > Pt/BaTiO 3 (BTO). With respect to morphology, we observed an MOR enhancement of our 3D spherical motifs, as compared with either ultra-small or cubic control samples. Moreover, the Pt/CTO sample yielded not only improved mass and specific activity values but also a greater stability and durability, as compared with both commercial TiO 2 nanoparticle standards and precursor TiO 2 templates.

59 BASIC BIOLOGICAL SCIENCES↗

Correlation function studies for snow and ice

The random medium model is used to characterize snow and ice fields in the interpretation of active and passive microwave remote sensing data. A correlation function is used to describe the random permittivity fluctuations with the associated mean and variance and correlation lengths; and several samples are investigated to determine typical correlation functions for snow and ice. It is shown that correlation functions are extracted directly from appropriate ground truth data, and an exponential correlation function is observed for snow and ice with lengths corresponding to the actual size of ice particles or air bubbles. Thus, given that a medium has spatially stationary statistics and a small medium, the random medium model can interpret remote sensing data where theoretical parameters correspond to actual physical parameters of the terrain.

Vallese, F.↗

Th-230 - U-238 series disequilibrium of the Olkaria basalts Gregory Rift Valley, Kenya: Petrogenesis

Strong mixing trends on a (Th-230/Th-232) versus Th diagram show that the basalts are mixed magmas which have undergone interaction with the crust. Instantaneous Th/U ratios are less than time integrated ones but these exceed the Th/U ratios in the MORB and OIB sources. This indicates that the mantle may have undergone some metasomatic fluxing, crustal contamination of the basalts will also enhance these ratios. Early activity on the Akira plain is represented by early basalts and hawaiites. The early basalt samples are known to predate the earliest comendites. The most recent phase of activity is represented by another cinder cone 40-50 m high being feldspar and clinopyroxene phyric. Inclusions which occur in the comendites vary in size and distribution. The largest and most porphyritic are the trachytes (up to 40 cm) with alkali feldspar phases up to 6 mm and small pyroxenes in the ground mass. The second set of inclusions are smaller (up to 10 cm) and are largely aphyric. The distribution of the inclusions are not uniform, the Broad Acres (C5) lavas contain 2-5 percent. The size of the inclusions decrease from south to north, as does the abundance of the trachytic inclusions. The major element variations in the Naivasha basalts, hawaiites and magmatic inclusions are discussed.

Black, S.↗

Occurrence and Stabilization of Colloidal and Nanoparticulate Mercury Sulfides (Final Technical Report)

Mercury contaminates the soils at several DOE sites, including facilities at Oak Ridge, Hanford, and the Savannah River site. This contamination poses a hazard to surrounding ecosystems because the metal continuously leaches to local streams and rivers and is converted to methylmercury, a neurotoxic and highly bioaccumulative form of mercury. Anaerobic microorganisms such as sulfate-reducing and iron-reducing bacteria are primarily responsible for methylmercury production in the environment. The remediation of mercury-contaminated ecosystems is difficult because we lack a poor understanding of the factors that control mercury mobilization and methylation in sediments. The objective of this project was to investigate the forms of mercury that are bioavailable to methylating bacteria. This project specifically focused on mercury-sulfide-organic matter geochemistry and the persistence of nanoparticulate HgS that are intermediates of metal sulfide precipitation reactions occurring in the presence of natural organic matter (NOM). This work involved two specific aims: 1) Investigate environmental processes that facilitate the occurrence and stabilization of nanoparticulate HgS and other metal sulfides; and (2) Investigate the occurrence of nanoparticulate HgS in sediments of the East Fork Poplar Creek, located in near a DOE facility in Oak Ridge, TN. For specific aim #1, laboratory experiments were performed to simulate metal sulfide precipitation processes in environmental systems. A key experimental variable was the type and concentration of NOM that was used in the experiments simulating metal sulfide precipitation. The characteristics of the metal sulfide nanoparticles were analyzed using a suite of analytical techniques that included X-ray based methods at DOE synchrotron facilities (National Synchrotron Light Source and Advanced Photon Source). Field sampling in Oak Ridge was performed in collaboration with colleague at the Oak Ridge National Laboratory. The experimental results demonstrated that a combination of photon scattering and spectroscopic tools can be used to characterize nanoscale structures that form during heterogeneous metal-sulfide-NOM precipitation. ZnS precipitation experiments with a low molecular weight thiol indicated the importance of specific metal-ligand binding for altering the kinetics of ZnS nanoparticle aggregation and monomer growth. However, NOM consists of macromolecular structures that may interact with nano-HgS in non-specific ways (e.g. electrosteric surface interactions. The methylation potential of nano-HgS in sediments may depend on size and crystallinity of nanoparticles and interactions with NOM. For specific aim #2, sediment core samples were collected in three separate locations along the East Fork Poplar Creek in Oak Ridge, TN in June 2011. Analysis of the sediment samples indicated that more than 90% of the Hg in the porewater is associated with colloidal material. These results highlight the importance of small particles for dominating Hg speciation in sediment porewater. Overall, this project highlights the importance of nanoscale phases for mercury speciation. Nanoparticulate HgS will be of particular importance in the anaerobic habitats of methylating microorganisms. In previous models for Hg bioavailability, nanoparticles of Hg are presumed to be unavailable for methylating bacteria. However, recent work shows that nanoparticles contribute to bioavailable Hg that greatly exceeds the contribution from larger HgS microparticles (Zhang et al., 2012). The work performed for this exploratory project lays the groundwork for further study on nanoscale mercury sulfide geochemistry and contributes to an improved understanding of Hg bioavailability and methylation potential.

54 ENVIRONMENTAL SCIENCES↗

(U)SAXS characterization of porous microstructure of chert: insights into organic matter preservation

This study characterizes the microstructure and mineralogy of 132 (ODP sample), 1000 and 1880 million-year-old chert samples. By using ultra-small-angle X-ray scattering (USAXS), wide-angle X-ray scattering and other techniques, the preservation of organic matter (OM) in these samples is studied. The scarce microstructural data reported on chert contrast with many studies addressing porosity evolution in other sedimentary rocks. The aim of this work is to solve the distribution of OM and silica in chert by characterizing samples before and after combustion to pinpoint the OM distribution inside the porous silica matrix. The samples are predominantly composed of alpha quartz and show increasing crystallite sizes up to 33 ± 5 nm (1σ standard deviation or SD). In older samples, low water abundances (~0.03%) suggest progressive dehydration. (U)SAXS data reveal a porous matrix that evolves over geological time, including, from younger to older samples, (1) a decreasing pore volume down to 1%, (2) greater pore sizes hosting OM, (3) decreasing specific surface area values from younger (9.3 ± 0.1 m 2 g -1 ) to older samples (0.63 ± 0.07 m 2 g -1 , 1σ SD) and (4) a lower background intensity correlated to decreasing hydrogen abundances. The pore-volume distributions (PVDs) show that pores ranging from 4 to 100 nm accumulate the greater volume fraction of OM. Raman data show aromatic organic clusters up to 20 nm in older samples. Raman and PVD data suggest that OM is located mostly in mesopores. Observed structural changes, silica–OM interactions and the hydro­phobicity of the OM could explain the OM preservation in chert.

(U)SAXS↗

Structural Reorganization of Noncellulosic Polymers Observed In Situ during Dilute Acid Pretreatment by Small-Angle Neutron Scattering

Production of second-generation bioethanol from lignocellulosic biomass requires pretreatment to open the plant cell wall structure and improve enzyme access. Many different thermochemical pretreatments have been extensively developed and employed, but the exact nature of plant cell wall recalcitrance and the most efficient and economical approach to alter plant cell wall structure via pretreatment still remain elusive. In order to understand the role of noncellulosic switchgrass polymers on the overall efficiency of pretreatment, the structural evolution of the noncellulosic polymers of the plant cell wall were investigated during dilute acid pretreatment (DAP) by employing in-situ small-angle neutron scattering (in-situ SANS). In this study, we observed real-time structural changes not possible to observe by any other technique. To deconvolute the structural evolution of lignin and hemicellulose polymers during DAP, native switchgrass (NATV), and isolated holocellulose (HOLO) and cellulose (CELL) fractions from NATV were studied. Our results show that aggregate particles first appear around 80 °C for NATV and HOLO samples. Further, due to the low temperature and pretreatment severity condition, these particles are likely derived from hemicellulose. The formations of much larger aggregate particles, only observed in the NATV sample, were attributed to lignin. For the HOLO sample, as the temperature and pretreatment severity condition increased, hemicellulose-derived aggregate particle sizes increased, suggesting this process was the nucleation and early stage formation of pseudolignin particles. Consistent with our interpretation of structural evolutions in NATV and HOLO samples, no formation of aggregate particles was observed in CELL samples for the entire duration of the pretreatment. These results suggest that not only lignin but also hemicellulose can form aggregate particles within plant cell walls during pretreatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global geologic map of asteroid (101955) Bennu indicates heterogeneous resurfacing in the past 500,000 years

Global geologic maps are useful tools for efficient interpretation of a planetary body, and they provide global context for the diversity and evolution of the surface. We used data acquired by the OSIRIS-REx spacecraft to create the first global geologic map of the near-Earth asteroid (101955) Bennu. As this is the first geologic map of a small, non-spherical, rubble-pile asteroid, we discuss the distinctive mapping challenges and best practices that may be useful for future exploration of similar asteroids, such as those to be visited with the Hera and Janus missions. By mapping on two centimeter-scale global image mosaics (2D projected space) and a centimeter-scale global shape model (3D space), we generated three input maps respectively describing Bennu’s shape features, geologic features, and surface texture. Based on these input maps, we defined two geologic units: the Smooth Unit and the Rugged Unit. The units are differentiated primarily on the basis of surface texture, concentrations of boulders, and the distributions of lineaments, mass movement features, and craters. They are bounded by several scarps. The Rugged Unit contains abundant boulders and signs of recent mass movement. It also has fewer small(<20 m), putatively fresh craters than the Smooth Unit, suggesting that such craters have been erased in the former. Based on these geologic indicators, we infer that the Rugged Unit has the younger surface of the two. Differential crater size-frequency distributions and the distribution of the freshest craters suggest that both unit surfaces formed ~10–65 million years ago, when Bennu was located in the Main Asteroid Belt, and the Smooth Unit has not been significantly resurfaced in the past 2 million years. Meanwhile, the Rugged Unit has experienced resurfacing within the past ~500,000years during Bennu’s lifetime as a near-Earth asteroid. The geologic units are consistent with global diversity in slope, surface roughness, normal albedo, and thermal emission spectral characteristics. The site on Bennu where the OSIRIS-REx mission collected a regolith sample is located in the Smooth Unit, in a small crater nested within a larger one. So although the Smooth Unit is an older surface than the Rugged Unit, the impact-crater settingindicates that the material sampled was recently exposed. Several similarities are apparent between Bennu and asteroid (162173) Ryugu from a global geologic perspective, including two geologic units distinguishable by variations in the number density of boulders, as well as in other datasets such as brightness.

E R Jawin↗

Further Investigation of Employing Binomial Small-Sample Estimator in M-LROM

Reduced Order Modeling (ROM) is essential to enable the identification of the so-called active parameter subspace for high-fidelity models that are employed in computationally intensive analyses such as uncertainty characterization and data assimilation. Given the complexity of models for real-world application, brute force methods to identify the active subspace may not be computationally efficient to realize the value of reduction, since in realistic problems, the full scale model must be executed few hundred times, which may be computationally taxing. Earlier work has introduced a new approach, which represents a further development of our previous work on ROM, referred to as multi-level ROM (M-LROM). In M-LROM, approximate models are employed to identify the active subspace instead of the full scale model. In this summary, we present our recent results on M-LROM, specifically in regard to the construction of reliable error metrics that can bound the errors resulting from the reduction. The idea is to employ a lower fidelity model, i.e., an approximation of the full scale or high fidelity model, to build the active subspace, and via the proposed error bound, the quality of the active subspace can be determined by calculating an upper bound on the errors resulting from the reduction. For this approach to be practical, the error bound must be realistic, implying that it should be close to the actual errors expected, and not much higher. Therefore, an important goal of this summary is to compare different error bounds in terms of the relation between their estimated bounds to the actual errors. In particular, we compare the proposed error bound which numerically construct an estimator using binomial distribution to the widely used small-sample norm estimators. Unlike most probabilistic error estimation methods, the proposed bound employs a binomial distribution to construct the estimator whereas in other conventional methods, the distribution employed is typically the standard Gaussian, normalized over a unit hypersphere as initially proposed by the seminal paper of J. Dixon on random matrix theory. Previous publications provided evidence that the binomial distribution calculates more realistic bounds as compared to the Gaussian distribution, which tends to be conservative, implying that the bound is much bigger than the actual errors, which diminishes the value of the reduction, as now must increase the size of the active subspace to reach the desired upper bound on the reduction errors. This work has provided a rigorous approach by which error bounds on the reduction errors can be established without having to execute the high fidelity model during the construction of the active subspace. The high fidelity model is only need to test the error bound using small number of executions. This methodology paves the way for a more efficient application of ROM to complex engineering models, where standard ROM techniques may be too taxing to optimally realize their value in practical situations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

New Detections of Radio Minihalos in Cool Cores of Galaxy Clusters

Cool cores of some galaxy clusters exhibit faint radio minihalos. Their origin is unclear, and their study has been limited by their small number. We undertook a systematic search for minihalos in a large sample of X-ray luminous clusters with high-quality radio data. In this article, we report four new minihalos (A 478, ZwCl 3146,RXJ 1532.9+3021, and A 2204) and five candidates found in the reanalyzed archival Very Large Array observations.The radio luminosities of our minihalos and candidates are in the range of 102325 W Hz1 at 1.4 GHz, which is consistent with these types of radio sources. Their sizes (40160 kpc in radius) are somewhat smaller than those of previously known minihalos. We combine our new detections with previously known minihalos, obtaining a total sample of 21 objects, and briefly compare the cluster radio properties to the average X-ray temperature and the total masses estimated from Planck.We find that nearly all clusters hosting minihalos are hot and massive. Beyond that, there is no clear correlation between the minihalo radio power and cluster temperature or mass (in contrast with the giant radio halos found in cluster mergers, whose radio luminosity correlates with the cluster mass). Chandra X-ray images indicate gas sloshing in the cool cores of most of our clusters, with minihalos contained within the sloshing regions in many of them. This supports the hypothesis that radio-emitting electrons are reaccelerated by sloshing. Advection of relativistic electrons by the sloshing gas may also play a role in the formation of the less extended minihalos.

X-ray Polarization↗

Polymorphism influences singlet fission rates in tetracene thin films

Here, we report the effect of crystal structure and crystallite grain size on singlet fission (SF) in polycrystalline tetracene, one of the most widely studied SF and organic semiconductor materials. SF has been comprehensively studied in one polymoprh (Tc I), but not in the other, less stable polymorph (Tc II). Using carefully controlled thermal evaporation deposition conditions and high sensitivity ultrafast transient absorption spectroscopy, we found that for large crystallite size samples, SF in nearly pure Tc II films is significantly faster than SF in Tc I films. We also discovered that crystallite size has a minimal impact on the SF rate in Tc II films, but a significant influence in Tc I films. Large crystallites exhibit SF times of 125 ps and 22 ps in Tc I and Tc II, respectively, whereas small crystallites have SF times of 31 ps and 33 ps. Our results demonstrate first, that attention must be paid to polymorphism in obtaining a self-consistent rate picture for SF in tetracene and second, that control of polymorphism can play a significant role towards achieving a mechanistic understanding of SF in polycrystalline systems. In this latter context we show that conventional theory based on non-covalent tetracene couplings is insufficient, thus highlighting the need for models that capture the delocalized and highly mobile nature of excited states in elucidating the full photophysical picture.

14 SOLAR ENERGY↗

Parallel continuous simulated tempering and its applications in large-scale molecular simulations

In this paper, we introduce a parallel continuous simulated tempering (PCST) method for enhanced sampling in studying large complex systems. It mainly inherits the continuous simulated tempering (CST) method in our previous studies [C. Zhang and J. Ma, J. Chem. Phys. 130, 194112 (2009); C. Zhang and J. Ma, J. Chem. Phys. 132, 244101 (2010)], while adopts the spirit of parallel tempering (PT), or replica exchange method, by employing multiple copies with different temperature distributions. Differing from conventional PT methods, despite the large stride of total temperature range, the PCST method requires very few copies of simulations, typically 2–3 copies, yet it is still capable of maintaining a high rate of exchange between neighboring copies. Furthermore, in PCST method, the size of the system does not dramatically affect the number of copy needed because the exchange rate is independent of total potential energy, thus providing an enormous advantage over conventional PT methods in studying very large systems. The sampling efficiency of PCST was tested in two-dimensional Ising model, Lennard-Jones liquid and all-atom folding simulation of a small globular protein trp-cage in explicit solvent. The results demonstrate that the PCST method significantly improves sampling efficiency compared with other methods and it is particularly effective in simulating systems with long relaxation time or correlation time. We expect the PCST method to be a good alternative to parallel tempering methods in simulating large systems such as phase transition and dynamics of macromolecules in explicit solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Confronting simulations of optically thick gas in massive halos with observations at z = 2-3

Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high-resolution 'zoom-in' simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M {sub vir} ∼ 2 × 10{sup 11} to 4 × 10{sup 12} M {sub ☉} at z ∼ 2. We find that both the mass fraction of cool (T ≤ 3 × 10{sup 4} K) gas and the covering fraction of optically thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f {sub c} ∼ 0.05-0.2, with significant scatter among halos. Our simulations of massive halos (M {sub vir} ≥ 10{sup 12} M {sub ☉}) underpredict the covering fraction of optically thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M {sub vir} ≥ 5 × 10{sup 11} M {sub ☉}) hosting star-forming galaxies, the predicted covering factor agrees with observations; however, current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale autocorrelation function of optically thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z ∼ 2-3, as well as the characteristic size and covering factor of the CGM.

79 ASTRONOMY AND ASTROPHYSICS↗