Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular geometry”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular-Scale Insights into the Heterogeneous Interactions between an m -Terphenyl Isocyanide Ligand and Noble Metal Nanoparticles

The structural and chemical properties of metal nanoparticles are often dictated by their interactions with molecular ligand shells. These interactions are highly material-specific and can vary significantly even among elements within the same group or materials with similar crystal structure. In this study, we surveyed the heterogeneous interactions between an m-terphenyl isocyanide ligand and Au and Ag nanoparticles (NPs) at the single-molecule limit. Specifically, we found that the ligation behavior with this molecule differs significantly between that of Au and AgNPs. Surface-enhanced Raman spectroscopy measurements revealed unique enhancement factors for two molecular vibrational modes between two metal surfaces, indicating different ligand binding geometries. Molecular-level characterization using scanning tunneling microscopy allowed us to directly visualize these variations between Ag and Au surfaces, which we assign as two distinct binding mechanisms. This molecular-scale visualization provides clear insights into the different ligand–metal interactions as well as the chemical behavior and spectroscopic characteristics of isocyanide-functionalized NPs.

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

Computational optimal transport for molecular spectra: The fully continuous case

Computational optimal transport is used to analyze the difference between pairs of continuous molecular spectra. It is demonstrated that transport distances which are derived from this approach may be a more appropriate measure of the difference between two continuous spectra than more familiar measures of distance under many common circumstances. Associated with the transport distances is the transport map which provides a detailed analysis of the difference between two molecular spectra and is a key component of our study of quantitative differences between two continuous spectra. The use of optimal transport for comparing molecular spectra is developed in detail here with a set of model spectra, so that the discussion is self-contained. The difference between the transport distance and more common definitions of distance is elucidated for some well-chosen examples and it is shown where transport distances may be very useful alternatives to standard definitions of distance. The transport distance between a theoretical and experimental electronic absorption spectrum for SO 2 is studied and it is shown how the theoretical spectrum can be modified to fit the experimental spectrum better adjusting the theoretical band origin and the resolution of the theoretical spectrum. In conclusion, this analysis includes the calculation of transport maps between the theoretical and experimental spectra suggesting future applications of the methodology.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Fused-ring isomerism modulates molecular packing and device performance in non-halogenated organic solar cells

Subtle changes in molecular backbone geometry impact intermolecular interactions and performance of organic solar cells. Here, three isomeric small-molecule acceptors (NaO1, NaO2, and NaO3) are investigated to reveal how different fused-ring configurations control molecular packing, electronic coupling, and film formation. Structural and spectroscopic analyses show that the linearly fused NaO1 forms a compact three-dimensional packing network with large and balanced electronic couplings (>24 meV) across multiple directions, while the more curved analogues exhibit excessive crystallization and phase segregation. In-situ optical measurements demonstrate that NaO1 promotes fast and continuous structural evolution during film formation, resulting in smooth morphology and homogeneous phase distribution. These structural and dynamic advantages facilitate efficient charge generation and transport, accompanied by reduced non-radiative energy loss, ultimately achieving an efficiency of 20.07% for non-halogenated ternary devices. Our findings highlight how fused-ring isomerism decisively governs structure–packing–performance relationships in organic solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Relationship between molecular structure and corrugations in self-assembled polypeptoid nanosheets revealed by cryogenic electron microscopy

Designing conformationally dynamic molecules that self-assemble into predictable nanostructures remains an important unmet challenge. This paper describes how atomic-scale cryogenic transmission electron microscopy (cryo-TEM) can be used to explore the relationship between molecular structure and self-assembly of block copolymers. We examined sheetlike micelles formed in water using a series of diblock copolypeptoids with the same hydrophilic block and three distinct crystalline hydrophobic blocks. Our cryo-TEM images revealed all the structures share nansoscale features, but differ in their intermolecular packing geometries. Different molecular arrangements, parallel and antiparallel V-shaped crystal motifs, were revealed by two-dimensional atomic-scale through-plane images. However, images from tilted samples revealed an unexpected feature when the hydrophobic polypeptoid block comprised phenyl rings with substituted bromine atoms at the para position. The nanosheets contained atomic-scale corrugations that were absent in the other systems which comprised unsubstituted aliphatic and aromatic side chains. We hypothesize that these corrugations are due to the dipolar characteristics of the brominated phenyl group and interactions between this group and water molecules. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
NASA NTRS1997

Tensegrity: the architectural basis of cellular mechanotransduction

Physical forces of gravity, hemodynamic stresses, and movement play a critical role in tissue development. Yet, little is known about how cells convert these mechanical signals into a chemical response. This review attempts to place the potential molecular mediators of mechanotransduction (e.g. stretch-sensitive ion channels, signaling molecules, cytoskeleton, integrins) within the context of the structural complexity of living cells. The model presented relies on recent experimental findings, which suggests that cells use tensegrity architecture for their organization. Tensegrity predicts that cells are hard-wired to respond immediately to mechanical stresses transmitted over cell surface receptors that physically couple the cytoskeleton to extracellular matrix (e.g. integrins) or to other cells (cadherins, selectins, CAMs). Many signal transducing molecules that are activated by cell binding to growth factors and extracellular matrix associate with cytoskeletal scaffolds within focal adhesion complexes. Mechanical signals, therefore, may be integrated with other environmental signals and transduced into a biochemical response through force-dependent changes in scaffold geometry or molecular mechanics. Tensegrity also provides a mechanism to focus mechanical energy on molecular transducers and to orchestrate and tune the cellular response.

Review↗
DOE OSTI2025

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗
DOE OSTI2022

Analytic gradients for compressed multistate pair-density functional theory

Photochemical reactions often involve states that are closely coupled due to near degeneracies, for example by proximity to conical intersections. Therefore, a multistate method is used to accurately describe these states; for example, ordinary perturbation theory is replaced by quasidegenerate perturbation theory. Multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to approximate the full dynamical correlation energy of strongly correlated systems, and we recently proposed compressed multistate pair-density functional theory (CMS-PDFT) to treat closely coupled states. In the present paper, we report the implementation of analytic gradients for CMS-PDFT in both OpenMolcas and PySCF, and here we illustrate the use of these gradients by applying the method to the excited states of formaldehyde and phenol.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

X-ray spectroscopic identification of strain and structure-based resonances in a series of saturated carbon-cage molecules: Adamantane, twistane, octahedrane, and cubane

Novel nanocarbons such as fullerenes, nanotubes, graphene, and nanodiamond reside at the cutting edge of nanoscience and technology. Along with chemical functionalization, geometric constraints (such as extreme curvature in nanotubes or defects within or at the surfaces of diamond nanoparticles) significantly alter the electronic states of the nanocarbon material. Understanding the effects of steric strain on the electronic structure is critical to developing nanoelectronic applications based on these materials. This paper presents a fundamental study of how strain affects the electronic structure in a benchmark series of some fundamental saturated carbon cage compounds. Adamantane, C10H16, the smallest diamondoid and arguably the smallest nanodiamond crystallite, has carbon atoms essentially commensurate with diamond lattice positions and possesses by far the least molecular strain of this series. Twistane also is a C10H16 isomer but the fixed cyclohexane twist conformation of the central ring introduces additional strain into the cage. Octahedrane [(CH)12] and cubane [(CH)8] are considerably more strained, culminating in cubane where carbon–carbon bonds lie either parallel or orthogonal to one another. Using gas-phase near-edge x-ray absorption fine structure spectroscopy to probe the unoccupied electronic states, we observe two major progressions across this series. First, a broad C–C σ* resonance in the absorption splits into two more narrow and intense resonances with increasing strain. Second, the first manifold of states previously associated with tertiary C–H σ* in the diamondoid series appears to broaden and shift to lower energy. This feature is more than twice as intense in cubane than in octahedrane, even though these two molecules have only tertiary carbons, with the chemical formula (CH)x. The spectral differences are entirely due to the shape of the molecules; in particular, in cubane, the features arise from a high degree of p-p interaction between parallel C–C bonds. In contrast to the conventional wisdom that near-edge x-ray absorption is primarily an atomically localized spectroscopy, molecular shape and associated strain lead to the dominant features in spectra acquired from this fundamental series of carbon cage structures.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Cross section for Ly-alpha emission by electron impact on methane

Utilizing Lyman-alpha emission cross sections for H2 as secondary standards, absolute values of Lyman-alpha emission cross sections for CH4 have been obtained for electron impact energies varying from threshold to 100 eV. A crossed electron beam-molecular beam geometry was employed and the Lyman-alpha radiation was detected at 90 deg and 45 deg with respect to the incident electron beam by a solar blind photomultiplier in tandem with an oxygen filter. The results are compared with previous measurements. Appreciable differences among the various experimental data are found.

Orient, O. J.↗
DOE OSTI2026

Validated ligand geometries for macromolecular refinement restraints and molecular-mechanics force fields

In macromolecular structure refinement, the low observation-to-parameter ratio and the lack of high-resolution data are countered by using a priori information in the form of restraints. Having accurate geometries of the chemical entities in the sample is paramount for generating accurate chemical restraints and, therefore, accurate macromolecular structures. In particular, it is desirable to have accurate restraints for known and novel ligand entities. Quantum mechanics (QM) can minimize the energy of a ligand by adjusting its geometry, and these geometries can be used to generate restraints for macromolecular refinement. This article describes a library of approximately 37 000 small molecules extracted from the Chemical Component Dictionary in the Protein Data Bank and minimized by density-functional QM. The library includes restraint files for use in crystallography or cryo-EM refinement, along with files suitable for molecular-dynamics simulation. Because the geometries are validated using the Cambridge Structural Database, the restraints library provides users with both functional restraints and minimized geometries. This work also provides procedures for generating new and accurate restraints.

Amber↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Geometry-complete perceptron networks for 3D molecular graphs

Abstract Motivation The field of geometric deep learning has recently had a profound impact on several scientific domains such as protein structure prediction and design, leading to methodological advancements within and outside of the realm of traditional machine learning. Within this spirit, in this work, we introduce GCPNet, a new chirality-aware SE(3)-equivariant graph neural network designed for representation learning of 3D biomolecular graphs. We show that GCPNet, unlike previous representation learning methods for 3D biomolecules, is widely applicable to a variety of invariant or equivariant node-level, edge-level, and graph-level tasks on biomolecular structures while being able to (1) learn important chiral properties of 3D molecules and (2) detect external force fields. Results Across four distinct molecular-geometric tasks, we demonstrate that GCPNet’s predictions (1) for protein–ligand binding affinity achieve a statistically significant correlation of 0.608, more than 5%, greater than current state-of-the-art methods; (2) for protein structure ranking achieve statistically significant target-local and dataset-global correlations of 0.616 and 0.871, respectively; (3) for Newtownian many-body systems modeling achieve a task-averaged mean squared error less than 0.01, more than 15% better than current methods; and (4) for molecular chirality recognition achieve a state-of-the-art prediction accuracy of 98.7%, better than any other machine learning method to date. Availability and implementation The source code, data, and instructions to train new models or reproduce our results are freely available at https://github.com/BioinfoMachineLearning/GCPNet.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

A molecular dynamics study of freezing in a confined geometry

The dynamics of freezing of a Lennard-Jones liquid in narrow channels bounded by molecular walls is studied by computer simulation. The time development of ordering is quantified and a novel freezing mechanism is observed. The liquid forms layers and subsequent in-plane ordering within a layer is accompanied by a sharpening of the layer in the transverse direction. The effects of channel size, the methods of quench, the liquid-wall interaction and the roughness of walls on the freezing mechanism are elucidated. Comparison with recent experiments on freezing in confined geometries is presented.

Ma, Wen-Jong↗
DOE OSTI2026

Electron–Ion Covariance Reveals Geometry-Resolved Inner-Shell Spectra in CS 2 Photodissociation

The chemical shifts of inner-shell atomic orbitals are highly sensitive to chemical bonding and molecular structure. In evolving systems, however, the spectra of distinct molecular species and geometries overlap, obscuring the underlying chemical dynamics. Here we demonstrate the use of electron–ion covariance analysis to combine the structural sensitivity of Coulomb explosion imaging with inner-shell spectroscopy, yielding geometry-specific spectra of transient and product species. We apply this approach to the excited state dynamics of CS 2 probed by ionization above the S 2p edge. Electron–ion covariance with time- and momentum-selected S + and S 2+ ions isolates distinct S 2p photoelectron spectra for ground-state CS 2 , bent photoexcited CS 2 , the CS photoproduct, and bare atomic sulfur ─ species whose spectra overlap strongly in the channel-averaged measurement. Clear chemical shifts are observed in the covariance photoelectron spectrum for each of these species, all of which are consistent with high-level calculations. Here, by extracting the atomic S contribution to the photoelectron spectrum in a finely time-resolved manner we can disentangle this contribution to the overall time-resolved photoelectron spectrum as the photodissociation proceeds. These results demonstrate the promise of electron–ion covariance as a general approach to geometry-resolved inner-shell spectroscopy, opening a route to tracking structural evolution through chemical shifts in complex photoexcited molecules.

Ionization↗
NASA NTRS1997

Molecular Hydrogen in Diffuse Interstellar Clouds of Arbitrary Three-Dimensional Geometry

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen in diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photodissociation rate of H2 in such clouds was computed with a 26 ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H2 line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 micron continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 micron continuum intensity to H I column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be the largest. Using a standard chi square analysis to assess the goodness of fit of our models, we find (at the 60 a level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be much smaller than its extent in the plane of the sky. If D is the distance to G236 + 39, and given standard assumptions about the rate of grain-catalyzed H2 formation, we find that the cloud has an extent along the line of sight that is 0.9 +/- 0.1 times its mean extent projected onto the plane of the sky and a gas density of 53 +/- 8 (100 pc/D) H nuclei/cc and is illuminated by a radiation field of 1.1 +/- 0.2 (100 pc/D) times the mean interstellar radiation field. The derived 100 micron emissivity per nucleon is 1.13 +/- 0.06 x 10(exp -20) MJy/sr sq cm.

Spaans, Marco↗
NASA NTRSDate not supplied

Rational Design of Nanoplasmonic Array Geometries for Biosensing

Background: Molecular diagnostics provide early and accurate diagnosis, which is essential for the prevention and treatment of infectious as well as chronic diseases. These tests are designed to detect disease-specific bioanalytes such as nucleic acid (DNA or RNA) or protein (antigens, antibodies) biomarkers. In the context of infectious disease diagnosis, nucleic acid-based detection methods are known to provide more specific and sensitive results. Here, the presence of a unique sequence belonging to the pathogenic genomic material is targeted to identify species, organism, genera and/or antimicrobial resistant gene markers. The majority of the common nucleic acid based diagnostic techniques require amplification (polymerase chain reaction, isothermal amplification etc.) of the pathogenic genetic material prior to detection impacting diagnostic speed, complexity, and cost thereby limiting ease of use. Thus, the development of simplified nucleic acid-based diagnostics that can be even used in resource-poor settings may hugely benefit patients across the globe. Nanopath is a molecular diagnostics company utilizing a solid-state nanosensor to enable sequence-specific detection of target nucleic acids without the need of amplification. These nanostructures enable ultra-sensitive biomarker detection using geometric, feature-dependent properties highly dependent on the local dielectric environment, allowing them to be sensitive to low concentration binding events. This paper describes an application of this approach to provide highly relevant clinical information within a single doctor’s office visit. Intro: The Nanopath team is in collaboration with NASA (National Aeronautics and Space Administration) and NIST (National Institute of Standards and Technology) to push the bounds of the fundamental physics associated with their biosensing platform. The ability of metals to support electromagnetic surface waves gives rise to surface plasmons when optically illuminated. This property, and its strong sensitivity to changes in the local refractive index, allows for the use of metal nanoparticles as ultra-sensitive transducers. In prior work by members of this team, ensembles of randomly oriented nanoparticles (i.e., colloidal nanorods dispersed on chip) were employed for sequence-specific nucleic acid sensing (1-3). While these particle sensors have the advantage of rapid fabrication, they suffer from low sensitivity and quality factor due to the random particle dispersity. In contrast, in this study we employ ordered array nanoparticle ensembles which can be used to improve sensor sensitivity and figure-of-merit. Study Methods Overview: In this talk, we detail the results of sensing experiments and computational simulations to outline a rational design of the structure of these plasmonic nanoparticle arrays for biomolecular sensing. Through simulation and experiment, we iteratively tailor nanostructure dimension to provide high quality signal and large resonance shifts upon modeled nucleic acid binding. In particular, full-wave electromagnetic simulations were conducted using Lumerical photonic simulation software in which periodic boundary conditions were applied in the x- and y- dimensions for each of the nanoplasmonic sensor geometries. To simulate the resonance response to changes in the bulk solution in contact with the sensor surface, the refractive index of the surrounding media was changed appropriately. Nucleic acid hybridization events were modeled using either using spherical structures approximating the relevant radius of genomic material as estimated by polymer models, or as conformal layers with the known refractive indices for nucleic acids. On the basis of initial simulations, nanosensors were fabricated using traditional electron-beam lithography protocols at NIST. To evaluate consensus between simulations and experiments, bulk sensing experiments were carried out in which the resonance peaks were obtained by submerging the sensors in refractive index standards. Key nanosensor characteristics including resonance peak locations, resonance peak shifts as a function of refractive index, and figure of merit (FOM) of extinction curves were examined between the experimental and simulation results prior to proceeding with simulations on additional geometries and more complex solution conditions, and further device fabrication. This iterative process is repeated toward a rational design of nanoplasmonic array geometries for biosensing optimizing response for targeted disease detection. In summary, this study puts forth a methodology for rational design and characterization of regularly spaced nanoparticle arrays for optics-based biosensing. The results of this study will allow for more informed design of nanostructure geometries towards sequence-specific nucleic acid detection. These improved designs have the potential to improve clinical sensitivity and limit-of-detection across disease indication.

sensor↗