Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First order phase transitions”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Commensuration Effects in Layered Nanoparticle Solids.

We have developed a hierarchical nanoparticle transport simulator (HINTS) and adapted it to study commensuration effects in two classes of nanoparticle (NP) solids: (1) a bilayer NP solid (BNS) with an energy offset and (2) a BNS as part of a field effect transistor (FET). HINTS integrates the ab initio characterization of single NPs with the phonon-assisted tunneling transition model of the NP-NP transitions into a kinetic Monte Carlo based simulation of the charge transport in NP solids. First, we studied a BNS with an interlayer energy offset Delta, possibly caused by a fixed electric field. Our results include the following. First, in the independent energy-offset model, we observed the emergence of commensuration effects when scanning the electron filling factor (FF) across integer values. These commensuration effects were profound as they reduced the mobility by several orders of magnitude. We analyzed these commensuration effects in a five-dimensional parameter space, as a function of the on-site charging energy E-C, energy offset Delta, disorder D, the electron FF, and temperature k(B)T. We demonstrated the complexity of our model by showing that at integer FFs commensuration effects are present in some regions of the parameter space, while they vanish in other regions, thus defining distinct dynamical phases of the model. We determined the phase boundaries between these dynamical phases. Second, using these results as a foundation, we shifted our focus to the experimentally much-studied NP-FETs. NP-FETs are also characterized by an interlayer energy offset Delta, which, in contrast to our first model, is set by the gate voltage V-G and thereby related to the electron FF. We repeated many of our simulations and again demonstrated the emergence of commensuration effects and distinct dynamical phases in these NP-FETs. Notably, the commensuration effects in the NP-FETs showed many similarities to those in the independent energy-offset BNS.

Qu, Luman↗
DOE OSTI2022

Solvent-Driven Transformation of Zn/Cd 2+ -Deoxycholate Assemblies

Deoxycholic acid (DOC) is a unique, biologically derived surfactant with facial amphiphilicity that has been exploited, albeit minimally, in supramolecular assembly of materials. Here, we present the synthesis and structural characterization of three hybrid metal (Zn 2+ and Cd 2+ )-DOC compounds. Analysis by single-crystal Xray diffraction reveals the many interactions that are possible between these facial surfactants and the influence of solvent molecules that drive the assembly of materials. These structures are the first metal- DOC complexes besides those obtained from alkali and alkaline earth metals. We isolated polymeric chains of both Cd and Zn (Zn poly -DOC and Cd poly -DOC) from water. Major interactions between DOC molecules in these phases are hydrophobic in nature. Cd poly -DOC exhibits unique P1 symmetry, with complete interdigitation of the amphiphiles between neighboring polymeric chains. Zn 4 -DOC, obtained from methanol dissolution of Zn poly -DOC, features the OZn 4 tetrahedron, widely known in basic zinc acetate and MOF-5 (metal organic framework). We document a solvent-driven, roomtemperature transition between Zn poly -DOC and Zn 4 -DOC (in both directions) by scanning and transmission electron microscopies in addition to small-angle X-ray scattering, powder X-ray diffraction, and infrared spectroscopy. These studies show the methanoldriven transition of Zn poly -DOC to Zn 4 -DOC occurs via an intermediate with no long-range order of the Zn 4 clusters, indicating the strongest interactions driving assembly are intramolecular. On the contrary, water-driven solid-to-solid transformation from Zn 4 -DOC to Zn poly -DOC exhibits crystal-to-crystal transformation. Zn poly -DOC is robust, easy to synthesize, and comprised of biologically benign components, so we demonstrate dye absorption as a proxy for water treatment applications. It favors absorption of positively charged dyes. These studies advance molecular level knowledge of the supramolecular assembly of facial surfactants that can be exploited in the design of organic-inorganic hybrid materials. This research also highlights the potential of solvent for tuning supramolecular assembly processes, leading to new hybrid materials featuring facial surfactants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Revealing the competition between charge density wave and superconductivity in CsV 3 Sb 5 through uniaxial strain

In this paper we report the effect of uniaxial strain ε applied along the crystalline a axis on the newly discovered kagome superconductor CsV 3 Sb 5 . At ambient conditions, CsV 3 Sb 5 shows a charge-density wave (CDW) transition at T CDW = 94.5K and superconducts below T c = 3.34K. In our paper, when the uniaxial strain ε is varied from –0.90% to 0.90%,T c monotonically increases by ~33% from 3.0 to 4.0 K, giving rise to the empirical relation T c (ε) = 3.4 + 0.56ε + 0.12ε 2 . On the other hand, for ε changing from –0.76% to 1.26%,TCDW decreases monotonically by ~10% from 97.5 to 87.5 K with T CDW (ε) = 94.5 – 4.72ε – 0.60ε 2 . The opposite response of T c and T CDW to the uniaxial strain suggests strong competition between these two orders. Comparison with hydrostatic pressure measurements indicate that it is the change in the c axis that is responsible for these behaviors of the CDW and superconducting transitions, and that the explicit breaking of the sixfold rotational symmetry by strain has a negligible effect. Combined with our first-principles calculations and phenomenological analysis, we conclude that the enhancement in Tc with decreasing c is caused primarily by the suppression of T CDW , rather than strain-induced modifications in the bare superconducting parameters. We propose that the sensitivity of T CDW with respect to the changes in the c axis arises from the impact of the latter on the trilinear coupling between the $M$$^{+}_{1}$ and the $L$$^{–}_{2}$ phonon modes associated with the CDW. Altogether, our paper reveals that the c-axis lattice parameter, which can be controlled by both pressure and uniaxial strain, is a powerful tuning knob for the phase diagram of CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Flare Statistics for Young Stars from a Convolutional Neural Network Analysis of TESS Data

All-sky photometric time-series missions have allowed for the monitoring of thousands of young (t age < 800 Myr) stars in order to understand the evolution of stellar activity. In this work, we developed a convolutional neural network (CNN), stella, specifically trained to find flares in Transiting Exoplanet Survey Satellite (TESS) short-cadence data. We applied the network to 3200 young stars in order to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars t age > 50 Myr across all temperatures T eff ≥ 4000 K, while stars from 2300 ≤ T eff < 4000 K show no evolution across 800 Myr. Stars of T eff ≤ 4000 K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4%–7% more atmosphere over the first 1 Gyr. stella is an open-source Python toolkit hosted on GitHub and PyPI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2020

Flare Statistics for Young Stars from a Convolutional Neural Network Analysis of TESS Data

All-sky photometric time-series missions have allowed for the monitoring of thousands of young (t(age) < 800 Myr) stars in order to understand the evolution of stellar activity. Here, we developed a convolutional neural network (CNN), stella, specifically trained to find flares in Transiting Exoplanet Survey Satellite (TESS) short-cadence data. We applied the network to 3200 young stars in order to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars t(age) > 50 Myr across all temperatures T(eff) ≥ 4000 K, while stars from 2300 ≤ T(eff) < 4000 K show no evolution across 800 Myr. Stars of T(eff) ≤ 4000 K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4%–7% more atmosphere over the first 1 Gyr. stella is an open-source Python toolkit hosted on GitHub and PyPI.

Adina D. Feinstein↗
DOE OSTI2021

HJ-aggregates of donor–acceptor–donor oligomers and polymers

A vibronic exciton model is developed to account for the spectral signatures of HJ-aggregates of oligomers and polymers containing donor–acceptor–donor (DAD) repeat units. In (DAD) N π-stacks, J-aggregate-promoting intrachain interactions compete with H-aggregate promoting interchain interactions. The latter includes Coulombic coupling, which arises from “side-by-side” fragment transition dipole moments as well as intermolecular charge transfer (ICT), which is enhanced in geometries with substantial overlap between donors on one chain and acceptors on a neighboring chain. J-behavior is dominant in single (DAD)N chains with enhanced intrachain order as evidenced by an increased red-shift in the low-energy absorption band along with a heightened A 1 /A 2 peak ratio, where A 1 and A 2 are the oscillator strengths of the first two vibronic peaks in the progression sourced by the symmetric quinoidal–aromatic vibration. By contrast, the positive H-promoting interchain Coulomb interactions operative in aggregates cause the vibronic ratio to attenuate, similar to what has been established in H-aggregates of homopolymers such as P3HT. An attenuated A 1 /A 2 ratio can also be caused by H-promoting ICT which occurs when the electron and hole transfer integrals are out-of-phase. In this case, the A 1 peak is red-shifted, in contrast to conventional Kasha H-aggregates. With slight modifications, the ratio formula derived previously for P3HT aggregates is shown to apply to (DAD) N aggregates as well, allowing one to determine the effective free-exciton interchain coupling from the A 1 /A 2 ratio. Applications are made to polymers based on 2T-DPP-2T and 2T-BT-2T repeat units, where the importance of the admixture of the excited acceptor state in the lowest energy band is emphasized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

Noncoherent DTTLs for Symbol Synchronization

Noncoherent data-transition tracking loops (DTTLs) have been proposed for use as symbol synchronizers in digital communication receivers. [Communication- receiver subsystems that can perform their assigned functions in the absence of synchronization with the phases of their carrier signals ( carrier synchronization ) are denoted by the term noncoherent, while receiver subsystems that cannot function without carrier synchronization are said to be coherent. ] The proposal applies, more specifically, to receivers of binary phase-shift-keying (BPSK) signals generated by directly phase-modulating binary non-return-to-zero (NRZ) data streams onto carrier signals having known frequencies but unknown phases. The proposed noncoherent DTTLs would be modified versions of traditional DTTLs, which are coherent. The symbol-synchronization problem is essentially the problem of recovering symbol timing from a received signal. In the traditional, coherent approach to symbol synchronization, it is necessary to establish carrier synchronization in order to recover symbol timing. A traditional DTTL effects an iterative process in which it first generates an estimate of the carrier phase in the absence of symbol-synchronization information, then uses the carrier-phase estimate to obtain an estimate of the symbol-synchronization information, then feeds the symbol-synchronization estimate back to the carrier-phase-estimation subprocess. In a noncoherent symbol-synchronization process, there is no need for carrier synchronization and, hence, no need for iteration between carrier-synchronization and symbol- synchronization subprocesses. The proposed noncoherent symbolsynchronization process is justified theoretically by a mathematical derivation that starts from a maximum a posteriori (MAP) method of estimation of symbol timing utilized in traditional, coherent DTTLs. In that MAP method, one chooses the value of a variable of interest (in this case, the offset in the estimated symbol timing) that causes a likelihood function of symbol estimates over some number of symbol periods to assume a maximum value. In terms that are necessarily oversimplified to fit within the space available for this article, it can be said that the mathematical derivation involves a modified interpretation of the likelihood function that lends itself to noncoherent DTTLs. The proposal encompasses both linear and nonlinear noncoherent DTTLs. The performances of both have been computationally simulated; for comparison, the performances of linear and nonlinear coherent DTTLs have also been computationally simulated. The results of these simulations show that, among other things, the expected mean-square timing errors of coherent and noncoherent DTTLs are relatively insensitive to window width. The results also show that at high signal-to-noise ratios (SNRs), the performances of the noncoherent DTTLs approach those of their coherent counterparts at, while at low SNRs, the noncoherent DTTLs incur penalties of the order of 1.5 to 2 dB.

Simon, Marvin↗
DOE OSTI2024

Pressure dependence of intermediate-range order and elastic properties of glassy Baltic amber

Amber is a unique example of a fragile glass that has been extensively aged below its glass transition temperature, thus reaching a state that is not accessible under normal experimental conditions. In this paper, we studied the medium-range order of Baltic amber by x-ray diffraction (XRD) at high pressures. The pressure dependences of the low-angle XRD intensity between 0 and 5 Å -1 were measured from 0 to 7.3 GPa by the energy-dispersive XRD. The first diffraction peak at 1.1 Å -1 and ambient pressure has a doublet structure consisting of the first sharp diffraction peak (FSDP) at 1.05 Å -1 and the second feature at 1.40 Å -1 . The peak position and the width of the FSDP increase as the pressure increases, while the intensity of the FSDP decreases. Below P 0 = 2.4 GPa, the rapid increase of the FSDP peak position was observed, while above P 0 , the gradual increase was observed. Below P 0 , voids and holes in a relatively low-density state are suppressed, whereas above P 0 , the suppression becomes mild. Such a change suggests the crossover from the low- to high-density state at P 0 . There is a close correlation between the pressure dependence of XRD and previously reported sound velocity results. The correlation between the mean-square fluctuation of the shear modulus on the nanometer scale and fragility in amber and other glass formers is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Theoretical Investigation of Solvent Effects on the Hydrodeoxygenation of Propionic Acid over a Ni(111) Catalyst Model

The effect of two solvents, liquid water and 1,4-dioxane, has been studied from first principles on the hydrodeoxygenation of propionic acid over a Ni (111) catalyst surface model. A mean-field microkinetic model was developed to investigate these effects at a temperature of 473 K. Under all reaction conditions, a decarbonylation mechanism is favored significantly over a decarboxylation pathway. Although no significant solvent effects were observed on the decarbonylation rate, a substantial solvent stabilization of two key surface intermediates in the decarboxylation mechanism, CH 3 CCOO and CH 3 CHCOO, lead to a notable increase of the decarboxylation rate by two orders of magnitude in liquid water and by one order of magnitude in liquid 1,4-dioxane. Furthermore, a significant solvent stabilization of the transition state of C-H bond cleavage of the α-carbon of CH 3 CHCO, relative to the stabilization of the C-C bond cleavage of the α-carbon of CH 3 CHCO, leads to a change in dominant pathway in the liquid phase environments. Finally, a sensitivity analysis shows that the C-OH bond cleavage of propionic acid and C-C bond cleavage of the α-carbon of CH 3 CHCO are the most rate controlling states in the gas phase. In contrast, in solvents the dehydrogenation of CH 3 CHCO becomes the most influential step. This shift in rate controlling state is attributed to the solvent effect on the dehydrogenation of CH 3 CHCO, which is facilitated in aqueous phase. Altogether, it is likely that the investigated (111) facet of Ni is not active for the hydrodeoxygenation of propionic acid in neither the gas nor liquid phase and other Ni facets or phases must be responsible for the experimentally observed kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

The Structure of Ice Nanoclusters: Implications for Interstellar Ice Grains

Crystalline ice nanoclusters and thin-films of pure and impure water ice were deposited and studied in order to evaluate the extent to which surface-related effects control bulk properties. In pure water ice nanoclusters and thin-films of impure water ice, the cubic to hexagonal phase transformation occurs at lower temperatures than in thin-film deposits of pure water ice. In laboratory-grown crystalline ice nanoclusters, approx. 20 nm diameter, a significant proportion of water molecules exists in surface and near-surface environments which have an amorphous or nearly amorphous character. These disordered regions, which are highly reactive, serve to promote transformations or reactions which would otherwise be kinetically hindered. Likewise, dilute impurities such as methanol, sequestered to the surface of thin-films on crystallization, introduce defects into the ice network, thereby allowing sluggish structural transitions to proceed. These structural effects, which are surface phenomena, are believed to play an important role in promoting molecular reactions known to occur within interstellar ice grains in cold molecular clouds, where the first organic compounds are formed.

Delzeit, L.↗
DOE OSTI2021

Antiphase boundary migration as a diffusion mechanism in a P3 sodium layered oxide

The electrode materials of Na-ion batteries must exhibit fast ion transport, a property that is often negatively impacted when mobile cations order. Many layered Na (and K) intercalation compounds adopt the P3 crystal structure, which tends to stabilize cation-vacancy orderings over a wide range of Na concentrations. The ordered phases have been predicted to consist of domains that are periodically separated by antiphase boundaries (APBs) to accommodate variations in Na concentration. The mechanisms and rates with which Na diffuses within such ordered phases of P3 are not understood. We have conducted a first-principles study of Na diffusion in the P3 structure, using Na x CoO 2 as a model system. We find that most simple nearest-neighbor hops in Na-vacancy ordered phases in P3 result in unstable configurations and are therefore forbidden. We identify an alternative mechanism of Na transport that is mediated by APB migration via the formation and expansion of kinks along domain boundaries. The limiting kinetic barriers for this mechanism are quite low (0.03 to 0.30 eV). Vacancy defects are not required, and furthermore have formation energies that are higher than those of the APB kink defects. Furthermore, our results suggest that APB migration is a fundamental diffusion mechanism in the P3 structure.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗
DOE OSTI2022

Structural magnetic glassiness in the spin ice Dy 2 Ti 2 O 7

The origin and nature of glassy dynamics presents one of the central enigmas of condensed-matter physics across a broad range of systems ranging from window glass to spin glasses. The spin-ice compound Dy 2 Ti 2 O 7 , which is perhaps best known as hosting a three-dimensional Coulomb spin liquid with magnetically charged monopole excitations, also falls out of equilibrium at low temperature. How and why it does so remains an open question. Based on an analysis of low-temperature diffuse neutron-scattering experiments employing different cooling protocols alongside recent magnetic noise studies, combined with extensive numerical modeling, we argue that upon cooling, the spins freeze into what may be termed a “structural magnetic glass,” without an a priori need for chemical or structural disorder. Specifically, our model indicates the presence of frustration on two levels, first producing a near-degenerate constrained manifold inside which phase ordering kinetics is in turn frustrated. A remarkable feature is that monopoles act as sole annealers of the spin network and their pathways and history encode the development of glass dynamics, allowing the glass formation to be visualized. Our results suggest that spin ice Dy 2 Ti 2 O 7 provides one prototype of magnetic glass formation specifically and a setting for the study of kinetically constrained systems more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Unprecedented non-hysteretic superelasticity of [001]-oriented NiCoFeGa single crystals

Superelasticity associated with the martensitic transformation has found a broad range of engineering applications. However, the intrinsic hysteresis and temperature sensitivity of the first-order phase transformation significantly hinder the usage of smart metallic components in many critical areas. In this paper, we report a large superelasticity up to 15.2% strain in [001]-oriented NiCoFeGa single crystals, exhibiting non-hysteretic mechanical responses, a small temperature dependence and high-energy-storage capability and cyclic stability over a wide temperature and composition range. In situ synchrotron X-ray diffraction measurements show that the superelasticity is correlated with a stress-induced continuous variation of lattice parameter accompanied by structural fluctuation. Neutron diffraction and electron microscopy observations reveal an unprecedented microstructure consisting of atomic-level entanglement of ordered and disordered crystal structures, which can be manipulated to tune the superelasticity. The discovery of the large elasticity related to the entangled structure paves the way for exploiting elastic strain engineering and development of related functional materials. NiCoFeGa single crystals exhibit large non-hysteretic superelasticity over broad temperature and composition ranges. It is attributed to the continuous phase transition with applied stress, which is related to the fluctuation of entangled ordered and disordered crystal structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetic properties of the Shastry-Sutherland lattice material BaNd 2 ZnO 5

We investigate the physical properties of the Shastry-Sutherland lattice material BaNd 2 ZnO 5 . Neutron diffraction, magnetic susceptibility, and specific heat measurements reveal antiferromagnetic order below 1.65 K. The magnetic order is found to be a 2–Q magnetic structure with the magnetic moments lying in the Shastry-Sutherland lattice planes comprising the tetragonal crystal structure of BaNd 2 ZnO 5 . The ordered moment for this structure is 1.9(1) μB per Nd ion. Inelastic neutron scattering measurements reveal that the crystal field ground state doublet is well separated from the first excited state at 8 meV. The crystal field Hamiltonian is determined through simultaneous refinement of models with both the LS coupling and intermediate coupling approximations to the inelastic neutron scattering and magnetic susceptibility data. The ground state doublet indicates that the magnetic moments lie primarily in the basal plane with magnitude consistent with the size of the determined ordered moment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Architecture and Dynamics of Kepler's Multi-Transiting Planet Systems.III. Comprehensive Investigation Using All Four Years of Kepler Mission Data

We study the orbital architectures of planetary systems orbiting within ∼1 AU of their stars by analyzing the ensemble of Kepler systems having two or more planet candidates. We use data from the entire Kepler mission, and in many cases we apply improved analysis techniques (e.g., replacing histograms by top-hat Kernel Density Estimators that avoid the loss of information resulting from choosing a particular phase for the bin boundaries) to extend and enhance the studies of Lissauer et al. (2011, ApJS 197, 8) and Fabrycky et al. (2014, ApJ 790, 146).These data show ~ 1700 transiting planet candidates in > 600 multiple-planet systems, far more than were available for our previous two studies. The increased numbers and better information about planetary radii and the properties of stellar hosts made possible by Gaia DR2 allow more statistically-robust analyses of the entire ensemble of Kepler multis as well as independent analyses of subsets of the population. We are thus able to contrast the dynamical configurations of small and large planets, short-period and longer-period planets, and planets orbiting various types of host stars. We reinforce our previous findings that most pairs of planets within the same system are neither in nor near low-order mean motion resonances and that there is a substantial excess of planets having period ratios slightly larger than those of first-order mean-motion resonances. However, neglecting three systems whose planets are locked in 3- body resonances and summing over all first-order mean motion resonances, the deficit of planet pairs with period ratios just narrow of resonance is as large as the excess of planets wide of resonance (within statistical uncertainties), suggesting that overall there is no overall excess of planet pairs in the vicinity of resonance. Other aspects of our study, including estimates of the typical relative inclinations of planetary orbits and their variations as functions of orbital period, planet sizes and stellar properties, are in progress, with results expected to be available for presentation by the time of the conference.

Lissauer, Jack↗
NASA NTRS2014

Measurements and Modeling of Total Solar Irradiance in X-class Solar Flares

The Total Irradiance Monitor (TIM) from NASA's SOlar Radiation and Climate Experiment can detect changes in the total solar irradiance (TSI) to a precision of 2 ppm, allowing observations of variations due to the largest X-class solar flares for the first time. Presented here is a robust algorithm for determining the radiative output in the TIM TSI measurements, in both the impulsive and gradual phases, for the four solar flares presented in Woods et al., as well as an additional flare measured on 2006 December 6. The radiative outputs for both phases of these five flares are then compared to the vacuum ultraviolet (VUV) irradiance output from the Flare Irradiance Spectral Model (FISM) in order to derive an empirical relationship between the FISM VUV model and the TIM TSI data output to estimate the TSI radiative output for eight other X-class flares. This model provides the basis for the bolometric energy estimates for the solar flares analyzed in the Emslie et al. study.

Sun: chromosphere↗