Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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,027 records · Page 57

SupportU: Smart UAS Program for the Population by Offering Resources and Tools to the Unhoused

Global warming, challenging economic conditions, the opioid epidemic, and other widespread problems, have impacted people globally, particularly over the past several years. Preventable diseases like the common cold and the effects of heat stroke have become increasingly prevalent due to these issues. It is estimated that 150 million people of the world’s population are unhoused globally, with many dwelling in unsafe and unsanitary conditions while lacking access to basic hygienic items and other essentials. Unsanitary conditions coupled with this lack of access exacerbates and prolongs health problems and harms quality of life. Traditional methods of aid, such as homeless shelters and meal programs, face numerous challenges such as having limited reach and resources. To address these problems, the Smart UAS Program for the Population by Offering Resources and Tools to the Unhoused (SUPPORT U) utilizes Uncrewed Aircraft Systems (UAS) to deliver resources to the unhoused and those in need of basic aid, prior to and during extreme temperature conditions, and after natural disasters in rural, suburban, and urban areas. The UA is envisioned to be an autonomous aircraft capable of efficiently distributing essential supplies, including blankets, water, food, and medicine. The UAS fleet relies on advanced navigation and communication technologies to accurately identify unhoused people and efficiently and safely distribute materials to them. This will be done through a machine learning algorithm. By focusing on identified “homeless clusters”, places where unhoused individuals are concentrated, the UAS network increases access to critical resources, thereby helping to mitigate some of the external risks to the health of unhoused individuals. The UA can also be used during crises, such as by transporting supplies to medical tents that are stationed in difficult-to-reach areas suffering from natural disasters.

Samuel Beard↗

Spin-free Dirac-Coulomb calculations augmented with a perturbative treatment of spin-orbit effects at the Hartree-Fock level

A perturbative approach to compute second-order spin-orbit (SO) corrections to a spin-free Dirac-Coulomb Hartree-Fock (SFDC-HF) calculation is suggested. The proposed scheme treats the difference between the DC and SFDC Hamiltonian as perturbation and exploits analytic second-derivative techniques. In addition, a cost-effective scheme for incorporating relativistic effects in high-accuracy calculations is suggested consisting of a SFDC coupled-cluster treatment augmented by perturbative SO corrections obtained at the HF level. Benchmark calculations for the hydrogen halides HX, X = F-At as well as the coinage-metal fluorides CuF, AgF, and AuF demonstrate the accuracy of the proposed perturbative treatment of SO effects on energies and electrical properties in comparison with the more rigorous full DC treatment. Furthermore, we present, as an application of our scheme, results for the electrical properties of AuF and XeAuF.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Study of Particle Rotation Effect in Gas-Solid Flows using Direct Numerical Simulation with a Lattice Boltzmann Method

A new and efficient direct numerical method with second-order convergence accuracy was developed for fully resolved simulations of incompressible viscous flows laden with rigid particles. The method combines the state-of-the-art immersed boundary method (IBM), the multi-direct forcing method, and the lattice Boltzmann method (LBM). First, the multi-direct forcing method is adopted in the improved IBM to better approximate the no-slip/no-penetration (ns/np) condition on the surface of particles. Second, a slight retraction of the Lagrangian grid from the surface towards the interior of particles with a fraction of the Eulerian grid spacing helps increase the convergence accuracy of the method. An over-relaxation technique in the procedure of multi-direct forcing method and the classical fourth order Runge-Kutta scheme in the coupled fluid-particle interaction were applied. The use of the classical fourth order Runge-Kutta scheme helps the overall IB-LBM achieve the second order accuracy and provides more accurate predictions of the translational and rotational motion of particles. The preexistent code with the first-order convergence rate is updated so that the updated new code can resolve the translational and rotational motion of particles with the second-order convergence rate. The updated code has been validated with several benchmark applications. The efficiency of IBM and thus the efficiency of IB-LBM were improved by reducing the number of the Lagragian markers on particles by using a new formula for the number of Lagrangian markers on particle surfaces. The immersed boundary-lattice Boltzmann method (IBLBM) has been shown to predict correctly the angular velocity of a particle. Prior to examining drag force exerted on a cluster of particles, the updated IB-LBM code along with the new formula for the number of Lagrangian markers has been further validated by solving several theoretical problems. Moreover, the unsteadiness of the drag force is examined when a fluid is accelerated from rest by a constant average pressure gradient toward a steady Stokes flow. The simulation results agree well with the theories for the short- and long-time behavior of the drag force. Flows through non-rotational and rotational spheres in simple cubic arrays and random arrays are simulated over the entire range of packing fractions, and both low and moderate particle Reynolds numbers to compare the simulated results with the literature results and develop a new drag force formula, a new lift force formula, and a new torque formula. Random arrays of solid particles in fluids are generated with Monte Carlo procedure and Zinchenko's method to avoid crystallization of solid particles over high solid volume fractions. A new drag force formula was developed with extensive simulated results to be closely applicable to real processes over the entire range of packing fractions and both low and moderate particle Reynolds numbers. The simulation results indicate that the drag force is barely affected by rotational Reynolds numbers. Drag force is basically unchanged as the angle of the rotating axis varies.

36 MATERIALS SCIENCE↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Iodosylbenzene Coordination Chemistry Relevant to Metal–Organic Framework Catalysis

Hypervalent iodine compounds formally feature expanded valence shells at iodine. These reagents are broadly used in synthetic chemistry due to the ability to participate in well-defined oxidation-reduction processes and because the ligand-exchange chemistry intrinsic to the hypervalent center allows hypervalent iodine compounds to be applied to a broad array of oxidative substrate functionalization reactions. We have recently developed methods to generate these compounds from O 2 that are predicated on diverting reactive intermediates of aldehyde autoxidation towards the oxidation of aryl iodides. Coupling the aerobic oxidation of aryliodides with catalysts that effect C–H bond oxidation would provide a strategy to achieve aerobic C–H oxidation chemistry. In this Forum article we discuss the aspects of hypervalent iodine chemistry and bonding that render this class of reagents attractive lynchpins for aerobic oxidation chemistry. We then discuss the oxidation processes relevant to the aerobic preparation of 2-(tertbutylsulfonyl) iodosylbenzene, which is a popular hypervalent iodine reagent for use with porous metal-organic framework-based catalysts because it displays significantly enhanced solubility as compared with unsubstituted iodosylbenzene. We demonstrate that popular synthetic methods to this reagent often provide material that displays unpredictable disproportionation behavior due to the presence of trace impurities. We provide a revised synthetic route that avoids impurities common in the reported methods and provides access to material that displays predictable stability. Lastly, we describe the coordination chemistry of hypervalent iodine compounds with metal clusters relevant to MOF chemistry and discuss the potential implications of this coordination chemistry to catalysis in MOF scaffolds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-E-T-565: RAdiation Resistance of Cancer CElls Using GEANT4 DNA: RACE

Purpose: The objective of the RACE project is to develop a comparison between Monte Carlo simulation using the Geant4-DNA toolkit and measurements of radiation damage on 3D melanoma and chondrosarcoma culture cells coupled with gadolinium nanoparticles. We currently expose the status of the developments regarding simulations. Methods: Monte Carlo studies are driven using the Geant4 toolkit and the Geant4-DNA extension. In order to model the geometry of a cell population, the opensource CPOP++ program is being developed for the geometrical representation of 3D cell populations including a specific cell mesh coupled with a multi-agent system. Each cell includes cytoplasm and nucleus. The correct modeling of the cell population has been validated with confocal microscopy images of spheroids. The Geant4 Livermore physics models are used to simulate the interactions of a 250 keV X-ray beam and the production of secondaries from gadolinium nanoparticles supposed to be fixed on the cell membranes. Geant4-DNA processes are used to simulate the interactions of charged particles with the cells. An atomistic description of the DNA molecule, from PDB (Protein Data Bank) files, is provided by the so-called PDB4DNA Geant4 user application we developed to score energy depositions in DNA base pairs and sugar-phosphate groups. Results: At the microscopic level, our simulations enable assessing microscopic energy distribution in each cell compartment of a realistic 3D cell population. Dose enhancement factors due to the presence of gadolinium nanoparticles can be estimated. At the nanometer scale, direct damages on nuclear DNA are also estimated. Conclusion: We successfully simulated the impact of direct radiations on a realistic 3D cell population model compatible with microdosimetry calculations using the Geant4-DNA toolkit. Upcoming validation and the future integration of the radiochemistry module of Geant4-DNA will propose to correlate clusters of ionizations with in vitro experiments. All those developments will be released publicly. This work was supported by grants from Plan Cancer 2009-2013 French national initiative managed by INSERM (Institut National de la Sante et de la Recherche Medicale)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Consistent inclusion of triple substitutions within a coupled cluster based static quantum embedding theory

We have previously proposed the MPCC static embedding framework for quantum chemistry that self-consistently couples a high-level coupled cluster (CC) treatment of the fragment (active region) with a lower level, Møller–Plesset perturbation treatment of the environment. Our initial implementation was limited to single and double (SD) substitutions, with CCSD for the fragment and first-order perturbative SD amplitudes for the environment. Here, we extend the MPCC embedding treatment to triple substitutions, which is essential for achieving chemical accuracy in energy differences. To this end, we employ a CCSDT solver for the fragment subsystem. For the environment subsystem, we construct a perturbative estimate of the triples amplitudes, explicitly accounting for feedback from all fragment amplitudes. The resulting approach is denoted MPCCSDT(pt). We further introduce a more complete formulation in which feedback from the environment amplitudes to the fragment amplitudes is also included. This scheme involves an iterative treatment of the environment triples amplitudes and is denoted MPCCSDT(it). In addition, we assess the accuracy of the previously proposed low-level method by introducing a modified low-level approach that incorporates a lowest-order treatment of selected long-range effects, including spin fluctuations and charge polarization. All resulting approaches may be viewed as post-CCSD(T) methods. We therefore consider test cases for which CCSD(T) exhibits substantial deviations from CCSDT. These include (i) single- and triple-bond stretching in F 2 and N 2 , (ii) bond dissociation energies of selected molecules from the W4-11 dataset, and (iii) total atomization energies of transition metal hydrides. Our results demonstrate that inclusion of triples amplitudes at the fragment level alone is insufficient; a perturbative treatment of the environment triples amplitudes is required. For many energy-difference applications, feedback from the environment triples amplitudes to the fragment amplitudes is not essential, but it does play a role in the very challenging CoH and FeH molecules. A very interesting finding from our study is that in some challenging cases, we need an improved (second-order) perturbative method for the SD amplitudes, going beyond the first-order one used in our earlier work. In conclusion, considering both cost and accuracy, the MP2CCSDT(pt) model is the most promising for future applications among the candidates considered here.

Shee, Avijit [University of California, Berkeley, ↗

Heat transport in confined strongly coupled two-dimensional dust clusters

Dusty plasmas are a model system for studying strong correlation. The dust grains’ size of a few micro-meters and their characteristic oscillation frequency of a few hertz allow for an investigation of many-particle effects on an “atomic” level. In this article, we model the heat transport through an axially confined 2D dust cluster from the center to the outside. The system behaves particularly interesting since heat is not only conducted within the dust component but also transferred to the neutral gas. Fitting the analytical solution to the radial temperature profiles obtained in molecular dynamics simulations allows to determine the heat conductivity k. The heat conductivity is found to be constant over a wide range of coupling strengths even including the phase transition from solid to liquid here, as it was also found in extended systems by Nosenko et al.[Phys. Rev. Lett. 100, 025003 (2008)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the implementation of the spherical collapse model for dark energy models

In this work we review the theory of the spherical collapse model and critically analyse the aspects of the numerical implementation of its fundamental equations. By extending a recent work by [1], we show how different aspects, such as the initial integration time, the definition of constant infinity and the criterion for the extrapolation method (how close the inverse of the overdensity has to be to zero at the collapse time) can lead to an erroneous estimation (a few per mill error which translates to a few percent in the mass function) of the key quantity in the spherical collapse model: the linear critical overdensity δ{sub c}, which plays a crucial role for the mass function of halos. We provide a better recipe to adopt in designing a code suitable to a generic smooth dark energy model and we compare our numerical results with analytic predictions for the EdS and the ΛCDM models. We further discuss the evolution of δ{sub c} for selected classes of dark energy models as a general test of the robustness of our implementation. We finally outline which modifications need to be taken into account to extend the code to more general classes of models, such as clustering dark energy models and non-minimally coupled models.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of cloud properties and atmospheric stability with MODIS

The work accomplished with the MAS (MODIS Airborne Simulator) during the TOGA-COARE (Tropical Ocean Global Atmosphere - Coupled Ocean Atmosphere Response Experiment) from January through June 1993 included cloud investigations with longwave CO2 channels (out to 13.9 microns) that reconfigured spectral channels. The comparison data was used with HIRS (high resolution Infrared Radiation Sounder) data from FIRE (First ISSCP Regional Experiment), CEPEX, and SCAR (Smoke Clouds, Aerosol and Radiation) data for the development of a two layer cloud CO2 slicing algorithm and the development of multi-spectral cloud detection algorithms (specifically tri-spectral) for cloud identification. Regarding the latter, a clustering analysis was added to improve cloud identification within the scatter diagram of tri-spectral (8, 11, and 12 microns) brightness temperatures. Via the regressions of the 8 - 11 and the 11 - 12 micron brightness temperature differences versus precipitable water, the problem of clear sky identification in the scatter diagram was further addressed. Comparison of the tri-spectral technique and visible channel reflectance ratioing indicate the tri-spectral infrared method distinguishes cloud phase more consistently in mixed cloud scenes rather than single phase cloud scenes where visible ratioing performs best.

Menzel, Paul↗

Magnetic cluster expansion model for random and ordered magnetic face-centered cubic Fe-Ni-Cr alloys

A Magnetic Cluster Expansion model for ternary face-centered cubic Fe-Ni-Cr alloys has been developed, using DFT data spanning binary and ternary alloy configurations. Using this Magnetic Cluster Expansion model Hamiltonian, we perform Monte Carlo simulations and explore magnetic structures of alloys over the entire range of compositions, considering both random and ordered alloy structures. In random alloys, the removal of magnetic collinearity constraint reduces the total magnetic moment but does not affect the predicted range of compositions where the alloys adopt low-temperature ferromagnetic configurations. During alloying of ordered fcc Fe-Ni compounds with Cr, chromium atoms tend to replace nickel rather than iron atoms. Replacement of Ni by Cr in ordered alloys with high iron content increases the Curie temperature of the alloys. This can be explained by strong antiferromagnetic Fe-Cr coupling, similar to that found in bcc Fe-Cr solutions, where the Curie temperature increase, predicted by simulations as a function of Cr concentration, is confirmed by experimental observations. In random alloys, both magnetization and the Curie temperature decrease abruptly with increasing chromium content, in agreement with experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ambient vapor sampling and selective cluster formation for the trace detection of tributyl phosphate via atmospheric flow tube mass spectrometry

In addition to serving as an f-element ligand and rare-earth method complexing agent, tributyl phosphate (TBP) is a compound containing core functional groups that mimic those routinely found in degradation products from industrial processes. Because detection of trace quantities of TBP can provide insight into the routes of contamination and degradation in the environment, there is a need to develop methods capable of detecting trace quantities of TBP. Vapor detection at atmospheric pressure is one approach that is both sensitive and rapid. We present here the use of atmospheric flow tube mass spectrometry (AFT-MS) for the ambient vapor sampling of TBP from headspace of ppb-level solutions in methanol. Gas phase clustering reactions were to enhance detection levels via the addition of small quantities of the dopants diethylamine (DEA), triethylamine (TEA), and pinacolyl methylphosphonate (PMP) in the vapor stream. Detection of the TBP vapor emanating from these solutions demonstrated a linear range for the protonated TBP species of 1-1000 ppb from solution. The clusters of TBP with DEA, TEA, and PMP each yielded linear ranges of 1-250 ppb for TBP in solution. Despite smaller linear ranges, the addition of these dopant species added a layer of analytical selectivity and reduced variability in signals from quality control samples. These data were obtained using an AFT source coupled to a linear ion trap mass spectrometer, which demonstrates the applicability of trapping systems to the AFT ionization technique while monitoring positive ions.

Ambient vapor sampling, Dielectric barrier dischar↗

First-principles studies on molecular beam epitaxy growth of GaAs 1-x Bi x

We investigate the molecular beam epitaxy (MBE) growth of GaAs 1-x Bi x film using density functional theory with spin-orbit coupling to understand the growth of this film, especially the mechanisms of Bi incorporation. We study the stable adsorption structures and kinetics of the incident molecules (As₂ molecule, Ga atom, Bi atom, and Bi₂ molecule) on the (2 x 1)-Ga sub ||Bi surface and a proposed q (1 x 1)-Ga sub ||AsAs surface has a quasi-(1 x 1) As layer above the Ga-terminated GaAs substrate and a randomly oriented As dimer layer on top. We obtain the desorption and diffusion barriers of the adsorbed molecules and also the reaction barriers of three key processes related to Bi evolution, namely, Bi incorporation, As/Bi exchange, and Bi clustering. The results help explain the experimentally observed dependence of Bi incorporation on the As/Ga ratio and growth temperature. Furthermore, we find that As₂ exchange with Bi of the (2 x 1)-Ga sub ||Bi surface is a key step controlling the kinetics of the Bi incorporation. Finally, we explore two possible methods to enhance the Bi incorporation, namely, replacing the MBE growth mode from codeposition of all fluxes with a sequential deposition of fluxes and applying asymmetric in-plane strain to the substrate.

Luo, Guangfu↗

Theoretical spectroscopic investigations of HNS{sup q} and HSN{sup q} (q = 0, +1, −1) in the gas phase

We performed accurate ab initio investigations of the geometric parameters and the vibrational structure of neutral HNS/HSN triatomics and their singly charged anions and cations. We used standard and explicitly correlated coupled cluster approaches in connection with large basis sets. At the highest levels of description, we show that results nicely approach those obtained at the complete basis set limit. Moreover, we generated the three-dimensional potential energy surfaces (3D PESs) for these molecular entities at the coupled cluster level with singles and doubles and a perturbative treatment of triple excitations, along with a basis set of augmented quintuple-zeta quality (aug-cc-pV5Z). A full set of spectroscopic constants are deduced from these potentials by applying perturbation theory. In addition, these 3D PESs are incorporated into variational treatment of the nuclear motions. The pattern of the lowest vibrational levels and corresponding wavefunctions, up to around 4000 cm{sup −1} above the corresponding potential energy minimum, is presented for the first time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect-induced magnetism in cobalt-doped ZnO epilayers

We used a synergic Co-edge X-ray absorption spectroscopy (XAS) and density functional theory calculations approach to perform a study of defects which could account for the room temperature ferromagnetism of ZnCoO, an oxide of great potential interest in semiconductor spintronics. Our results suggest that a key role is played by specific defect complexes in which O vacancies are located close to the Co atoms. Extended defects such as Co clusters have a marginal function, although we observe their formation at the epilayer surface under certain growth conditions. We also show preliminary results of the study of hydrogen-induced defects in ZnCoO epilayers deliberately hydrogen irradiated via a Kaufman source. Hydrogen was in fact predicted to mediate a ferromagnetic spin-spin interaction between neighboring magnetic impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Filaments in and between galaxy clusters at low and mid-frequency with the SKA

Context. Understanding the magnetised Universe is a major challenge in modern astrophysics, and cosmic magnetism has been acknowledged as one of the key scientific drivers of the most ambitious radio instrument ever planned, the Square Kilometre Array. Aims. With this work, we aim to investigate the potential of the Square Kilometre Array and its precursors and pathfinders in the study of magnetic fields in galaxy clusters and filaments through diffuse synchrotron radio emission. Galaxy clusters and filaments of the cosmic web are indeed unique laboratories in which to investigate turbulent fluid motions and large-scale magnetic fields in action, and much of what is known about magnetic fields in galaxy clusters comes from sensitive radio observations. Methods. Based on cosmological magneto-hydrodynamic simulations, we predict radio properties (total intensity and polarisation) of a pair of galaxy clusters connected by a cosmic-web filament. Results. We use our theoretical expectations to explore the potential of polarimetric observations to study large-scale structure magnetic fields in the frequency ranges 50–350 MHz and 950–1760 MHz. We also present predictions for galaxy cluster polarimetric observations with the Square Kilometre Array precursors and pathfinders, such as the LOw frequency ARray 2.0 and the MeerKAT+ telescope. Conclusions. Our findings point out that polarisation observations are particularly powerful for the study of large-scale magnetic fields, since they are not significantly affected by confusion noise. The unprecedented sensitivity and spatial resolution of the intermediatefrequency radio telescopes make them the favourite instruments for the study of these sources through polarimetric data, potentially allowing us to understand if the energy density of relativistic electrons is in equipartition with the magnetic field or rather coupled with the thermal gas density. Our results show that low-frequency instruments also represent a precious tool to study diffuse synchrotron emission in total intensity and polarisation.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulations of the dissociation of small helium clusters with ab initio molecular dynamics in electronically excited states

The dynamics resulting from electronic excitations of helium clusters were explored using ab initio molecular dynamics. The simulations were performed with configuration interaction singles and adiabatic classical dynamics coupled to a state-following algorithm. 100 different configurations of He{sub 7} were excited into the 2s and 2p manifold for a total of 2800 trajectories. While the most common outcome (90%) was complete fragmentation to 6 ground state atoms and 1 excited state atom, 3% of trajectories yielded bound, He {sub 2}{sup *}, and <0.5% yielded an excited helium trimer. The nature of the dynamics, kinetic energy release, and connections to experiments are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Source characteristics and geological implications of the January 2016 induced earthquake swarm near Crooked Lake, Alberta

On 2016 January 12, an intraplate earthquake with an initial reported local magnitude (M L ) of 4.8 shook the town of Fox Creek, Alberta. While there were no reported damages, this earthquake was widely felt by the local residents and suspected to be induced by the nearby hydraulic-fracturing (HF) operations. We determine the earthquake source parameters using moment tensor inversions, and then detect and locate the associated swarm using a waveform cross-correlation based method. The broad-band seismic recordings from regional arrays suggest a moment magnitude (M) 4.1 for this event, which is the largest in Alberta in the past decade. Similar to other recent M ~ 3 earthquakes near Fox Creek, the 2016 January 12 earthquake exhibits a dominant strike-slip (strike = 184°) mechanism with limited non-double-couple components (~22 per cent). This resolved focal mechanism, which is also supported by forward modelling and P-wave first motion analysis, indicates an NE–SW oriented compressional axis consistent with the maximum compressive horizontal stress orientations delineated from borehole breakouts. Further detection analysis on industry-contributed recordings unveils 1108 smaller events within 3 km radius of the epicentre of the main event, showing a close spatial-temporal relation to a nearby HF well. The majority of the detected events are located above the basement, comparable to the injection depth (3.5 km) on the Duvernay shale Formation. The spatial distribution of this earthquake cluster further suggests that (1) the source of the sequence is an N–S-striking fault system and (2) these earthquakes were induced by an HF well close to but different from the well that triggered a previous (January 2015) earthquake swarm. Reactivation of pre-existing, N-S oriented faults analogous to the Pine Creek fault zone, which was reported by earlier studies of active source seismic and aeromagnetic data, are likely responsible for the occurrence of the January 2016 earthquake swarm and other recent events in the Crooked Lake area.

58 GEOSCIENCES↗