Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse Optimization”

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 1,063 records · Page 59

Optimization of Source Material for in-situ Arsenic Doping via Vapor Transport Deposition of CdTe Films

In-situ doping of CdTe/CdSeTe films with arsenic is performed using pre-doped source material prepared by high pressure Bridgman (HPB) melt growth. Arsenic level in the source material is varied from 1018-1020 cm-3. Correspondingly, dopant incorporation in the films prepared using vapor transport deposition (VTD) varied from 3x1017 cm-3–1x1018cm-3. In this range, dopant activation is found to inversely correlate with the dopant incorporation, suggesting compensation. The results from this study indicate film properties can be adjusted by source material.

41 EE - Solar Energy Technologies Office (EE-4S)↗

SU-E-T-143: Automatic PBS Treatment Planning for Prostate Cancer

Purpose: This study presents a novel idea of using automated planning technique for parallel-opposed prostate plans with pencil beam scanning proton therapy. Methods: A random 5 clinical prostate patients treated with parallel-opposed modulated-scanning proton beam in our proton center were selected in this study. All of the 5 cases were re-planned with in-house developed python scripts integrated into RayStation 4.0 clinical version. This automatic tool was designed to perform patient DICOM data import, external structure contours, CT density assignment, prescription to PTV assignment, beam placement, optimization, and dose calculation. Results: The mean time to generate a complete treatment plan was 1 minute per plan. For the automatically generated plans, 5 out of 5 plans (100%) were deemed clinically acceptable by one experienced radiation oncologist. In addition, the automated plans were overall dosimetrically equivalent to the clinical plans when scored for target coverage, rectum and bladder dose. Mann-Whitney U tests did not show significant dosimetric differences for normal tissue structures between the clinical and automated plans. Conclusion: We have developed a robust and automated method for fully inversed planned paralleled-opposed prostate proton treatment planning. This efficient method can be readily integrated into clinical practice. The tool produces clinically acceptable plans using only manually contoured anatomic structures and simulation CT data set. We anticipate that this tool will not only improve patient access to high-quality proton treatment by simplifying the planning process and will also reduce the effort and cost of incorporating more advanced planning into clinical practice in the near future.

60 APPLIED LIFE SCIENCES↗

Photoredox Catalysis with Spin Magnetic Field Effects: A Scheme for Chiral Resolution

Quantities that break both mirror symmetry and time-reversal symmetry, such as the orbital angular momentum, are known to connect molecular chirality with an applied magnetic field. This concept has led to observations such as magneto-chiral dichroism and chirality-induced spin selectivity (CISS). However, being small, these effects often require additional amplification procedures such as flow chemistry to achieve bulk enantioseparation. In this work, we demonstrate how the magnetic field effect on photogenerated radical pairs, which also breaks time-reversal symmetry, can be harnessed for enantiopurification. Fundamental to this process is the collective decay of the singlet and triplet radical-pair states made possible by an applied magnetic field. Because opposite enantiomers exhibit spin-orbit coupling matrix elements of opposite signs, the singlet and triplet decay channels interfere constructively in one enantiomer. Meanwhile, molecules of the other enantiomer are funnelled into the first enantiomer through excited-state chirality inversion, achieving (dynamic kinetic) chiral resolution. Using an axially chiral binaphthyl derivative and a borane photosensitizer as prototype, we predict an appreciable enantiomeric excess (e.e.) of 90% to be possible at steady state, attained within hundreds of milliseconds when irradiated by a laser. Importantly, our analytical results showcase regimes of perfect enantioselectivity (100% e.e.), accessible by further chemical optimization of the photosensitizer for which general strategies are discussed. Altogether, this work illustrates a so-far untapped but powerful control knob for photoredox catalysis based on spin chemistry principles.

Magnetic properties↗

Conjugation effect of amine molecules in non-aqueous Mg redox flow batteries

Nonaqueous magnesium redox flow batteries (Mg RFBs) are attractive for low-cost, high-energy-density and long-cycle-life stationary energy storage applications. However, state-of-the-art cathode redox-active molecules suffer from low solubility and low redox potential. Herein, we screened a range of cathode redox-active molecules and identified amine molecules as optimal to couple with the Mg anode. The properties of amine derivatives and their performances were collected to establish the correlation between molecular structures and electrochemical performances. The redox potential and solubility of these amine molecules were influenced by the π-conjugated and non-conjugated structures of amine derivatives. Density functional theory (DFT) simulations and the inverse aromatic fluctuation index (FLU −1 ) verified that conjugation had an important role in stabilizing the molecule and increasing its redox potential. Notably, tris[4-(diethylamino)phenyl]amine (TDPA) achieved the highest theoretical energy density (∼120 Wh L −1 ) due to its high solubility (∼0.9 M) and voltage (∼2.5 V vs. Mg/Mg 2+ ). We also demonstrated that ether solvents were crucial for stable, high-solubility catholytes, while bulk anions did not affect the redox potential of these p-type molecules. In a Mg-amine RFB configuration, the battery delivered 2.50 V, a specific discharge capacity of 106.5 mAh g −1 , an initial coulombic efficiency of 90.74%, and a capacity retention of 93.88% after 150 cycles.

Qin, Yunan [University of Utah, Salt Lake City, UT↗

High-throughput spin-bath characterization of spin defects in semiconductors

Detailed knowledge of the local environments of spin defects in semiconductors, such as nitrogenvacancy (NV) centers in diamond or divacancies in silicon carbide, is crucial for optimizing control and entanglement protocols in quantum sensing and information applications. However, at present a direct experimental characterization of individual defect environments is not scalable, as conventional spin-bath measurements are time consuming and difficult to automate. Achieving high-throughput characterization requires short experiments to probe the spin bath. However, with fewer and noisier measurements, the inverse problem of recovering spin-bath properties from measured data becomes ill posed, with multiple spin baths having a high likelihood of yielding the same data. In this work, we present a set of computational tools to resolve the ill-posed inverse problem of recovering the atomic positions and hyperfine couplings of random nuclei surrounding spin defects from sparse, noisy experimental coherence data, which can be obtained in hours. Here, we use a trans-dimensional Bayesian approach that incorporates ab initio data to yield full posterior distributions over nuclear spin environments, enabling robust recovery from limited data. We also provide practical tools and guidelines to determine the limits of detectability for hyperfine couplings under specific dynamical decoupling sequences and sampling conditions. In addition, we demonstrate how the tools developed here, in combination with ab initio simulations of spin baths, can guide the design of efficient experimental protocols for application-specific high-throughput screening. To showcase the utility of our approach, we apply it to design fast dynamical decoupling experiments to characterize the spin baths often individual NV centers in diamond. While the primary focus is on accelerating spin-bath characterization of spin defects, this Bayesian approach also lays the foundation for digital-twin studies of spin defects, where a virtual model of the spin-defect system evolves in real time with ongoing experimental measurements. Together, the set of tools we designed and applied paves the way for scalable deployment of spin defects in semiconductors for quantum sensing and information applications.

Bayesian methods↗

Considerations for Isochronous Data Services over the Proximity-1 Space Link

Future mission concepts for robotic and human explorations will involve a high level of real time control/monitoring operations such as tele-operation for spacecraft rendezvous and surface mobile platforms carrying life-support equipments. The timely dissemination of voice, command, and real-time telemetry for monitoring and coordination purposes is critical for mission success. It is envisioned that future missions will require a network infrastructure capable of supporting isochronous data services. The CCSDS Proximity-1 Space Link Protocol1 could be used to provide isochronous service over the surface-to-Earth relay as well as "beyond-the-horizon" communications between distant Lunar or Mars surface elements. This paper will analyze the latency, jitter, and throughput performance of the Proximity-1 protocol for isochronous applications. In particular we will focus on constrained scenarios where the protocol operates in full-duplex mode, carrying isochronous traffic in one direction and error-controlled traffic in the other direction. We analyze the impact of the strict priority scheme in Proximity-1 on delay jitter and the impact of the isochronous traffic on the efficiency of the reliable data transfer in the other direction, and discuss methods for performance optimization. In general, jitter performance is driving by relative loading of isochronous traffic on the forward link compared to the acknowledgement traffic. Under light loading condition, the upper-bound of the delay jitter is the transmission duration of an acknowledgement frame on the forward link; for higher loading scenarios, the maximum jitter is scaled up by the inverse of the residual bandwidth, i.e., the spare capacity available in the forward link to carry isochronous traffic.

isochronous communications↗

SU-F-J-38: Dose Rates and Preliminary Evaluation of Contouring Similarity Metrics Using 4D Cone Beam CT

Purpose: 4D imaging modalities require detailed characterization for clinical optimization. The On-Board Imager mounted on the linear accelerator was used to investigate dose rates in a tissue mimicking phantom using 4D-CBCT and assess variability of contouring similarity metrics between 4D-CT and 4D-CBCT retrospective reconstructions. Methods: A 125 kVp thoracic protocol was used. A phantom placed on a motion platform simulated a patient’s breathing cycle. An ion chamber was affixed inside the phantom’s tissue mimicking cavities (i.e. bone, lung, and soft tissue). A sinusoidal motion waveform was executed with a five second period and superior-inferior motion. Dose rates were measured at six ion chamber positions. A preliminary workflow for contouring similarity between 4D-CT and 4D-CBCT was established using a single lung SBRT patient’s historical data. Average intensity projection (Ave-IP) and maximum intensity projection (MIP) reconstructions generated offline were compared between the 4D modalities. Similarity metrics included Dice similarity coefficient (DSC), Hausdorff distance, and center of mass (COM) deviation. Two isolated lesions were evaluated in the patient’s scans: one located in the right lower lobe (ITVRLL) and one located in the left lower lobe (ITVLLL). Results: Dose rates ranged from 2.30 (lung) to 5.18 (bone) E-3 cGy/mAs. For fixed acquisition parameters, cumulative dose is inversely proportional to gantry speed. For ITVRLL, DSC were 0.70 and 0.68, Hausdorff distances were 6.11 and 5.69 mm, and COM deviations were 1.24 and 4.77 mm, for Ave-IP and MIP respectively. For ITVLLL, DSC were 0.64 and 0.75, Hausdorff distances were 10.74 and 8.00 mm, and COM deviations were 7.55 and 4.3 mm, for Ave-IP and MIP respectively. Conclusion: While the dosimetric output of 4D-CBCT is low, characterization is necessary to assure clinical optimization. A basic workflow for comparison of simulation and treatment 4D image-based contours was established. This work was partially supported by a Research Scholar Grant (RSG-15-137-01-CCE) from the American Cancer Society.

60 APPLIED LIFE SCIENCES↗

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗

Strontium vapour laser with ionisation and recombination mechanisms of inversion formation

We have experimentally and numerically investigated a multiwavelength strontium vapour laser which simultaneously emits both in the IR spectrum range on the self-terminating transitions in strontium atoms and ions (λ=6.456 and ∼3μm for Sr I, and λ≈1μm for Sr II) and in the visible range on the recombination transition in strontium ions ( λ=0.4305 μm for Sr II). A number of regularities in the behaviour of the output characteristics of lasing, which allow for a targeted search for optimal excitation conditions of the multiwavelength strontium vapour lasers, are established. (lasers, their applications and other topics in quantum electronics)

42 ENGINEERING↗

Optimal parallel solution of sparse triangular systems

A method for the parallel solution of triangular sets of equations is described that is appropriate when there are many right-handed sides. By preprocessing, the method can reduce the number of parallel steps required to solve Lx = b compared to parallel forward or backsolve. Applications are to iterative solvers with triangular preconditioners, to structural analysis, or to power systems applications, where there may be many right-handed sides (not all available a priori). The inverse of L is represented as a product of sparse triangular factors. The problem is to find a factored representation of this inverse of L with the smallest number of factors (or partitions), subject to the requirement that no new nonzero elements be created in the formation of these inverse factors. A method from an earlier reference is shown to solve this problem. This method is improved upon by constructing a permutation of the rows and columns of L that preserves triangularity and allow for the best possible such partition. A number of practical examples and algorithmic details are presented. The parallelism attainable is illustrated by means of elimination trees and clique trees.

Alvarado, Fernando L.↗

A novel approach to interface high-Q Fabry–Pérot resonators with photonic circuits

The unique benefits of Fabry–Pérot resonators as frequency-stable reference cavities and as an efficient interface between atoms and photons make them an indispensable resource for emerging photonic technologies. To bring these performance benefits to next-generation communications, computation, and time-keeping systems, it will be necessary to develop strategies to integrate compact Fabry–Pérot resonators with photonic integrated circuits. In this paper, we demonstrate a novel reflection cancellation circuit that utilizes a numerically optimized multi-port polarization-splitting grating coupler to efficiently interface high-finesse Fabry–Pérot resonators with a silicon photonic circuit. This circuit interface produces a spatial separation of the incident and reflected waves, as required for on-chip Pound–Drever–Hall frequency locking, while also suppressing unwanted back reflections from the Fabry–Pérot resonator. Using inverse design principles, we design and fabricate a polarization-splitting grating coupler that achieves 55% coupling efficiency. This design realizes an insertion loss of 5.8 dB for the circuit interface and more than 9 dB of back reflection suppression, and we demonstrate the versatility of this system by using it to interface several reflective off-chip devices.

Optics↗

Emittance Measurements at FAST Facility

The FAST facility at Fermilab recently been commissioned has demonstrated the generation of electron beam within a wide range of parameter (energy, charge) suitable for accelerator-science and beam-physics experiments. This accelerator consists of a photo-electron gun, injector, ILC-type cryomodules, and multiple downstream beam-lines. It will mainly serve as injector for the upcoming Integrable Optical Test Accelerator (IOTA). At the same time we will also carry out a LINAC based intense gamma ray experiment based on the Inverse Compton scattering. It is essential to understand the beam emittance for both experiments. A number of techniques are used to characaterizing the beam emittance including slit based method and quad scan method. An on-line emittance measurement based on multi-slit method is developed so the emittance measured will be immediately available to support further beam optimization. In this report we will present the results from the emittance studies using this tool. We will also present the emittance measurement based on quads scan technique for the high energy beam line.

43 PARTICLE ACCELERATORS↗

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗

Optimizing the use of remotely sensed atmospheric data through four-dimensional data assimilation

Results with the GLA interactive forecast-retrieval-assimilation system, which directly uses satellite radiance data in the assimilation cycle, are presented. Forecasts made from initial conditions derived with this system are shown to be more skillful than those whose initial conditions came from a comparable assimilation using statistically based temperature soundings, or ones using no sounding data at all. Examples are also given of additional parameters which are produced through the physical inversion of the satellite radiance data during the interactive assimilation cycle. These are potentially useful for analysis purposes.

Pfaendtner, J.↗

Fabrication and accelerated long-term stability test of asymmetrical hollow fiber-supported thin film oxygen separation membrane

Asymmetrical hollow fiber-supported thin film membrane may provide microstructural advantages for air separation and oxygen production. The fabrication of such a membrane is usually very difficult, particularly the sintering behaviors and thermal expansion coefficient (TEC) mismatch among multiple layers. This directly affects the reliability and long-term stability of the membrane. Here, in this research, the sintering behaviors of a set of simple oxides are systematically studied and ZnO is identified as a material component for composite substrate La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF)-ZnO. The LSCF-ZnO ratios are further optimized by trading off several factors, e.g., matching the sintering behaviors and TECs of the LSCF-ZnO composite with those of LSCF thin film separation layer while enhancing the substrate strength. With the identified LSCF-ZnO ratio, the LSCF-ZnO composite hollow fiber substrate precursor is formed through slurry spinning process with porous microstructure being formed via phase inversion process. The thin film LSCF dense separation layer is then dip-coated on the substrate followed by co-sintering process. A thin, porous PrBaCo(Fe 0.6 Zr 0.2 Y 0.2 )O 5+δ (PBCFZY) catalyst layer is fabricated on the dense LSCF separation layer with dip-coating and sintering process, forming an asymmetrical membrane device LSCF-ZnO/LSCF/PBCFZY. Oxygen permeation test of the device is systematically conducted, and the fundamental mechanisms are analyzed. An accelerated long-term stability test of the membrane is also conducted (~550h, 46 thermal cycling loads), demonstrating excellent robustness and durability. The device is characterized and analyzed before and after the test. By replacing a considerable amount of high cost LSCF with low cost ZnO in the substrate, it not only enhances the stability but also reduces the capital cost of the membrane for practical applications.

42 ENGINEERING↗

Less is More: Bigger Data from Compressive Measurements

Compressive sensing approaches are beginning to take hold in (scanning) transmission electron microscopy (S/TEM) [1,2,3]. Compressive sensing is a mathematical theory about acquiring signals in a compressed form (measurements) and the probability of recovering the original signal by solving an inverse problem [4]. The inverse problem is underdetermined (more unknowns than measurements), so it is not obvious that recovery is possible. Compression is achieved by taking inner products of the signal with measurement weight vectors. Both Gaussian random weights and Bernoulli (0,1) random weights form a large class of measurement vectors for which recovery is possible. The measurements can also be designed through an optimization process. The key insight for electron microscopists is that compressive sensing can be used to increase acquisition speed and reduce dose. Building on work initially developed for optical cameras, this new paradigm will allow electron microscopists to solve more problems in the engineering and life sciences. We will be collecting orders of magnitude more data than previously possible. The reason that we will have more data is because we will have increased temporal/spatial/spectral sampling rates, and we will be able ability to interrogate larger classes of samples that were previously too beam sensitive to survive the experiment. For example consider an in-situ experiment that takes 1 minute. With traditional sensing, we might collect 5 images per second for a total of 300 images. With compressive sensing, each of those 300 images can be expanded into 10 more images, making the collection rate 50 images per second, and the decompressed data a total of 3000 images [3]. But, what are the implications, in terms of data, for this new methodology? Acquisition of compressed data will require downstream reconstruction to be useful. The reconstructed data will be much larger than traditional data, we will need space to store the reconstructions during analysis, and the computational demands for analysis will be higher. Moreover, there will be time costs associated with reconstruction. Deep learning [5] is an approach to address these problems. Deep learning is a hierarchical approach to find useful (for a particular task) representations of data. Each layer of the hierarchy is intended to represent higher levels of abstraction. For example, a deep model of faces might have sinusoids, edges and gradients in the first layer; eyes, noses, and mouths in the second layer, and faces in the third layer. There has been significant effort recently in deep learning algorithms for tasks beyond image classification such as compressive reconstruction [6] and image segmentation [7]. A drawback of deep learning, however, is that training the model requires large datasets and dedicated computational resources (to reduce training time to a few days). A second issue is that deep learning is not user-friendly and the meaning behind the results is usually not interpretable. We have shown it is possible to reduce the data set size while maintaining model quality [8] and developed interpretable models for image classification [9], but the demands are still significant. The key to addressing these problems is to NOT reconstruct the data. Instead, we should design computational sensors that give answers to specific problems. A simple version of this idea is compressive classification [10], where the goal is to classify signal type from a small number of compressed measurements. Classification is a much simpler problem than reconstruction, so 1) much fewer measurements will be necessary, and 2) these measurements will probably not be useful for reconstruction. Other simple examples of computational sensing include determining object volume or the number of objects present in the field of view [11].

Stevens, Andrew↗

Aperture‐Synthesis Radar Imaging With Compressive Sensing for Ionospheric Research

Abstract Inverse methods involving compressive sensing are tested in the application of two‐dimensional aperture‐synthesis imaging of radar backscatter from field‐aligned plasma density irregularities in the ionosphere. We consider basis pursuit denoising, implemented with the fast iterative shrinkage thresholding algorithm, and orthogonal matching pursuit (OMP) with a wavelet basis in the evaluation. These methods are compared with two more conventional optimization methods rooted in entropy maximization (MaxENT) and adaptive beamforming (linearly constrained minimum variance or often “Capon's Method.”) Synthetic data corresponding to an extended ionospheric radar target are considered. We find that MaxENT outperforms the other methods in terms of its ability to recover imagery of an extended target with broad dynamic range. Fast iterative shrinkage thresholding algorithm performs reasonably well but does not reproduce the full dynamic range of the target. It is also the most computationally expensive of the methods tested. OMP is very fast computationally but prone to a high degree of clutter in this application. We also point out that the formulation of MaxENT used here is very similar to OMP in some respects, the difference being that the former reconstructs the logarithm of the image rather than the image itself from basis vectors extracted from the observation matrix. MaxENT could in that regard be considered a form of compressive sensing.

Hysell, D. L.↗

Ensuring thermodynamic consistency with invertible coarse-graining

Coarse-grained models are a core computational tool in theoretical chemistry and biophysics. A judicious choice of a coarse-grained model can yield physical insights by isolating the essential degrees of freedom that dictate the thermodynamic properties of a complex, condensed-phase system. The reduced complexity of the model typically leads to lower computational costs and more efficient sampling compared with atomistic models. Designing “good” coarse-grained models is an art. Generally, the mapping from fine-grained configurations to coarse-grained configurations itself is not optimized in any way; instead, the energy function associated with the mapped configurations is. In this work, we explore the consequences of optimizing the coarse-grained representation alongside its potential energy function. We use a graph machine learning framework to embed atomic configurations into a low-dimensional space to produce efficient representations of the original molecular system. Because the representation we obtain is no longer directly interpretable as a real-space representation of the atomic coordinates, we also introduce an inversion process and an associated thermodynamic consistency relation that allows us to rigorously sample fine-grained configurations conditioned on the coarse-grained sampling. We illustrate that this technique is robust, recovering the first two moments of the distribution of several observables in proteins such as chignolin and alanine dipeptide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗