Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactions and forces”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

The Green Sahara: Climate Change, Hydrologic History and Human Occupation

Archaeology can provide insight into interactions of climate change and human activities in sensitive areas such as the Sahara, to the benefit of both disciplines. Such analyses can help set bounds on climate change projections, perhaps identify elements of tipping points, and provide constraints on models. The opportunity exists to more precisely constrain the relationship of natural solar and climate interactions, improving understanding of present and future anthropogenic forcing. We are beginning to explore the relationship of human occupation of the Sahara and long-term solar irradiance variations synergetic with changes in atmospheric-ocean circulation patterns. Archaeological and climate records for the last 12 K years are gaining adequate precision to make such comparisons possible. We employ a range of climate records taken over the globe (e.g. Antarctica, Greenland, Cariaco Basin, West African Ocean cores, records from caves) to identify the timing and spatial patterns affecting Saharan climate to compare with archaeological records. We see correlation in changing ocean temperature patterns approx. contemporaneous with drying of the Sahara approx. 6K years BP. The role of radar images and other remote sensing in this work includes providing a geographically comprehensive geomorphic overview of this key area. Such coverage is becoming available from the Japanese PALSAR radar system, which can guide field work to collect archaeological and climatic data to further constrain the climate change chronology and link to models. Our initial remote sensing efforts concentrate on the Gilf Kebir area of Egypt.

remote sensing↗
NASA NTRS2015

Predicting the Mineral Composition of Dust Aerosols: Representing Key Processes - Part 1

Soil dust aerosols created by wind erosion are typically assigned globally uniform physical and chemical properties within Earth system models, despite known regional variations in the mineral content of the parent soil. Mineral composition of the aerosol particles is important to their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Here, aerosol mineral composition is derived by extending a method that provides the composition of a wet-sieved soil. The extension accounts for measurements showing significant differences between the mineral fractions of the wetsieved soil and the emitted aerosol concentration. For example, some phyllosilicate aerosols are more prevalent at silt sizes, even though they are nearly absent at these diameters in a soil whose aggregates are dispersed by wet sieving. We calculate the emitted mass of each mineral with respect to size by accounting for the disintegration of soil aggregates during wet sieving. These aggregates are emitted during mobilization and fragmentation of the original undispersed soil that is subject to wind erosion. The emitted aggregates are carried far downwind from their parent soil. The soil mineral fractions used to calculate the aggregates also include larger particles that are suspended only in the vicinity of the source. We calculate the emitted size distribution of these particles using a normalized distribution derived from aerosol measurements. In addition, a method is proposed for mixing minerals with small impurities composed of iron oxides. These mixtures are important for transporting iron far from the dust source, because pure iron oxides are more dense and vulnerable to gravitational removal than most minerals comprising dust aerosols. A limited comparison to measurements from North Africa shows that the model extensions result in better agreement, consistent with a more extensive comparison to global observations as well as measurements of elemental composition downwind of the Sahara, as described in companion articles.

aerosols↗
NASA NTRS2015

Systematic Relationships Between Lidar Observables and Sizes And Mineral Composition Of Dust Aerosols

The physical and chemical properties of soil dust aerosol particles fundamentally affect their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates on the surface of dust particles, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Lidar measurements, such as extinction-to-backscatter, color and depolarization ratios, are frequently used to distinguish between aerosol types with different physical and chemical properties. The chemical composition of aerosol particles determines their complex refractive index, hence affecting their backscattering properties. Here we present a study on how dust aerosol backscattering and depolarization properties at wavelengths of 355, 532 and 1064 nm are related to size and complex refractive index, which varies with the mineral composition of the dust. Dust aerosols are represented by collections of spheroids with a range of prolate and oblate aspect ratios and their optical properties are obtained using T-matrix calculations. We find simple, systematic relationships between lidar observables and the dust size and complex refractive index that may aid the use of space-based or airborne lidars for direct retrieval of dust properties or for the evaluation of chemical transport models using forward simulated lidar variables. In addition, we present first results on the spatial variation of forward-simulated lidar variables based on a dust model that accounts for the atmospheric cycle of eight different mineral types plus internal mixtures of seven mineral types with iron oxides, which was recently implemented in the NASA GISS Earth System ModelE2.

Dust↗
DOE OSTI2020

Monte Carlo Global QCD Analyses of the Pion Parton Distribution Functions

As the lightest hadron, the pion presents itself as a dichotomy. While being the pseudo Goldstone boson associated with chiral symmetry breaking, it is simultaneously regarded as the lightest pseudoscalar meson typically composed of constituent quark and antiquark. As the chiral symmetry shares the same global symmetry as QCD, studying the pion and its internal structure is key in understanding one of the four fundamental forces of nature, namely the quantum chromodynamic strong interaction. We study the pion's parton distribution functions (PDFs), which are universal quantities that describe the structure of the pion in terms of its constituent quarks, antiquarks, and gluons. Our goal is to discover what the available data reveal for these universal quantities, namely the PDFs. Through the use of factorization theorems and perturbative QCD, we use Monte Carlo (MC) methods to extract the PDFs from the available Drell-Yan (DY) and leading neutron (LN) data. While the DY process involves two hadrons colliding, the LN electroproduction involves an electron beam incident on a target nucleon. We can constrain well the pion PDFs at large momentum fraction using DY data, and at low momentum fraction using LN data. We also use threshold resummation in the DY process to gather predictable higher order terms associated with the soft gluon radiation because their contributions to the cross section are nontrivial. In this dissertation, we parametrize the pion PDFs and fit those parameters to the available data. We make use of Bayesian inference and use MC techniques to sample the parameter space. Three sets of results are presented. We first extract the pion PDFs by fitting the DY and LN data. Then, we include transverse momentum dependent DY data in conjunction with the DY and LN data to extract the pion PDFs. Finally, we apply various methods of threshold resummation to the DY cross section and extract pion PDFs. The pion PDFs presented here are the first global QCD analyses performed as well as the first MC extracted pion PDFs, which are at the forefront of the pion PDF community.

Barry, Patrick↗
DOE OSTI2026

Surface Reconstruction in Hydrated Amphiphilic Block Copolymer Thin Films Probed by Fluid Cell Atomic Force Microscopy

In many thin film materials, nuanced interplays of interfacial energies control the surface morphology and rearrangement. This work evaluates polymer−solvent interactions and solvent-driven surface reconstructions via ex situ and in situ fluid cell Atomic Force Microscopy (fc-AFM) analysis of amphiphilic block copolymer (BCP) thin films upon exposure to deionized (DI) water. We examine the differences in surface morphology, whole-film swelling, and force response in thin films of polystyrene-block-poly(ethylene oxide) (PS-b-PEO) and polystyrene- block-poly[(allyl glycidyl ether)-co-(ethylene oxide)] (PS-b- P[AGE-co-EO]) processed into standing-up cylinder morphologies perpendicular to a silicon substrate (⊥C). Using Amplitude Modulation AFM (AM-AFM) and Amplitude-Phase Distance (APD) force spectroscopy, this work probes the mechanoresponsive nature of the dynamic surface layers of these films, unveiling surface layer stratification and surface chain rearrangement via minimal tip−sample stimulation. To help rationalize the observed reconfigurations, the energetic driving forces were estimated using the harmonic mean approximations of interfacial energies. Given the nonionizable nature of the minority P(AGE-co-EO) block and the energetic driving forces for chain mobility, this work shows how the elimination of unfavorable PS−water interfaces drives chain rearrangement and coverage of the PS surface by chains of the hydrophilic block. This work highlights considerations for increasing the heterogeneity and complexity of BCP thin films via random blocks and how those changes to local interfacial energies may drive larger scale film morphology reconstructions, with broader implications for tuning interface hydrophilicity.

Copolymers↗
DOE OSTI2022

Tidal frequencies and quasiperiodic subsurface water level variations dominate redox dynamics in a salt marsh system

Salt marshes are hotspots of nutrient processing en route to sensitive coastal environments. Whilst our understanding of these systems has improved over the years, we still have limited knowledge of the spatiotemporal variability of critical biogeochemical drivers within salt marshes. Sea-level rise will continue to force change on salt marsh functioning, highlighting the urgency of filling this knowledge gap. Our study was conducted in a central California estuary experiencing extensive marsh drowning and relative sea-level rise, making it a model system for such an investigation. Here we instrumented three marsh positions subjected to different degrees of tidal inundation (6.7%, 8.9%, and 11.2% of the time for the upper, middle, and lower marsh positions, respectively), providing locations with varied biogeochemical characteristics and hydrological interactions at the site. We continuously monitored redox potential (Eh) at depths of 0.1, 0.3, and 0.5 m, subsurface water levels (WL), and temperature at 0.7 m depth at each marsh position. To understand how drivers of subsurface biogeochemical processes fluctuate across tidal cycles, we used wavelet analyses to explain the interactions between Eh and WL. We found that tidal forcing significantly affects key drivers of biogeochemical processes by imparting controls on Eh variability, likely driving subsurface hydro-biogeochemistry of the salt marsh. Wavelet coherence showed that the Eh-WL relationship is nonlinear, and their lead–lag relationship is variable. We found that precipitation events perturb Eh at depth over timescales of hours, even though WL shows relatively minimal change during events. Furthermore, this work highlights the importance of high frequency in situ measurements, such as Eh, to help explain factors that govern subsurface biogeochemistry and hydrological processes in salt marshes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
NASA NTRS2013

Damping Effects of Drogue Parachutes on Orion Crew Module Dynamics

Currently, simulation predictions of the Orion Crew Module (CM) dynamics with drogue parachutes deployed are under-predicting the amount of damping as seen in free-flight tests. The Apollo Legacy Chute Damping model has been resurrected and applied to the Orion system. The legacy model has been applied to predict CM damping under drogue parachutes for both Vertical Spin Tunnel free flights and the Pad Abort-1 flight test. Comparisons between the legacy Apollo prediction method and test data are favorable. A key hypothesis in the Apollo legacy drogue damping analysis is that the drogue parachutes' net load vector aligns with the CM drogue attachment point velocity vector. This assumption seems reasonable and produces good results, but has never been quantitatively proven. The wake of the CM influences the drogue parachutes, which makes performance predictions of the parachutes difficult. Many of these effects are not currently modeled in the simulations. A forced oscillation test of the CM with parachutes was conducted in the NASA LaRC 20-Ft Vertical Spin Tunnel (VST) to gather additional data to validate and refine the Apollo legacy drogue model. A second loads balance was added to the original Orion VST model to measure the drogue parachute loads independently of the CM. The objective of the test was to identify the contribution of the drogues to CM damping and provide additional information to quantify wake effects and the interactions between the CM and parachutes. The drogue parachute force vector was shown to be highly dependent on the CM wake characteristics. Based on these wind tunnel test data, the Apollo Legacy Chute Damping model was determined to be a sufficient approximation of the parachute dynamics in relationship to the CM dynamics for preliminary entry vehicle system design. More wake effects should be included to better model the system. These results are being used to improve simulation model fidelity of CM flight with drogues deployed, which has been identified by the project as key to a successful Orion Critical Design Review.

Aubuchon, Vanessa V.↗
NASA NTRS1988

The GISS global climate-middle atmosphere model. II - Model variability due to interactions between planetary waves, the mean circulation and gravity wave drag

The variability which arises in the GISS Global Climate-Middle Atmosphere Model on two time scales is reviewed: interannual standard deviations, derived from the five-year control run, and intraseasonal variability as exemplified by statospheric warnings. The model's extratropical variability for both mean fields and eddy statistics appears reasonable when compared with observations, while the tropical wind variability near the stratopause may be excessive possibly, due to inertial oscillations. Both wave 1 and wave 2 warmings develop, with connections to tropospheric forcing. Variability on both time scales results from a complex set of interactions among planetary waves, the mean circulation, and gravity wave drag. Specific examples of these interactions are presented, which imply that variability in gravity wave forcing and drag may be an important component of the variability of the middle atmosphere.

Rind, D.↗
NASA NTRS2016

Damping Effects of Drogue Parachutes on Orion Crew Module Dynamics

Because simulations of the Orion Crew Module (CM) dynamics with drogue parachutes deployed were under-predicting the amount of damping seen in free-flight tests, an attach-point damping model was applied to the Orion system. A key hypothesis in this model is that the drogue parachutes' net load vector aligns with the CM drogue attachment point velocity vector. This assumption seems reasonable and has historically produced good results, but has never been experimentally verified. The wake of the CM influences the drogue parachutes, which makes performance predictions of the parachutes difficult. Many of these effects are not currently modeled in the simulations. A forced oscillation test of the CM with parachutes was conducted in the NASA LaRC 20-Ft Vertical Spin Tunnel (VST) to gather additional data to validate and refine the attach-point damping model. A second loads balance was added to the original Orion VST model to measure the drogue parachute loads independently of the CM. The objective of the test was to identify the contribution of the drogues to CM damping and provide additional information to quantify wake effects and the interactions between the CM and parachutes. The drogue parachute force vector was shown to be highly dependent on the CM wake characteristics. Based on these wind tunnel test data, the attach-point damping model was determined to be a sufficient approximation of the parachute dynamics in relationship to the CM dynamics for preliminary entry vehicle system design. More wake effects should be included to better model the system.

Aubuchon, Vanessa V.↗
DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

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

The migration of fluid droplets and their interactions in a thermal gradient

When materials are processed in free fall, buoyant forces will be substantially reduced. Thus, the buoyant migration of droplets and bubbles which normally occurs on earth is expected to be overshadowed by migration due to other mechanisms in space processing. In particular, capillary forces on droplets due to the variation of interfacial tension around their periphery will play a significant role in governing their motion in space. While such interfacial tension gradients can be caused by thermal, compositional, and/or electrical gradients in the continuous phase, thermal gradients are convenient to use in controlled experimentation. On earth, due to interference from buoyant effects, it is difficult to study thermocapillary migration in sufficient detail. Also, the effects of a thermal gradient on the interactions among droplets are hard to study on Earth. Thus, an orbital facility for conducting experiments on the migration and interactions of fluid droplets in a continuous phase due to the action of a thermal gradient appears attractive.

Subramanian, R. S.↗
DOE OSTI2023

Engineering Supramolecular Hybrid Architectures with Directional Organofluorine Bonds

Understanding how chemical modifications alter the atomic‐scale organization of materials is of fundamental importance in materials engineering and the target of considerable efforts in computational prediction. Incorporating covalent and noncovalent interactions in designing crystals while “piggybacking” on the driving force of molecular self‐assembly has augmented efforts to understand the emergence of complex structures using directed synthesis. In this work, microcrystalline powders of the silver 2‐, 3‐, and 4‐fluorobenzenethiolates are prepared and their structures are resolved by small‐molecule serial femtosecond X‐ray crystallography. These three compounds enable the emergence and role of supramolecular synthons in the crystal structures of 3D metal‐organic chalcogenolates to be examined. The unique divergence in their optoelectronic, morphological, and structural behaviors is assessed. The extent of CHF interactions and their influence on the structure and the observed trends in the thermal stability of the crystals are quantified through theoretical calculations and thermogravimetric analysis.

Kotei, Patience A.↗
NASA NTRS1992

Force production mechanisms of a tangential jet on bodies at high alpha

This work explores the mechanisms by which tangential slot blowing creates forces on a body at a high angle of attack. The study is conducted numerically by solving the three-dimensional, compressible-flow Navier-Stokes equations. A tangent-ogive cylinder configuration is used with the blowing slot located both on the nose and on the cylindrical part of the body. The angle of attack used is 30 deg., the Mach no. is 0.2 and the Reynolds no., based on diameter, is 52,000. Several conclusions were made concerning the physical mechanisms by which the jet interacts with the ambient flowfield to produce a side force: (1) A centrifugal force component is created at the wall due to the momentum of the jet being forced to follow the curvature of the surface. (2) A large amount of vorticity is added to the flowfield by the jet. In the region of the slot, the vorticity has the effect of inducing circulation around the body. Downstream of the slot, the vorticity alters the strength of the nose vortices. (3) The position of the nose vortices can be altered to the jet changing the location of separation. And (4), the jet has the ability to excite unstable behavior producing a global change in the character of the flow.

Font, G. I.↗
NASA NTRS1992

Propulsion simulator for magnetically suspended wind tunnel models

Simulation of propulsion induced aerodynamic forces and moments, which arise as a result of interactions between propulsive jets and the free stream, is one of the most desired capabilities in magnetic suspension wind tunnels. The feasibility of generating exhaust jets of appropriate characteristics onboard magnetically suspended models was examined. Four concepts of remotely operated propulsion simulators was considered. Three conceptual designs involving conventional technologies such as compressed gas cylinders, liquid monopropellants, and solid propellants were developed. The fourth concept, a laser assisted thruster, which can potentially simulate both inlet and exhaust flows, was found to require very high power levels. The measurement of aerodynamic forces and/or moments were demonstrated, including the effects of exhaust jets, in MSBS wind tunnels. Two propulsion simulator models were developed, a small scale and a large scale unit, both using compressed, liquified CO2 as propellant. Results of the wind tunnel tests with the small scale simulator are discussed.

Joshi, P. B.↗
DOE OSTI2025

Evidence for the existence and ecological relevance of fast-cycling mineral-associated organic matter

Longstanding theories and models classify mineral-associated organic matter as the large ( ~ 60%) but slow-cycling and persistent portion of soil organic matter. Strong physico-chemical interactions and diffusion limitations restrict the turnover of mineral-associated organic matter, allowing carbon and nitrogen bound therein to persist in soil for as long as centuries to millennia. However, mineral-associated organic matter is a chemically and functionally diverse pool with a substantial portion cycling at relatively fast (i.e., minutes to years) timescales. Despite a growing body of evidence for the heterogenous and multi-pool nature of mineral-associated organic matter, we lack consensus on how to conceptualize and directly quantify fast-cycling mineral-associated organic matter and its ecological significance. We demonstrate that the dynamic qualities of fast-cycling mineral-associated organic matter vary based on 1) the chemistry of the mineral particles and organic matter, 2) the complex set of interactions between organic matter and the mineral matrix, and 3) the presence and strength of destabilizing forces that lead to decomposition or loss of mineral-associated organic matter (i.e., plant-microbe interactions, agricultural intensification, and climate change). Finally, we discuss potential implications and research opportunities for how we measure, manage, and model the dynamic subfraction of this otherwise persistent pool of soil organic matter.

Biological and medical sciences↗
NASA NTRS2018

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗
DOE OSTI2021

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗