Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conformal prediction”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural switching dynamically controls the doubly pseudoknotted Rous sarcoma virus–programmed ribosomal frameshifting element

A hallmark of retrovirus replication is the translation of two different polyproteins from one RNA through programmed –1 frameshifting. This is a mechanism in which the actively translating ribosome is induced to slip in the 5′ direction at a defined codon and then continues translating in the new reading frame. Programmed frameshifting controls the stoichiometry of viral proteins and is therefore under stringent evolutionary selection. Forty years ago, the first frameshifting stimulatory element was discovered in the Rous sarcoma virus. The ~120 nt RNA segment was predicted to contain a pseudoknot, but its 3D structure has remained elusive. Now, we have determined cryoEM and X-ray crystallographic structures of this classic retroviral element, finding that it adopts a butterfly-like double-pseudoknot fold. One “wing” contains a dynamic pyrimidine-rich helix, observed crystallographically in two conformations and in a third conformation via cryoEM. The other wing encompasses the predicted pseudoknot, which interacts with a second unexpected pseudoknot through a toggle residue, A2546. This key purine switches conformations between structural states and tunes the stability of interacting residues in the two wings. We find that its mutation can modulate frameshifting by as much as 50-fold, likely by altering the relative abundance of different structural states in the conformational ensemble of the RNA. Taken together, our structure–function analyses reveal how a dynamic double pseudoknot junction stimulates frameshifting by taking advantage of conformational heterogeneity, supporting a multistate model in which high Shannon entropy enhances frameshifting efficiency.

Science & Technology - Other Topics↗
DOE OSTI2020

Polarization whorls from M87* at the event horizon telescope

The event horizon telescope (EHT) is anticipated to soon produce polarimetric images of the supermassive black hole at the centre of the neighbouring galaxy M87. There are indications that this black hole is rapidly spinning. General relativity predicts that such a high-spin black hole has an emergent conformal symmetry near its event horizon. In this paper, we use this symmetry to analytically predict the polarized near-horizon emissions to be seen at the EHT and find a distinctive pattern of whorls aligned with the spin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA’s conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg 2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg 2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Modeling receptor flexibility in the structure-based design of KRAS G12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS G12C —sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS G12C . In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS G12C , we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Structural basis of peptidoglycan endopeptidase regulation

Most bacteria surround themselves with a cell wall, a strong meshwork consisting primarily of the polymerized aminosugar peptidoglycan (PG). PG is essential for structural maintenance of bacterial cells, and thus for viability. PG is also constantly synthesized and turned over; the latter process is mediated by PG cleavage enzymes, for example, the endopeptidases (EPs). EPs themselves are essential for growth but also promote lethal cell wall degradation after exposure to antibiotics that inhibit PG synthases (e.g., β-lactams). Thus, EPs are attractive targets for novel antibiotics and their adjuvants. However, we have a poor understanding of how these enzymes are regulated in vivo, depriving us of novel pathways for the development of such antibiotics. Here, we have solved crystal structures of the LysM/M23 family peptidase ShyA, the primary EP of the cholera pathogen Vibrio cholerae. Our data suggest that ShyA assumes two drastically different conformations: a more open form that allows for substrate binding and a closed form, which we predicted to be catalytically inactive. Mutations expected to promote the open conformation caused enhanced activity in vitro and in vivo, and these results were recapitulated in EPs from the divergent pathogens Neisseria gonorrheae and Escherichia coli. Our results suggest that LysM/M23 EPs are regulated via release of the inhibitory Domain 1 from the M23 active site, likely through conformational rearrangement in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

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

Experimental Evidence of Dioxole Unimolecular Decay Pathway for Isoprene-Derived Criegee Intermediates

Ozonolysis of isoprene, one of the most abundant volatile organic compounds emitted into the Earth’s atmosphere, generates two four-carbon unsaturated Criegee intermediates, methyl vinyl ketone oxide (MVK-oxide) and methacrolein oxide (MACR-oxide). The extended conjugation between the vinyl substituent and carbonyl oxide groups of these Criegee intermediates facilitates rapid electrocyclic ring closures that form 5-membered cyclic peroxides, known as dioxoles. This research discusses the first experimental evidence of this novel decay pathway, which is predicted to be the dominant atmospheric sink for specific conformational forms of MVK-oxide (anti) and MACR-oxide (syn) with the vinyl substituent adjacent to the terminal O atom. The resulting dioxoles are predicted to undergo rapid unimolecular decay to oxygenated hydrocarbon radical products, including acetyl, vinoxy, formyl, and 2-methyl-vinoxy radicals. In the presence of O 2 , these radicals rapidly react to form peroxy radicals (ROO), which quickly decay via carbon-centered radical intermediates (QOOH) to stable carbonyl products that are identified in this work. The carbonyl products are detected under thermal conditions (298 K, 10 torr He) using multiplexed photoionization mass spectrometry (MPIMS). The main products (and associated relative abundances) originating from unimolecular decay of anti-MVK-oxide and subsequent reaction with O 2 are formaldehyde (88 ± 5%), ketene (9 ± 1%) and glyoxal (3 ± 1%). Those identified from the unimolecular decay of syn-MACR-oxide and subsequent reaction of O 2 are acetaldehyde (37 ± 7%), vinyl alcohol (9 ± 1%), methylketene (2 ± 1%), and acrolein (52 ± 5%). In addition to the stable carbonyl products, the secondary peroxy chemistry also generates OH or HO 2 radical co-products.

09 BIOMASS FUELS↗
DOE OSTI2025

Unraveling trace anomaly of supradense matter via neutron star compactness scaling

The trace anomaly Δ ≡ 1/3 −𝑃/𝜖 =1/3 −𝜙 quantifies the possibly broken conformal symmetry in supradense matter under pressure 𝑃 at energy density 𝜖. Perturbative QCD (pQCD) predicts a vanishing Δ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) is still very uncertain. The extraction of Δ from NS observations strongly depends on the employed model for nuclear equation of state (EOS). Using the IPAD-TOV method based on an intrinsic and perturbative analysis of the dimensionless (IPAD) Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using 10 5 EOSs generated randomly with a metamodel in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness 𝜉 ≡ 𝐺⁡𝑀 NS /𝑅⁢𝑐 2 ≡ 𝑀 NS /𝑅 of a NS with mass 𝑀 NS and radius 𝑅 scales very accurately with $\bar{Π}$ c ≡ $Π$ c · (1 +18⁢X/25) ≡ X/(1 +3⁢X 2 +4⁢X) · (1 +18⁢X/25) where X ≡ 𝜙 c = 𝑃 c /𝜖 c is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly Δ c = 1/3 −X directly from the observational data of either the mass-radius or red-shift measurements. Finally, we then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model insensitively the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter.

nuclear astrophysics↗
DOE OSTI2026

SmileyLlama: modifying large language models for directed chemical space exploration

Here we show that large language models (LLMs) can be transformed via supervised fine-tuning of engineered prompts into SmileyLlama for exploring the chemical space of drug molecules. We benchmark SmileyLlama against pretrained LLMs and chemical language models trained from scratch for generating valid and novel drug-like molecules, and use direct preference optimization to both improve SmileyLlama’s adherence to a prompt and as part of the iMiner reinforcement learning framework to predict molecules with optimized three-dimensional conformations and high binding affinity to drug targets. By training an LLM to speak directly as a chemical language model, while retaining most of its natural language capabilities, we show that SmileyLlama can reliably generate molecules with user-specified properties rather than acting only as a chatbot with knowledge of chemistry or as a virtual assistant. While SmileyLlama is geared toward drug discovery, the supervised fine-tuning/direct preference optimization/LLM framework can be extended to other chemical, biological and materials applications.

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

Multipartite entanglement in the one-dimensional spin- 1 2 Heisenberg antiferromagnet

Multipartite entanglement refers to the simultaneous entanglement between multiple subsystems of a many-body quantum system. While multipartite entanglement can be difficult to quantify analytically, it is known that it can be witnessed through the quantum Fisher information (QFI), a quantity that can also be related to dynamical Kubo response functions. In this work, we first show that the finite temperature QFI can generally be expressed in terms of a static structure factor of the system, plus a correction that vanishes as T→0. This implies that the static structure factor witnesses multipartite entanglement near quantum critical points at temperatures below a characteristic energy scale of the system. Therefore, in systems with a known static structure factor, we can deduce finite temperature scaling of multipartite entanglement and low temperature entanglement depth without knowledge of the full dynamical response function of the system. This is particularly useful to study 1D quantum critical systems in which sub-power-law divergences can dominate entanglement growth, where the conventional scaling theory of the QFI breaks down. Furthermore, the 1D spin- 1 2 antiferromagnetic Heisenberg model is an important example of such a system, and we show that multipartite entanglement in the Heisenberg chain diverges nontrivially as ~ln(1/T) 3/2 . We verify these predictions with calculations of the QFI using conformal field theory and matrix product state simulations. Finally, we discuss the implications of our results for experiments to probe entanglement in quantum materials, comparing to neutron scattering data in KCuF 3 , a material well described by the Heisenberg chain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effects of HLA single chain trimer design on peptide presentation and stability

MHC class I “single-chain trimer” molecules, coupling MHC heavy chain, β 2 -microglobulin, and a specific peptide into a single polypeptide chain, are widely used in research. To more fully understand caveats associated with this design that may affect its use for basic and translational studies, we evaluated a set of engineered single-chain trimers with combinations of stabilizing mutations across eight different classical and non-classical human class I alleles with 44 different peptides, including a novel human/murine chimeric design. While, overall, single-chain trimers accurately recapitulate native molecules, care was needed in selecting designs for studying peptides longer or shorter than 9-mers, as single-chain trimer design could affect peptide conformation. In the process, we observed that predictions of peptide binding were often discordant with experiment and that yields and stabilities varied widely with construct design. We also developed novel reagents to improve the crystallizability of these proteins and confirmed novel modes of peptide presentation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1985

The application of diffusion theory to the analysis of hydrogen desorption data at 25 deg C

The application of diffusion theory to the analysis of hydrogen desorption data (coulombs of H2 desorbed versus time) has been studied. From these analyses, important information concerning hydrogen solubilities and the nature of the hydrogen distributions in the metal has been obtained. Two nickel base alloys, Rene' 41 and Waspaloy, and one ferrous alloy, 4340 steel, are studied in this work. For the nickel base alloys, it is found that the hydrogen distributions after electrolytic charging conforms closely to those which would be predicted by diffusion theory. For Waspaloy samples charged at 5,000 psi, it is found that the hydrogen distributions are essentially the same as those obtained by electrolytic charging. The hydrogen distributions in electrolytically charged 4340 steel, on the other hand, are essentially uniform in nature, which would not be predicted by diffusion theory. A possible explanation has been proposed. Finally, it is found that the hydrogen desorption is completely explained by the nature of the hydrogen distribution in the metal, and that the fast hydrogen is not due to surface and sub-surface hydride formation, as was originally proposed.

Danford, M. D.↗
NASA NTRS1986

New vistas in the determination of hydrogen in aerospace engine metal alloys

The application of diffusion theory to the analysis of hydrogen desorption data has been studied. From these analyses, important information concerning hydrogen solubilities and the nature of the hydrogen distributions in the metal has been obtained. Two nickel base alloys, Rene' 41 and Waspaloy, and one ferrous alloy, 4340 steel, were studied in this work. For the nickel base alloys, it was found that the hydrogen distributions after electrolytic charging conformed closely to those which would be predicted by diffusion theory. The hydrogen distributions in electrolytically charged 4340 steel, on the other hand, were essentially uniform in nature, which would not be predicted by diffusion theory. Finally, it has been found that the hydrogen desorption is completely explained by the nature of the hydrogen distribution in the metal, and that the 'fast' hydrogen is not due to surface and subsurface hydride formation, as was originally proposed.

Danford, M. D.↗
NASA NTRS1994

Equilibrium and Dynamics Properties of Poly(oxyethylene) Melts and Related Poly(alkylethers) from Simulations and Ab Initio Calculations

Molecular dynamics simulations of POE melts have been performed utilizing a potential force field parameterized to reproduce conformer energies and rotational energy barriers in dimethoxyethane as determined from ab initio electronic structure calculations. Chain conformations and dimensions of POE from the simulations were found to be in good agreement with predictions of a rotational isomeric state (RIS) model based upon the ab initio conformational. energies. The melt chains were found to be somewhat extended relative to chains at theta conditions. This effect will be discussed in light of neutron scattering experiments which indicate that POE chains are extended in the melt relative to theta solutions. The conformational characteristics of POE chains will also be compared with those of other poly (alkylethers), namely poly(oxymethylene), poly(oxytrimethylene) and poly(oxytetramethylene). Local conformational dynamics were found to be more rapid than in polymethylene. Calculated C-H vector correlation times were found to be in reasonable agreement with experimental values from C-13 NMR spin-lattice relaxation times. The influence of ionic salts on local conformations and dynamics will also be discussed.

Smith, Grant D.↗
DOE OSTI2022

Rapid Allylic 1,6 H-Atom Transfer in an Unsaturated Criegee Intermediate

Here, a novel allylic 1,6 hydrogen-atom-transfer mechanism is established through infrared activation of the 2-butenal oxide Criegee intermediate, resulting in very rapid unimolecular decay to hydroxyl (OH) radical products. A new precursor, Z/E -1,3-diiodobut-1-ene, is synthesized and photolyzed in the presence of oxygen to generate anew four-carbon Criegee intermediate with extended conjugation across the vinyl and carbonyl oxide groups that facilitates rapid allylic 1,6 H-atom transfer. A low-energy reaction pathway involving isomerization of 2-butenal oxide from a lower-energy ( tZZ ) conformer to a higher-energy ( cZZ ) conformer followed by 1,6 hydrogen transfer via a seven-membered ring transition state is predicted theoretically and shown experimentally to yield OH products. The low-lying ( tZZ ) conformer of 2-butenaloxide is identified based on computed an harmonic frequencies and intensities of its conformers. Experimental IR action spectra recorded in the fundamental CH stretch region with OH product detection by UV laser-induced fluorescence reveal a distinctive IR transition of the low-lying ( tZZ ) conformer at 2996 cm -1 that results in rapid unimolecular decay to OH products. Statistical RRKM calculations involving a combination of conformational isomerization and unimolecular decay via 1,6 H-transfer yield an effective decay rate k eff ( E ) on the order of 10 8 s -1 at ca. 3000 cm -1 in good accord with the experiment. Unimolecular decay proceeds with significant enhancement due to quantum mechanical tunneling. A rapid thermal decay rate of ca. 10 6 s -1 is predicted by master-equation modeling of 2-butenal oxide at 298 K, 1 bar. This novel unimolecular decay pathway is expected to increase the nonphotolytic production of OH radicals upon alkene ozonolysis in the troposphere.

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

Wind Prediction Accuracy for Air Traffic Management Decision Support Tools

The performance of Air Traffic Management and flight deck decision support tools depends in large part on the accuracy of the supporting 4D trajectory predictions. This is particularly relevant to conflict prediction and active advisories for the resolution of conflicts and the conformance with of traffic-flow management flow-rate constraints (e.g., arrival metering / required time of arrival). Flight test results have indicated that wind prediction errors may represent the largest source of trajectory prediction error. The tests also discovered relatively large errors (e.g., greater than 20 knots), existing in pockets of space and time critical to ATM DST performance (one or more sectors, greater than 20 minutes), are inadequately represented by the classic RMS aggregate prediction-accuracy studies of the past. To facilitate the identification and reduction of DST-critical wind-prediction errors, NASA has lead a collaborative research and development activity with MIT Lincoln Laboratories and the Forecast Systems Lab of the National Oceanographic and Atmospheric Administration (NOAA). This activity, begun in 1996, has focussed on the development of key metrics for ATM DST performance, assessment of wind-prediction skill for state of the art systems, and development/validation of system enhancements to improve skill. A 13 month study was conducted for the Denver Center airspace in 1997. Two complementary wind-prediction systems were analyzed and compared to the forecast performance of the then standard 60 km Rapid Update Cycle - version 1 (RUC-1). One system, developed by NOAA, was the prototype 40-km RUC-2 that became operational at NCEP in 1999. RUC-2 introduced a faster cycle (1 hr vs. 3 hr) and improved mesoscale physics. The second system, Augmented Winds (AW), is a prototype en route wind application developed by MITLL based on the Integrated Terminal Wind System (ITWS). AW is run at a local facility (Center) level, and updates RUC predictions based on an optimal interpolation of the latest ACARS reports since the RUC run. This paper presents an overview of the study's results including the identification and use of new large mor wind-prediction accuracy metrics that are key to ATM DST performance.

Cole, Rod↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Mechanosensitive Ion Channels in Bacteria: Functional Domains and Mechanisms of Gating

The past funding period was productive for the group. The progress in the mechanosensitive channel field was critically affected in the end of 1998 by the solution of the crystal structure of the mycobacterial homolog of MscL by our colleagues from Caltech. Having the structure of TbMscL in the closed state, we developed a detailed homology model of EcoMscL, and related the structural model with the wealth of functional phenomenology available for the E. coli version of the channel (EcoMscL). The biophysical properties of the open MscL helped to model the open conformation and infer the pathway for the entire gating transition. The following experiments provided strong support to the atomic model of the gating process, and allowed to make further predictions. The work has advanced our understanding of tension-driven conformational transitions in membrane-embedded mechanosensory proteins, determine major energetic contributions and set the stage for further exploration of the whole family of mechanosensitive channels. The results have been published in seven experimental and theoretical papers, with three other papers currently in press or in preparation.

Sukharev, Sergei↗