Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERMEDIATE COUPLING”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Synthesis, Electronic Spectroscopy, and Photochemistry of Methacrolein Oxide: A Four-Carbon Unsaturated Criegee Intermediate from Isoprene Ozonolysis

Ozonolysis of isoprene, one of the most abundant volatile organic compounds in the earth’s atmosphere, generates the four-carbon unsaturated methacrolein oxide (MACR-oxide) Criegee intermediate. The first laboratory synthesis and direct detection of MACR-oxide is achieved through reaction of photolytically generated, resonance-stabilized iodoalkene radicals with oxygen. MACR-oxide is characterized on its first π* ← π electronic transition using a ground-state depletion method. MACR-oxide exhibits a broad UV–visible spectrum peaked at 380 nm with weak oscillatory structure at long wavelengths ascribed to vibrational resonances. Complementary theory predicts two strong π* ← π transitions arising from extended conjugation across MACR-oxide with overlapping contributions from its four conformers. Electronic promotion to the 1 1 ππ* state agrees well with experiment, and results in nonadiabatic coupling and prompt release of O 1 D products observed as anisotropic velocity-map images. Furthermore, this UV–visible detection scheme will enable study of its unimolecular and bimolecular reactions under thermal conditions of relevance to the atmosphere.

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

Soft Landing of Bare Nanoparticles with Controlled Size, Composition, and Morphology

A kinetically-limited physical synthesis method based on magnetron sputtering and gas aggregation has been coupled with size-selection and ion soft landing to prepare bare metal nanoparticles on surfaces with controlled coverage, size, composition, and morphology. Employing atomic force microscopy (AFM) and scanning electron microscopy (SEM), it is demonstrated that the size and coverage of bare nanoparticles soft landed onto flat glassy carbon and silicon as well as stepped graphite surfaces may be controlled through size-selection with a quadrupole mass filter and the length of deposition, respectively. The bare nanoparticles are observed with AFM to bind randomly to the flat glassy carbon surface when soft landed at relatively low coverage (1012 ions). In contrast, on stepped graphite surfaces at intermediate coverage (1013 ions) the soft landed nanoparticles are shown to bind preferentially along step edges forming extended linear chains of particles. At the highest coverage (5 x 1013 ions) examined in this study the nanoparticles are demonstrated with both AFM and SEM to form a continuous film on flat glassy carbon and silicon surfaces. On a graphite surface with defects, however, it is shown with SEM that the presence of localized surface imperfections results in agglomeration of nanoparticles onto these features and the formation of neighboring depletion zones that are devoid of particles. Employing high resolution scanning transmission electron microscopy in the high angular annular dark field imaging mode (STEM-HAADF) and electron energy loss spectroscopy (EELS) it is demonstrated that the magnetron sputtering/gas aggregation synthesis technique produces single metal particles with controlled morphology as well as bimetallic alloy nanoparticles with clearly defined core-shell structure. Therefore, this kinetically-limited physical synthesis technique, when combined with ion soft landing, is a versatile complementary method for preparing a wide range of bare supported nanoparticles with selected properties that are free of the solvent, organic capping agents, and residual reactants present with nanoparticles synthesized in solution.

cluster↗
DOE OSTI2017

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗
DOE OSTI2017

The M33 Synoptic Stellar Survey. II. Mira Variables

We present the discovery of 1847 Mira candidates in the Local Group galaxy M33 using a novel semi-parametric periodogram technique coupled with a random forest classifier. The algorithms were applied to ∼2.4 × 10{sup 5} I -band light curves previously obtained by the M33 Synoptic Stellar Survey. We derive preliminary period–luminosity relations at optical, near-infrared, and mid-infrared wavelengths and compare them to the corresponding relations in the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Q-Band Electron-Nuclear Double Resonance Reveals Out-of-Plane Hydrogen Bonds Stabilize an Anionic Ubisemiquinone in Cytochrome bo 3 from Escherichia coli [Q-Band ENDOR Reveals Out-of-Plane Hydrogen Bonds Stabilize an Anionic Ubisemiquinone in Cytochrome bo 3 from Escherichia coli ]

The respiratory cytochrome bo 3 ubiquinol oxidase from E. coli has a high affinity ubiquinone binding site that stabilizes the one-electron reduced ubisemiquinone (SQ H ), which is a transient intermediate during the electron mediated reduction of O 2 to water. It is known that SQ H is stabilized by two strong hydrogen bonds from R71 and D75 to the ubiquinone carbonyl oxygen O1, and weak hydrogen bonds from H98 and Q101 to O4. In the current work, SQ H was investigated with orientation selective Q-band (~34 GHz) pulsed 1 H ENDOR spectroscopy on fully deuterated cyt bo 3 in an H 2 O solvent so that only exchangeable protons contribute to the observed ENDOR spectra. Simulations of the experimental ENDOR spectra provided the principal values and directions of the hyperfine (hfi) tensors for the two strongly coupled H-bond protons (H1 and H2). For H1, the largest principal component of the proton anisotropic hfi tensor T z' = 11.8 MHz, whereas for H2 T z' = 8.6 MHz. Remarkably, the data show that the direction of the H1 H-bond is nearly perpendicular to the quinone plane (~70° out of plane). The orientation of the second strong hydrogen bond, H2, is out of plane by about 25°. Equilibrium molecular dynamics (MD) simulations on a membrane-embedded model of the cyt bo 3 Q H site show that these H-bond orientations are plausible but do not distinguish which H-bond, from R71 or D75, is nearly perpendicular to the quinone ring. Here, density functional theory (DFT) calculations support that the distances and geometries of the H-bonds to the ubiquinone carbonyl oxygens, along with the measured proton anisotropic hfi couplings, are most compatible with an anionic (deprotonated) ubisemiquinone

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

Structural Characterization of the Reaction and Substrate Specificity Mechanisms of Pathogenic Fungal Acetyl-CoA Synthetases

Acetyl CoA synthetases (ACSs) are Acyl-CoA/NRPS/Luciferase (ANL) superfamily enzymes that couple acetate with CoA to generate acetyl CoA, a key component of central carbon metabolism in eukaryotes and prokaryotes. Normal mammalian cells are not dependent on ACSs, while tumor cells, fungi, and parasites rely on acetate as a precursor for acetyl CoA. Consequently, ACSs have emerged as a potential drug target. As part of a program to develop antifungal ACS inhibitors, we characterized fungal ACSs from five diverse human fungal pathogens using biochemical and structural studies. ACSs catalyze a two-step reaction involving adenylation of acetate followed by thioesterification with CoA. Our structural studies captured each step of these two half-reactions including the acetyl-adenylate intermediate of the first half-reaction in both the adenylation conformation and the thioesterification conformation and thus provide a detailed picture of the reaction mechanism. We also used a systematic series of increasingly larger alkyl adenosine esters as chemical probes to characterize the structural basis of the exquisite ACS specificity for acetate over larger carboxylic acid substrates. Consistent with previous biochemical and genetic data for other enzymes, structures of fungal ACSs with these probes bound show that a key tryptophan residue limits the size of the alkyl binding site and forces larger alkyl chains to adopt high energy conformers, disfavoring their efficient binding. Together, our analysis provides highly detailed structural models for both the reaction mechanism and substrate specificity that should be useful in designing selective inhibitors of eukaryotic ACSs as potential anticancer, antifungal, and antiparasitic drugs.

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

Field tunable magnetic transitions of CsCo 2 (MoO 4 ) 2 (OH): a triangular chain structure with a frustrated geometry

The sawtooth chain compound CsCo 2 (MoO 4 ) 2 (OH) is a complex magnetic system and here, we present a comprehensive series of magnetic and neutron scattering measurements to determine its magnetic phase diagram. The magnetic properties of CsCo 2 (MoO 4 ) 2 (OH) exhibit a strong coupling to the crystal lattice and its magnetic ground state can be easily manipulated by applied magnetic fields. There are two unique Co 2+ ions, base and vertex, with J bb and J bv magnetic exchange. The magnetism is highly anisotropic with the b-axis (chain) along the easy axis and the material orders antiferromagnetically at T N = 5 K. There are two successive metamagnetic transitions, the first at H c 1 = 0.2 kOe into a ferrimagnetic structure, and the other at H c 2 = 20 kOe to a ferromagnetic phase. Heat capacity measurements in various fields support the metamagnetic phase transformations, and the magnetic entropy value is intermediate between S = 3/2 and 1/2 states. The zero field antiferromagnetic phase contains vertex magnetic vectors (Co(1)) aligned parallel to the b-axis, while the base vectors (Co(2)) are canted by 34° and aligned in an opposite direction to the vertex vectors. The spins in parallel adjacent chains align in opposite directions, creating an overall antiferromagnetic structure. Further, at a 3 kOe applied magnetic field, adjacent chains flip by 180° to generate a ferrimagnetic phase. An increase in field gradually induces the Co(1) moment to rotate along the b-axis and align in the same direction with Co(2) generating a ferromagnetic structure. The antiferromagnetic exchange parameters are calculated to be J bb = 0.028 meV and J bv = 0.13 meV, while the interchain exchange parameter is considerably weaker at J ch = (0.0047/N ch ) meV. Our results demonstrate that the CsCo 2 (MoO 4 ) 2 (OH) is a promising candidate to study new physics associated with sawtooth chain magnetism and it encourages further theoretical studies as well as the synthesis of other sawtooth chain structures with different magnetic ions.

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

Voltage-induced magnetic domain evolution in a phase-change material

Applying voltage to metal–insulator transition (MIT) materials allows electrical actuation of the local electronic phase state. In MIT systems that have the electronic order coupled with the magnetic order, voltage switching of the electronic phase state can also enable the electrical manipulation of magnetic properties. In this work, we utilized x-ray magnetic circular dichroism photoemission electron microscopy (XMCD-PEEM) to investigate the control of magnetic domain configurations in ferromagnetic MIT electrical switches. For applied voltages above a threshold value, the XMCD-PEEM images show that the magnetic domains separate into two distinct regions: one with a high contrast (white/black), indicating well-defined micrometer-scale magnetic domains with a component of their magnetization aligned parallel/antiparallel to the x-ray helicity, and the other with different shades of intermediate contrast (gray). Significant changes in magnetic domain configurations upon voltage biasing were only observed in these gray regions. Furthermore, the voltage-induced magnetic domain separation was found to be bias polarity-dependent, with the gray regions expanding from the opposite sample edge when the applied voltage polarity was reversed. This polarity-dependent electrical control of magnetic domain configurations during the MIT switching opens alternative opportunities in memory applications for magnetic MIT switching materials.

42 ENGINEERING↗
DOE OSTI2013

THE IMPACT OF BARS ON DISK BREAKS AS PROBED BY S{sup 4}G IMAGING

We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6 {mu}m images from the Spitzer Survey of Stellar Structure in Galaxies. In particular, we have studied the structural properties of those disks with a broken or down-bending profile. We find that, on average, disks with a genuine single-exponential profile have a scale length and a central surface brightness which are intermediate to those of the inner and outer components of a down-bending disk with the same total stellar mass. In the particular case of barred galaxies, the ratio between the break and the bar radii (R{sub br}/R{sub bar}) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10{sup 10} M{sub Sun }, the distribution is bimodal, peaking at R{sub br}/R{sub bar} {approx} 2 and {approx}3.5. The first peak, which is the most populated one, is linked to the outer Lindblad resonance of the bar, whereas the second one is consistent with a dynamical coupling between the bar and the spiral pattern. For galaxies below 10{sup 10} M{sub Sun }, breaks are found up to {approx}10 R{sub bar}, but we show that they could still be caused by resonances given the rising nature of rotation curves in these low-mass disks. While not ruling out star formation thresholds, our results imply that radial stellar migration induced by non-axisymmetric features can be responsible not only for those breaks at {approx}2 R{sub bar}, but also for many of those found at larger radii.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Mathematical Modeling of Multi-Phenomena Anisotropic Systems: Ejection of Primary Aerosols during the Fast Pyrolysis of Biomass

This study introduces a novel particle model for biomass fast pyrolysis, incorporating an anisotropic cylindrical particle to address mass and energy transport coupled with aerosol ejection, which previous models have overlooked. The main contribution lies in developing a model that considers aerosol generation in anisotropic cylindrical particles for the first time, addressing bubbling dynamics and bursting within the liquid phase. The population balance equation describes bubble dynamics and aerosol formation, capturing phenomena like nucleation, growth, coalescence, and bursting. The model employs the method of moments with bubble volume as an internal variable, substantially reducing computational costs by eliminating dependence on this variable. Results highlight the significant impact of anisotropy and particle size on aerosol ejection: smaller, less elongated particles experience faster heating, quicker conversion, and the increased accumulation of the liquid intermediate phase. Specifically, 1 mm diameter particles yield higher concentrations of metaplast and bio-oil aerosols, exceeding 15%, compared to concentrations below 11% for 3 mm particles. This model provides insights into aerosol structure (volume, surface area), aiding in understanding aerosol reactivity at the reactor scale.

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

Identifying the missing link in catalyst transfer polymerization

Nickel-catalyzed catalyst transfer polycondensation (CTP) of thiophenes is an efficient strategy for the controlled synthesis of polythiophenes. However, a detailed view of its reaction mechanism has remained elusive with unresolved questions regarding the geometry and bonding of critical Ni(0) thiophene intermediates. Herein, we provide experimental and computational evidence of structurally characterized square planar η 2 -Ni(0)–thiophene species and their relevance to the mechanism of CTP. These results confirm the viability of C,C-η 2 bound intermediates in CTP of thiophenes, providing an electronic rationale for the stability of such species, and thus that such processes can proceed as living polymerizations. We further show that C,S-κ 2 species may also be relevant in nickel-catalyzed CTP of thiophenes, providing new avenues for exploitation and optimization.

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

Viscous Analysis of Pulsating Hydrodynamic Instability and Thermal Coupling Liquid-Propellant Combustion

A pulsating form of hydrodynamic instability has recently been shown to arise during liquid-propellant deflagration in those parameter regimes where the pressure-dependent burning rate is characterized by a negative pressure sensitivity. This type of instability can coexist with the classical cellular, or Landau form of hydrodynamic instability, with the occurrence of either dependent on whether the pressure sensitivity is sufficiently large or small in magnitude. For the inviscid problem, it has been shown that, when the burning rate is realistically allowed to depend on temperature as well as pressure, sufficiently large values of the temperature sensitivity relative to the pressure sensitivity causes like pulsating form of hydrodynamic instability to become dominant. In that regime, steady, planar burning becomes intrinsically unstable to pulsating disturbances whose wave numbers are sufficiently small. This analysis is extended to the fully viscous case, where it is shown that although viscosity is stabilizing for intermediate and larger wave number perturbations, the intrinsic pulsating instability for small wave numbers remains. Under these conditions, liquid-propellant combustion is predicted to be characterized by large unsteady cells along the liquid/gas interface.

Margolis, Stephen B.↗
DOE OSTI2025

Off-shell vertices in heavy particle effective theories and B → Dπℓν

We study the modifications to decay amplitudes in heavy to heavy semileptonic decays with multiple hadrons in the final state due to intermediate heavy hadrons being off-shell or having a finite width. Combining Heavy Hadron Chiral Perturbation Theory (HHχPT) with a BCFW on-shell factorization formula, we show that these effects induce O (1/M) corrections to the standard results computed in the narrow-width approximation and therefore are important in extracting form factors from data. A combination of perturbative unitarity, analyticity, and reparameterization invariance fully determine these corrections in terms of known Isgur-Wise functions without the need to introduce new form factors. In doing so, we develop a novel technique to compute the boundary term at complex infinity in the BCFW formula for theories with derivatively coupled scalars. While we have used the $\overline{B}$ → Dπℓν decay as an example, these techniques can generally be applied to effective field theories with (multiple) distinct reference vectors. Article PDF

Chiral Lagrangian↗
DOE OSTI2016

Compositional effects on the ignition of FACE gasolines [Compositional effects on the ignition of FACE gasoline fuels: experiments, surrogate fuel formulation, and chemical kinetic modeling]

As regulatory measures for improved fuel economy and decreased emissions are pushing gasoline engine combustion technologies towards extreme conditions (i.e., boosted and intercooled intake with exhaust gas recirculation), fuel ignition characteristics become increasingly important for enabling stable operation. Here, this study explores the effects of chemical composition on the fundamental ignition behavior of gasoline fuels. Two well-characterized, high-octane, non-oxygenated FACE (Fuels for Advanced Combustion Engines) gasolines, FACE F and FACE G, having similar antiknock indices but different octane sensitivities and chemical compositions are studied. Ignition experiments were conducted in shock tubes and a rapid compression machine (RCM) at nominal pressures of 20 and 40 atm, equivalence ratios of 0.5 and 1.0, and temperatures ranging from 650 to 1270 K. Results at temperatures above 900 K indicate that ignition delay time is similar for these fuels. However, RCM measurements below 900 K demonstrate a stronger negative temperature coefficient behavior for FACE F gasoline having lower octane sensitivity. In addition, RCM pressure profiles under two-stage ignition conditions illustrate that the magnitude of low-temperature heat release (LTHR) increases with decreasing fuel octane sensitivity. However, intermediate-temperature heat release is shown to increase as fuel octane sensitivity increases. Various surrogate fuel mixtures were formulated to conduct chemical kinetic modeling, and complex multicomponent surrogate mixtures were shown to reproduce experimentally observed trends better than simpler two- and three-component mixtures composed of n-heptane, iso-octane, and toluene. Measurements in a Cooperative Fuels Research (CFR) engine demonstrated that the multicomponent surrogates accurately captured the antiknock quality of the FACE gasolines. Simulations were performed using multicomponent surrogates for FACE F and G to reveal the underlying chemical kinetics linking fuel composition with ignition characteristics. Finally, a key discovery of this work is the kinetic coupling between aromatics and naphthenes, which affects the radical pool population and thereby controls ignition.

02 PETROLEUM↗
DOE OSTI2012

Inference of Causal Networks from Time-Varying Transcriptome Data via Sparse Coding

Temporal analysis of genome-wide data can provide insights into the underlying mechanism of the biological processes in two ways. First, grouping the temporal data provides a richer, more robust representation of the underlying processes that are co-regulated. The net result is a significant dimensional reduction of the genome-wide array data into a smaller set of vocabularies for bioinformatics analysis. Second, the computed set of time-course vocabularies can be interrogated for a potential causal network that can shed light on the underlying interactions. The method is coupled with an experiment for investigating responses to high doses of ionizing radiation with and without a small priming dose. From a computational perspective, inference of a causal network can rapidly become computationally intractable with the increasing number of variables. Additionally, from a bioinformatics perspective, larger networks always hinder interpretation. Therefore, our method focuses on inferring the simplest network that is computationally tractable and interpretable. The method first reduces the number of temporal variables through consensus clustering to reveal a small set of temporal templates. It then enforces simplicity in the network configuration through the sparsity constraint, which is further regularized by requiring continuity between consecutive time points. We present intermediate results for each computational step, and apply our method to a time-course transcriptome dataset for a cell line receiving a challenge dose of ionizing radiation with and without a prior priming dose. Our analyses indicate that (i) the priming dose increases the diversity of the computed templates (e.g., diversity of transcriptome signatures); thus, increasing the network complexity; (ii) as a result of the priming dose, there are a number of unique templates with delayed and oscillatory profiles; and (iii) radiation-induced stress responses are enriched through pathway and subnetwork studies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Global Ocean Primary Production Trends in the Modern Ocean Color Satellite Record (1998-2015)

Ocean primary production (PP), representing the uptake of inorganic carbon through photosynthesis, supports marine life and affects carbon exchange with the atmosphere. It is difficult to ascertain its magnitude, variability, and trends due to our inability to measure it directly at large scales. Yet it is paramount for understanding changes in marine health, fisheries, and the global carbon cycle. Using assimilation of ocean color satellite data into an ocean biogeochemical model, we estimate that global net ocean PP has experienced a small but significant decline −0.8 PgC/y (−2.1%) per decade (P < 0.05) in the 18-year satellite record from 1998 to 2015. This decline is associated with shallowing surface mixed layer depth (−2.4% per decade) and decreasing nitrate concentrations (−3.2% per decade). Relative contributions to PP by various types of ocean phytoplankton have changed, with decreases in production by intermediate-sized phytoplankton represented by chlorophytes (−14.3% per decade). This is partially compensated by increases from the unique, more nutrient-efficient, coccolithophores (8.4% per decade). Geographically, the North and Equatorial Indian Oceans are responsible for much of the decline in PP, falling 0.16 and 0.69 PgC/y per decade, respectively. Reduced production by large, fast-growing diatoms along with chlorophytes characterizes the decline here. In contrast, increases in PP are found in the North and North Central Pacific. The increases here are led by chlorophytes in the North Pacific and the small cyanobacteria in the North Central Pacific. These results suggest that the multi-decadal satellite observational record, coupled with an underlying representation of marine biodiversity in a model, can monitor the uptake of carbon by phytoplankton and that changes, although small, are occurring in the global oceans.

Watson W Gregg↗
DOE OSTI2015

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

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

Correlation functions on the Half-BPS Wilson loop: perturbation and hexagonalization

We compute correlation functions of protected primaries on the 1/2-BPS Wilson loop in N=4 super Yang-Mills theory at weak coupling. We first perform direct perturbative computation at one loop in the planar limit and present explicit formulae for general two-, three- and four-point functions. The results for two- and three-point functions as well as four-point functions in special kinematics are in perfect agreement with the localization computation performed in arXiv:1802.05201 . We then analyze the results in view of the integrability-based approach called “hexagonalization”, which was introduced previously to study the correlation functions in the absence of the Wilson loop. In this approach, one decomposes the correlator into fundamental building blocks called “hexagons”, and glues them back together by summing over the intermediate states. Through the comparison, we conjecture that the correlation functions on the Wilson loop can be computed by contracting hexagons with boundary states, where each boundary state represents a segment of the Wilson loop. As a byproduct, we make predictions for the large-charge asymptotics of the structure constants on the Wilson loop. Along the way, we refine the conjecture for the integrability-based approach to the general non-BPS structure constants on the Wilson loop, proposed originally in arXiv:1706.02989 .

1/N Expansion↗