SEARCH · Engineering Papers
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.
Our research focuses on designing all-reflective phase retarders (RPRs) to generate and maintain circularly polarized (CP) light for NSF OPAL and other multipetawatt laser facilities. Assuming that the input polarization is either s- or p-polarized, RPRs need to be used in an out-of-plane configuration. Here, we will discuss the impact of transporting and focusing optics on CP light.
Dynamical Lie algebras (DLAs) are a versatile tool for various topics that span from the expressibility-trainability of variational quantum algorithms (VQAs), to simulation of many body Hamiltonians. Quantum gates and most of the Hamiltonians of interest consist of local interactions; therefore, the analysis of all possible DLAs generated by 1- and 2-local operators is crucial for quantum simulation and VQAs on current hardware. Previously in [R. Wiersema et al ., npj Quantum Inf. 10 , 110 (2024)], we analyzed the DLAs on linear, circular and all-to-all topologies, and obtained results about their dimensions and algebraic structure. Here, in this work, we extend our analysis into any possible hardware topology and provide a classification of all DLAs generated by Pauli strings on any undirected interaction graph. Our results indicate that the DLAs depend solely on whether the connectivity or interaction graph is bipartite or not. In addition, we find that the non-trivial polynomially scaling DLAs appear only on 1D line or circle topologies, and all other DLAs have dimensions scaling exponentially with the system size. Together with the current VQA literature, our results imply that either the majority of VQAs are non-trainable, or we are yet to understand the role of DLAs on the trainability of VQAs.
Nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) films offer tremendous potential for diverse electronic applications. However, the absence of a reliable Ohmic contact at room temperature limits their practical integration and broader applicability. Here, in this study, we investigate the room temperature specific contact resistivity (ρ c ) of Ti/Pt/Au metal stack deposited on n-UNCD films grown on an intrinsic single crystal diamond substrate using a microwave plasma chemical vapor deposition system. We employ a circular transfer length model (c-TLM) and find the room temperature ρ c to be ∼4.67 × 10 −5 Ω cm 2 , which is among the lowest reported value for n-UNCD films. High temperature vacuum annealing conducted at 700 and 800 °C results in an initial improvement, followed by a minor degradation in ρ c values, respectively. The electrical contacts remain highly Ohmic for all measurements. Furthermore, cross-sectional transmission electron microscopy analysis suggests formation of conductive titanium carbide layer with no significant inter metallic diffusion. Overall, the electrical contacts demonstrate robust thermal stability, both of which are critical for attaining high-performance nanocrystalline diamond-based electronic devices.
Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.
Protein structure prediction algorithms, such as AlphaFold, have accelerated protein design and advanced the understanding of the relationship between amino acid sequence and protein structure. However, these algorithms are limited in their ability to predict the structures of conformationally dynamic, intrinsically disordered, and stimuli-responsive proteins. To evaluate sequence-to-structure predictions of such challenging proteins, we explored a class of conformationally dynamic, repeats-in-toxin (RTX) proteins. RTX proteins adopt intrinsically disordered conformations in the absence of calcium and undergo reversible folding into β-roll structures upon binding to calcium. RTX proteins are characterized by tandem repeats of the sequence GGXGXDXUX, in which X can be any amino acid and U is an aliphatic amino acid. We designed RTX sequence variants with global substitutions of nonconserved amino acids, tandem repeats of consensus sequences GGAGXDTLY, and tandem repeats of scrambled sequences GGAGXDTYL. AlphaFold2 and AlphaFold3 predicted that all of these RTX variants adopt β-roll structures, characteristic of wild-type RTX bound to calcium. However, modeling the predicted structures with molecular dynamics simulations and characterizing the protein variants with circular dichroism spectroscopy, small-angle x-ray scattering, and x-ray crystallography revealed that variants adopt diverse, sequence-dependent structures in the absence and presence of calcium. To better design proteins for applications in biotechnology and sustainability, it is critical to build predictive tools that consider intrinsically disordered protein states and validate these tools with multi-mode, multi-scale experimental data.
Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.
Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.
We present images of microholes drilled in diamond using a homemade femtosecond Yb : KGW laser. We use a femtosecond laser source emitting pulses with a duration of 230 fs at a wavelength of 1030 nm, the focusing spot size amounting to 8.9 μm. The effect of the pulse energy and the number of pulses on the microhole geometry (hole diameter, circularity, taper angle, and drilling quality) is evaluated. The obtained results demonstrate the feasibility of drilling of hourglass-shaped holes in a diamond sample, which have similar diameters at the hole entrance (92 μm) and exit (95 μm), but a much smaller diameter (28 μm) at a certain waist section inside the hole. (paper)
Based on the model of a bipolarised laser, which takes into account real orientations of absorbing and emitting dipoles of active centres in the unit cell of an Nd : YAG single crystal, we have shown that the appearance of low-frequency relaxation oscillations requires the participation of both competing pump channels through linearly and circularly polarised dipole transitions. (paper)
We have numerically simulated the process of measuring stress – strain states by the method of speckle-shearing interferometry using the phase-shift technique. A computer model with the possibility of setting its strain and roughness is developed, which includes a model of a diffusely reflecting test object corresponding to the characteristics of a real membrane, as well as a speckle interferometer model that allows speckle interferograms to be obtained for different speckle sizes and angles between interfering beams. The process of reconstructing the object surface topogram from model speckle interferograms by the phase-shift technique is implemented. Using the developed models, a two-dimensional shearogram are obtained, which is a derivative of the strain field of a circular membrane. Comparison of the results of numerical simulation with experimental data shows that the differences (rms deviations) do not exceed 0.02 μm. It is also shown that the error of interferogram reconstruction by the phase-shift technique increases significantly when the test object strains exceed 12 μm. (paper)
We consider the effect of ion mobility on the threshold of relativistic self-induced transparency under irradiation of thin plasma layers by circularly polarised laser radiation. An analytical model of the motion of ions during the removal of electrons from the layer surface by the ponderomotive force of laser radiation is constructed. The model is used to analyse the motion of probe electrons in the resulting electromagnetic field. It is shown that the higher the ion mobility and the longer the laser pulse, the more stable the plasma layer to longitudinal perturbations of a finite value and the higher the threshold of self-induced transparency. These conclusions are verified by one-dimensional numerical simulation of a complete system of kinetic equations for plasma and Maxwell’s equations. (paper)
The ability to precisely design large proteins with diverse shapes would enable applications ranging from the design of protein binders that wrap around their target to the positioning of multiple functional sites in specified orientations. We describe a protein backbone design method for generating a wide range of rigid fusions between helix-containing proteins and use it to design 75,000 structurally unique junctions between monomeric and homo-oligomeric de novo designed and ankyrin repeat proteins (RPs). Of the junction designs that were experimentally characterized, 82% have circular dichroism and solution small-angle X-ray scattering profiles consistent with the design models and are stable at 95 °C. Crystal structures of four designed junctions were in close agreement with the design models with rmsds ranging from 0.9 to 1.6 Å. Electron microscopic images of extended tetrameric structures and ~10-nm-diameter “L” and “V” shapes generated using the junctions are close to the design models, demonstrating the control the rigid junctions provide for protein shape sculpting over multiple nanometer length scales.
The properties of organic molecules can be influenced by magnetic fields, and these magnetic field effects are diverse. They range from inducing nuclear Zeeman splitting for structural determination in NMR spectroscopy to polaron Zeeman splitting organic spintronics and organic magnetoresistance. A pervasive magnetic field effect on an aromatic molecule is the aromatic ring current, which can be thought of as an induction of a circular current of π-electrons upon the application of a magnetic field perpendicular to the π-system of the molecule. While in NMR spectroscopy the effects of ring currents on the chemical shifts of nearby protons are relatively well understood, and even predictable, the consequences of these modified electronic states on the spectroscopy of molecules has remained unknown. In this work, we find that photophysical properties of model phthalocyanine compounds and their aggregates display clear magnetic field dependences up to 25 T, with the aggregates showing more drastic magnetic field sensitivities depending on the intermolecular interactions with the amplification of ring currents in stacked aggregates. These observations are consistent with ring currents measured in NMR spectroscopy and simulated in time-dependent density functional theory calculations of magnetic field-dependent phthalocyanine monomer and dimer absorption spectra. Here, we propose that ring currents in organic semiconductors, which commonly comprise aromatic moieties, may present new opportunities for the understanding and exploitation of combined optical, electronic, and magnetic properties.
Significance Reinforcement learning (RL) has been applied effectively in games and robotic manipulation. We demonstrate the effectiveness of RL in experimental fluid mechanics by applying it to reduce the drag of circular cylinders in turbulent flow, a canonical fluid–structure interaction problem. Although physics agnostic, RL managed to reduce the drag by 30 % or reach another specified optimum point very quickly. Following this discovery, we used high-fidelity simulations to probe the underlying physical mechanisms so that the discovered control techniques can be generalized to other similar flow problems. More broadly, RL-guided active control can lead to efficient exploration of additional flow-control strategies in experimental fluid mechanics, potentially paving the way for accelerating scientific discovery and different designs in flow-related engineering problems.
Locomotion of an organism interacting with an environment is the consequence of a symmetry-breaking action in space-time. Here we show a minimal instantiation of this principle using a thin circular sheet, actuated symmetrically by a pneumatic source, using pressure to change shape nonlinearly via a spontaneous buckling instability. This leads to a polarized, bilaterally symmetric cone that can walk on land and swim in water. In either mode of locomotion, the emergence of shape asymmetry in the sheet leads to an asymmetric interaction with the environment that generates movement––via anisotropic friction on land, and via directed inertial forces in water. Scaling laws for the speed of the sheet of the actuator as a function of its size, shape, and the frequency of actuation are consistent with our observations. The presence of easily controllable reversible modes of buckling deformation further allows for a change in the direction of locomotion in open arenas and the ability to squeeze through confined environments––both of which we demonstrate using simple experiments. Our simple approach of harnessing elastic instabilities in soft structures to drive locomotion enables the design of novel shape-changing robots and other bioinspired machines at multiple scales.
Significance More than 400 million tons of plastic waste is produced each year, the overwhelming majority of which ends up in landfills. Bioconversion strategies aimed at plastics have emerged as important components of enabling a circular economy for synthetic plastics, especially those that exhibit chemically similar linkages to those found in nature, such as polyesters. The enzyme system described in this work is essential for mineralization of the xenobiotic components of poly(ethylene terephthalate) (PET) in the biosphere. Our description of its structure and substrate preferences lays the groundwork for in vivo or ex vivo engineering of this system for PET upcycling.
In many organs, small openings across capillary endothelial cells (ECs) allow the diffusion of low–molecular weight compounds and small proteins between the blood and tissue spaces. These openings contain a diaphragm composed of radially arranged fibers, and current evidence suggests that a single-span type II transmembrane protein, plasmalemma vesicle-associated protein-1 (PLVAP), constitutes these fibers. Here, we present the three-dimensional crystal structure of an 89-amino acid segment of the PLVAP extracellular domain (ECD) and show that it adopts a parallel dimeric alpha-helical coiled-coil configuration with five interchain disulfide bonds. The structure was solved using single-wavelength anomalous diffraction from sulfur-containing residues (sulfur SAD) to generate phase information. Biochemical and circular dichroism (CD) experiments show that a second PLVAP ECD segment also has a parallel dimeric alpha-helical configuration—presumably a coiled coil—held together with interchain disulfide bonds. Overall, ~2/3 of the ~390 amino acids within the PLVAP ECD adopt a helical configuration, as determined by CD. We also determined the sequence and epitope of MECA-32, an anti-PLVAP antibody. Taken together, these data lend strong support to the model of capillary diaphragms formulated by Tse and Stan in which approximately ten PLVAP dimers are arranged within each 60- to 80-nm-diameter opening like the spokes of a bicycle wheel. Passage of molecules through the wedge-shaped pores is presumably determined both by the length of PLVAP—i.e., the long dimension of the pore—and by the chemical properties of amino acid side chains and N-linked glycans on the solvent-accessible faces of PLVAP.
Acylated peptides composed of glucagon-like peptide-1 receptor agonists modified with a fatty acid side chain are an important class of therapeutics for type 2 diabetes and obesity but are susceptible to an unusual physical instability in the presence of hydrophobic surfaces, i.e., spontaneous emulsification, also known as ouzo formation in practice. In this work, light scattering, small-angle X-ray scattering, and circular dichroism measurements are used to characterize the physical properties of the semaglutide colloidal phase, including size distribution, shape, secondary structure, internal structure, and internal composition, as a function of solution physico-chemical conditions. The existence and size of the colloids formed are successfully predicted by a classical Rayleigh model, which identifies the parameters controlling their size and formation. Colloid formation is found to be catalyzed by hydrophobic surfaces, and formation rates are modeled as an autocatalytic reaction, enabling the formation of a master curve for various surfaces that elucidates the mechanism. Surfaces differ due to differences in surface wettability, which can be correlated with Hansen solubility parameters. This work provides insights into this unusual colloidal phenomenon and guides the peptide synthesis process and drug product formulation in the pharmaceutical industry.