SEARCH · Search NASA
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.
Momentum accomodation coefficients obtained from low velocity molecule-surface interactions, for space probe drag calculations
This work reports on the detection of γ-ray emission coincident with the supernova remnant (SNR) G5.7–0.1 using data collected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The SNR is believed to be interacting with molecular clouds, based on 1720 MHz hydroxyl (OH) maser emission observations in its direction. This interaction is expected to provide targets for the production of γ-ray emission from π{sup 0}-decay. A γ-ray source was observed in the direction of SNR G5.7–0.1, positioned near the bright γ-ray source SNR W28. We model the emission from radio to γ-ray energies using a one-zone model. Following consideration of both π{sup 0}-decay and leptonically dominated emission scenarios for the MeV–TeV source, we conclude that a considerable component of the γ-ray emission must originate from the π{sup 0}-decay channel. Finally, constraints were placed on the reported ambiguity of the SNR distance through X-ray column density measurements made using XMM-Newton observations. We conclude G5.7–0.1 is a significant γ-ray source positioned at a distance of ∼3 kpc with luminosity in the 0.1–100 GeV range of L{sub γ} ≈ 7.4 × 10{sup 34} erg s{sup −1}.
Conjugated polymer (CP) materials that are considered as an attractive choice for various electronic and optoelectronic applications possess highly heterogeneous complex structures that can span several orders of magnitude in length scale. This is due to the intricate coupling effects of weak secondary interactions and entropic forces that determine molecular organization in bulk polymer materials. Understanding the influence of molecular interactions on the emergence and evolution of nano- and microscale structures is of paramount importance to the guided design and development of condensed CP materials. Such understanding is even more critical for the rational design of hierarchical systems away from equilibrium, where weak molecular interactions can guide the system along competing kinetic pathways toward local energy minima and metastable architectures. In this work, we studied a promising concept for accessing various kinetically stabilized semi-crystalline CP nanostructures that are formed in the process of controlled chain-growth Kumada catalyst-transfer polymerization initiated by a small-molecule catalytic species which can exist in equilibrium between free and aggregated states. Specifically, we chose a small-molecule perylenedicarboximide (PDCI) polymerization catalytic initiator which has a strong tendency toward reversible supramolecular assembly–disassembly that can be controlled by temperature. Addition of a thiophene monomer to a solution of the PDCI initiator starts the chain-growth polymerization process, which produces distinct nanoscale semi-crystalline polythiophene structures formed under predominant kinetic control. Furthermore, we demonstrated that, depending on the reaction temperature (which affects the fine balance between the position of PDCI assembly–disassembly equilibrium, rate of polymerization, and solvent–solute interactions for the growing polythiophene chains), a range of kinetically trapped hierarchically organized semi-crystalline CP systems with substantially varying optoelectronic properties could be obtained. In addition to spectroscopic and electron microscopic studies, in order to better reveal the structural aspects of the generated polymer nanoscale systems, we carried out a series of X-ray diffraction and neutron scattering experiments which indicated complex hierarchical organization in these kinetically stabilized CP assembled materials. We expect that this in situ polymerization-based approach can lead to a general way to expand access to various semi-crystalline CP nanostructured materials with broadly tunable electronic and optical properties.
We report the detection of {gamma}-ray emission coincident with three supernova remnants (SNRs) using data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. W41, G337.7-0.1, and MSH 17-39 are SNRs known to be interacting with molecular clouds, as evidenced by observations of hydroxyl (OH) maser emission at 1720 MHz in their directions and other observational information. SNR shocks are expected to be sites of cosmic-ray acceleration, and clouds of dense material can provide effective targets for production of {gamma}-rays from {pi}{sup 0}-decay. The observations reveal unresolved sources in the direction of G337.7-0.1 and MSH 17-39 and an extended source coincident with W41. We model their broadband emission (radio to {gamma}-ray) using a simple one-zone model, and after considering scenarios in which the MeV-TeV sources originate from either {pi}{sup 0}-decay or leptonic emission, we conclude that the {gamma}-rays are most likely produced through the hadronic channel.
Given the piecewise approach to modeling intermolecular interactions for force fields, they can be difficult to parametrize since they are fit to data like total energies that only indirectly connect to their separable functional forms. Furthermore, by neglecting certain types of molecular interactions such as charge penetration and charge transfer, most classical force fields must rely on, but do not always demonstrate, how cancellation of errors occurs among the remaining molecular interactions accounted for such as exchange repulsion, electrostatics, and polarization. In this work we present the first generation of the (many-body) MB-UCB force field that explicitly accounts for the decomposed molecular interactions commensurate with a variational energy decomposition analysis, including charge transfer, with force field design choices that reduce the computational expense of the MB-UCB potential while remaining accurate. We optimize parameters using only a single water molecule and water cluster data up through pentamers, with no fitting to condensed phase data, and we demonstrate that high accuracy is maintained when the force field is subsequently validated against conformational energies of larger water cluster data sets, radial distribution functions of the liquid phase, and the temperature dependence of thermodynamic and transport water properties. Finally, we conclude that MB-UCB is comparable in performance to MB-Pol but is less expensive and more transferable by eliminating the need to represent short-ranged interactions through large parameter fits to high order polynomials.
Two isostructural 3D lanthanide arenedisulfonate metal-organic frameworks (MOFs) [Ln(Hbidc)(nds){sub 0.5}(H{sub 2}O)]{sub n}(Ln=Eu(1), La(2)) have been successfully synthesized by the hydrothermal reaction of lanthanide oxide with 2,6-naphthalenedisulfonate sodium (Na{sub 2}nds) and an auxiliary ligand, 1H-benzimidazole-5,6-dicarboxylic acid (H{sub 3}bidc). The two complexes are both constructed from 2D [Ln(Hbidc)]{sup +} double layers pillared by nds{sup 2−} ligands to generate 3D (3, 8)-connected open-framework structures with 1D long narrow channels running along the a axis. From topological point of view, the 3D framework is a (3, 8)-connected tfz-d net. The weak interactions including N–H⋯O, O–H⋯O hydrogen bonds and π–π stacking are observed in 1. The 2D IR correlation spectroscopy was applied to study the molecular interactions induced by thermal perturbation. The emission spectra of 1 exhibit the characteristic transition of {sup 5}D{sub 0}→{sup 7}F{sub J}(J=0–4) of Eu(III). - Graphical abstract: Two isostructural 3D (3,8)-connected lanthanide arenedisulfonates were hydrothermally synthesized. The 2D IR correlation spectroscopy was applied to study the molecular interactions induced by thermal perturbation. Display Omitted - Highlights: • The first lanthanide arenedisulfonates incorporating fused-ring aromatic carboxylic acid. • Three-dimensional (3,8)-connected framework with tfz-d network topology. • The emission spectra of 1 exhibit the characteristic transition of {sup 5}D{sub 0}→{sup 7}F{sub J} (J=0–4) of Eu(III). • The 2D IR correlation spectroscopy was applied to study the molecular interactions.
Dexter energy transfer (DET), also known as spin forbidden energy transfer (e.g., 3 D 1 A → 1 D 3 A), is undergoing renewed interest because of its utility in energy science, photocatalysis, and biomedical imaging. The DET reactions enable exciton fission and fusion and have been recently used to facilitate nanoscale imaging of protein assemblies on cell surfaces. While DET and Förster energy transfer rates are often well described using Fermi’s golden rule, the donor–acceptor couplings for DET may be bridge-state-mediated and depend on the electronic structure of and interactions among the donor, bridge, and acceptor. We provide a simple approach to approximate DET couplings using a product of molecular orbital overlap factors. This approach may be used to guide the design of DET systems with targeted energy-transfer kinetics.
MicroRNAs (miRNAs) are small noncoding RNA molecules that regulate gene expression. Among these, members of the let-7 miRNA family control many cell-fate determination genes to influence pluripotency, differentiation, and transformation. Lin28 is a specific, posttranscriptional inhibitor of let-7 biogenesis. We report crystal structures of mouse Lin28 in complex with sequences from let-7d, let-7-f1, and let-7g precursors. The two folded domains of Lin28 recognize two distinct regions of the RNA and are sufficient for inhibition of let-7 in vivo. We also show by NMR spectroscopy that the linker connecting the two folded domains is flexible, accommodating Lin28 binding to diverse let-7 family members. Protein-RNA complex formation imposes specific conformations on both components that could affect downstream recognition by other processing factors. Our data provide a molecular explanation for Lin28 specificity and a model for how it regulates let-7.
Due to various applications enabled by diverse morphologies of self-assembled sequence-defined polymers, controlling the self-assembly of synthetic peptidomimetics into designed morphologies has emerged as a promising route for the development of bioinspired functional materials. Herein, we report morphological control over the assembly of a series of short peptoids, or poly-N-substituted glycines, that contain asymmetric hydrophobic domains. We demonstrate that the inherent flexibility of amphiphilic peptoid bilayers drives assembly polymorphism, resulting in the coexistence of nanosheets, twisted ribbons, and nanofibers three distinct morphologies. By tuning peptoid molecular interactions through variations in sequence design, solution pH, and temperature, we demonstrate precise control over the twisting and folding of peptoid bilayers, enabling the formation of well-defined nanosheets and nanohelices. Molecular dynamics simulations further unravel how the introduction of asymmetric hydrophobic domains enables the flexibility of peptoid bilayers and results in peptoid assembly polymorphism. By tuning peptoid molecular interactions through heating, we further demonstrate the transformation of nanosheets into nanohelices. We envision that our mechanistic investigation of peptoid assembly polymorphism provides a strong foundation for leveraging peptoid sequences and chemistries to achieve controlled molecular interactions, driving the creation of biomimetic materials with tailored morphologies and functionalities.
The supramolecular assembly of molecular uranyl species via halogen-oxo interactions and spectroscopic manifestations thereof are probed in the solid state.
GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.
A new bottom-up procedure for constructing non-conservative (dissipative and stochastic) interactions for dissipative particle dynamics (DPD) models is described and applied to perform hierarchical coarse-graining of a polar molecular liquid (nitromethane). The distant-dependent radial and shear frictions in functional-free form are derived consistently with a chosen form for conservative interactions by matching two-body force-velocity and three-body velocity-velocity correlations along the microscopic trajectories of the centroids of Voronoi cells (clusters), which represent the dissipative particles within the DPD description. The Voronoi tessellation is achieved by application of the K-means clustering algorithm at regular time intervals. Consistently with a notion of many-body DPD, the conservative interactions are determined through the multi-scale coarse-graining (MS-CG) method, which naturally implements a pairwise decomposition of the microscopic free energy. A hierarchy of MS-CG/DPD models starting with one molecule per Voronoi cell and up to 64 molecules per cell is derived. The radial contribution to the friction appears to be dominant for all models. As the Voronoi cell sizes increase, the dissipative forces rapidly become confined to the first coordination shell. For Voronoi cells of two and more molecules the time dependence of the velocity autocorrelation function becomes monotonic and well reproduced by the respective MS-CG/DPD models. A comparative analysis of force and velocity correlations in the atomistic and CG ensembles indicates Markovian behavior with as low as two molecules per dissipative particle. The models with one and two molecules per Voronoi cell yield transport properties (diffusion and shear viscosity) that are in good agreement with the atomistic data. The coarser models produce slower dynamics that can be appreciably attributed to unaccounted dissipation introduced by regular Voronoi re-partitioning as well as by larger numerical errors in mapping out the dissipative forces. The framework presented herein can be used to develop computational models of real liquids which are capable of bridging the atomistic and mesoscopic scales.
The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent Cα model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.
Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.
An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view towards describing the phenomenon. While our main focus lies on the computational aspect of the models, we summarize the common characteristics of biomolecular solvation (e.g., solvent structure, polarization, ion binding, and nonpolar behavior) in order to provide reasonable backgrounds to understand the solvation models.
An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.
The zero-multipole summation method has been developed to efficiently evaluate the electrostatic Coulombic interactions of a point charge system. This summation prevents the electrically non-neutral multipole states that may artificially be generated by a simple cutoff truncation, which often causes large amounts of energetic noise and significant artifacts. The resulting energy function is represented by a constant term plus a simple pairwise summation, using a damped or undamped Coulombic pair potential function along with a polynomial of the distance between each particle pair. Thus, the implementation is straightforward and enables facile applications to high-performance computations. Any higher-order multipole moment can be taken into account in the neutrality principle, and it only affects the degree and coefficients of the polynomial and the constant term. The lowest and second moments correspond respectively to the Wolf zero-charge scheme and the zero-dipole summation scheme, which was previously proposed. Relationships with other non-Ewald methods are discussed, to validate the current method in their contexts. Good numerical efficiencies were easily obtained in the evaluation of Madelung constants of sodium chloride and cesium chloride crystals.
We present the first high-energy X-ray (>10 keV) observations of the non-thermal filament G359.97-0.038 using the Nuclear Spectroscopic Telescope Array (NuSTAR). This filament is one of approximately 20 X-ray filaments of unknown origin located in the central 20 pc region in the Galactic Center near Sgr A*. Its NuSTAR and Chandra broadband spectrum is characterized by a single power law with Γ = 1.3 ± 0.3 that extends from 2 to 50 keV, with an unabsorbed luminosity of 1.3 × 10{sup 33} erg s{sup –1} (d/8 kpc){sup 2} in the 2-8 keV band. Despite possessing a cometary X-ray morphology that is typical of a pulsar wind nebula (PWN) in high-resolution Chandra imaging, our spatially resolved Chandra spectral analysis found no significant spectral softening along the filament as would be expected from particle synchrotron cooling. Coincident radio emission is detected using the Very Large Array at 5.5 and 8.3 GHz. We examine and subsequently discard a PWN or magnetic flux tube as the origin of G359.97-0.038. We use broadband spectral characteristics and a morphological analysis to show that G359.97-0.038 is likely an interaction site between the shell of Sgr A East and an adjacent molecular cloud. This is supported by CS molecular line spectroscopy and the presence of an OH maser.