Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiability”

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,603 records · Page 89

Stable calcium-free myocilin olfactomedin domain variants reveal challenges in differentiating between benign and glaucoma-causing mutations

Nonsynonymous gene mutations can be beneficial, neutral, or detrimental to the stability, structure, and biological function of the encoded protein, but the effects of these mutations are often not readily predictable. For example, the β-propeller olfactomedin domain of myocilin (mOLF) exhibits a complex interrelationship among structure(s), stability, and aggregation. Numerous mutations within mOLF are linked to glaucoma; the resulting variants are less stable, aggregation-prone, and sequestered intracellularly, causing cytotoxicity. In this study, we report the first stable mOLF variants carrying substitutions in the calcium-binding site that exhibit solution characteristics indistinguishable from those of glaucoma variants. Crystal structures of these stable variants at 1.8–2.0-Å resolution revealed features that we could not predict by molecular dynamics simulations, including loss of loop structure, helix unwinding, and a blade shift. Double mutants that combined a stabilizing substitution and a selected glaucoma-causing single-point mutant rescued in vitro folding and stability defects. In the context of full-length myocilin, secretion of stable single variants was indistinguishable from that of the WT protein, and the double mutants were secreted to varying extents. In summary, our finding that mOLF can tolerate particular substitutions that render the protein stable despite a conformational switch emphasizes the complexities in differentiating between benign and glaucoma-causing variants and provides new insight into the possible biological function of myocilin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes

Current Monte Carlo codes use one of three models: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S(α,β) model, depending on the neutron energy and the specific Monte Carlo code. This thesis addresses the consequences of using the free gas scattering model, which assumes that the neutron interacts with atoms in thermal motion in a monatomic gas in thermal equilibrium at material temperature, T. Most importantly, the free gas model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that the exact resonance scattering model is temperaturedependent, and neglecting the resonances in the lower epithermal range can under-predict resonance absorption due to the upscattering phenomenon mentioned above, leading to an over-prediction of keff by several hundred pcm. Existing methods to address this issue involve changing the neutron weights or implementing an extra rejection scheme in the free gas sampling scheme, and these all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame to continue the random walk of the neutron. The goal of this paper was to develop a sampling methodology that (1) accounted for the energydependent scattering cross sections in the collision analysis and (2) was performed in the laboratory frame,avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials (2nd and 4th order) to approximate the scattering cross section in Blackshaw’s equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using methods developed in this dissertation showed very close comparison to results using the reference Dopplerbroadened rejection correction (DBRC) scheme.

42 ENGINEERING↗

Measurements of Neutron Scattering from a Copper Sample Using a Quasi-Differential Method in the Region from 2 keV to 20 MeV

Neutron scattering from a copper sample was measured at Rensselaer Polytechnic Institute utilizing the quasi-differential method. The measurement spanned the energy range from 0.5 to 20 MeV using the high-energy scattering system and from 2 keV to 0.5 MeV using the new mid-energy scattering system. Copper was selected as a material of interest to measure due to large discrepancies between experiments and simulations of the Zeus benchmark. The Zeus benchmark consists of a copper reflected highly enriched uranium system, and the angular distribution of copper scattering was thought to potentially be the cause of the discrepancy. The copper measurements found differences in the scattering response particularly in the incident energy region from 1 to 2 MeV for the high-energy measurement and from 2 to 4 keV in the mid-energy system. These differences are particularly noticeable at angles near 90 deg in the high-energy system and back angles in the mid-energy system. Additionally, for ENDF/B-VIII.0 there is a large discrepancy at the forward angle in the energy range around 0.5 MeV. For these reasons, a new evaluation of copper scattering utilizing these results is recommended and perhaps could help to improve the agreement with the Zeus benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Element differentiation with a Hartmann- based X-ray phase imaging system

Significant efforts are currently ongoing in X-ray imaging to provide multimodal imaging systems, targeting better sensitivity and specificity for both biomedical or non-destructive testing applications. Knowing the elemental composition of specific structures, such as breast microcalcifications in mammography, would help to differentiate malign and benign tumours. Standard X-ray Phase Contrast Imaging techniques provide only qualitative information on elements with similar absorption properties. However, their chemical composition can be determined from the measurement of the phase as it is directly related to the optical index of elemental materials. We will show new experimental results obtained with an X-ray phase imaging system based on a Hartmann mask. Early data treatment succeeded in retrieving both the real and imaginary parts of the refractive index. The system demonstrates the capability to discriminate materials based on elemental composition.

36 MATERIALS SCIENCE↗

Negative differential mobility and trapping in active matter systems

Using simulations, we explore the average velocity as a function of applied drift force for active matter particles moving through a random obstacle array. We find that for low drift force, there is an initial flow regime where the mobility increases linearly with drive, while for higher drift forces a regime of negative differential mobility appears in which the velocity decreases with increasing drive due to the trapping of active particles behind obstacles. A fully clogged regime exists at very high drift forces when all the particles are permanently trapped behind obstacles. We discover for increasing activity that the overall mobility is nonmonotonic, with an enhancement of the mobility for small levels of activity and a decrease in mobility for large activity levels. We show how these effects evolve as a function of disk and obstacle density, active run length, drift force, and motor force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗

Harnessing mass differential confinement effects in magnetized rotating plasmas to address new separation needs

The ability to separate large volumes of mixed species based on atomic mass appears desirable for a variety of emerging applications with high societal impact. One possibility to meet this objective consists in leveraging mass differential effects in rotating plasmas. Beyond conventional centrifugation, rotating plasmas offer in principle additional ways to separate elements based on mass. Single ion orbits show that ion radial mass separation in a uniform magnetized plasma column can be achieved by applying a tailored electric potential profile across the column, or by driving a rotating magnetic field within the column. Furthermore, magnetic pressure and centrifugal effects can be combined in a non-uniform geometry to separate ions based on mass along the field lines. In conclusion, practical application of these separation schemes hinges on the ability to produce the desirable electric and magnetic field configuration within the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗

Revisiting dark energy models using differential ages of galaxies

In this work, we use a test based on the differential ages of galaxies for distinguishing the dark energy models. As proposed by Jimenez and Loeb in [1], relative ages of galaxies can be used to put constraints on various cosmological parameters. In the same vein, we reconstruct H {sub 0} {sub dt} / dz and its derivative ( H {sub 0} {sub d} {sup 2} {sup t} / dz {sup 2}) using a model independent technique called non-parametric smoothing . Basically, dt / dz is the change in the age of the object as a function of redshift which is directly linked with the Hubble parameter. Hence for reconstruction of this quantity, we use the most recent H ( z ) data. Further, we calculate H {sub 0} {sub dt} / dz and its derivative for several models like Phantom, Einstein de Sitter (EdS), ΛCDM, Chevallier-Polarski-Linder (CPL) parametrization, Jassal-Bagla-Padmanabhan (JBP) parametrization and Feng-Shen-Li-Li (FSLL) parametrization. We check the consistency of these models with the results of reconstruction obtained in a model independent way from the data. It is observed that H {sub 0} {sub dt} / dz as a tool is not able to distinguish between the ΛCDM, CPL, JBP and FSLL parametrizations but, as expected, EdS and Phantom models show noticeable deviation from the reconstructed results. Further, the derivative of H {sub 0} {sub dt} / dz for various dark energy models is more sensitive at low redshift. It is found that the FSLL model is not consistent with the reconstructed results, however, the ΛCDM model is in concordance with the 3σ region of the reconstruction at redshift z ≥ 0.3.

79 ASTRONOMY AND ASTROPHYSICS↗

CRPropa 3.1—a low energy extension based on stochastic differential equations

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel κ{sub ∥} and perpendicular κ{sub ⊥} diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

79 ASTRONOMY AND ASTROPHYSICS↗

The nature of the progression of drought stress drives differential metabolomic responses in Populus deltoides

Background and Aims: The use of woody crops for Quad-level (approx. 1 × 10 18 J) energy production will require marginal agricultural lands that experience recurrent periods of water stress. Populus species have the capacity to increase dehydration tolerance by lowering osmotic potential via osmotic adjustment. The aim of this study was to investigate how the inherent genetic potential of a Populus clone to respond to drought interacts with the nature of the drought to determine the degree of biochemical response. Methods: A greenhouse drought stress study was conducted on Populus deltoides ‘WV94’ and the resulting metabolite profiles of leaves were determined by gas chromatography–mass spectrometry following trimethylsilylation for plants subjected to cyclic mild (–0.5 MPa pre-dawn leaf water potential) drought vs. cyclic severe (–1.26 MPa) drought in contrast to well-watered controls (–0.1 MPa) after two or four drought cycles, and in contrast to plants subjected to acute drought, where plants were desiccated for up to 8 d.Key Results: The nature of drought (cyclic vs. acute), frequency of drought (number of cycles) and the severity of drought (mild vs. severe) all dictated the degree of osmotic adjustment and the nature of the organic solutes that accumulated. Whereas cyclic drought induced the largest responses in primary metabolism (soluble sugars, organic acids and amino acids), acute onset of prolonged drought induced the greatest osmotic adjustment and largest responses in secondary metabolism, especially populosides (hydroxycinnamic acid conjugates of salicin). Conclusions: The differential adaptive metabolite responses in cyclic vs. acute drought suggest that stress acclimation occurs via primary metabolism in response to cyclic drought, whereas expanded metabolic plasticity occurs via secondary metabolism following severe, acute drought. Here, the shift in carbon partitioning to aromatic metabolism with the production of a diverse suite of higher order salicylates lowers osmotic potential and increases the probability of post-stress recovery.

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of growth and development of sorghum genotypes with differential susceptibility to Striga hermonthica

Two sorghum varieties, Shanqui Red (SQR) and SRN39, have distinct levels of susceptibility to the parasitic weed Striga hermonthica, which have been attributed to different strigolactone composition within their root exudates. Root exudates of the Striga-susceptible variety Shanqui Red (SQR) contain primarily 5-deoxystrigol, which has a high efficiency for inducing Striga germination. SRN39 roots primarily exude orobanchol, leading to reduced Striga germination and making this variety resistant to Striga. The structural diversity in exuded strigolactones is determined by a polymorphism in the LOW GERMINATION STIMULANT 1 (LGS1) locus. Yet, the genetic diversity between SQR and SRN39 is broad and has not been addressed in terms of growth and development. Here, we demonstrate additional differences between SQR and SRN39 by phenotypic and molecular characterization. A suite of genes related to metabolism was differentially expressed between SQR and SRN39. Increased levels of gibberellin precursors in SRN39 were accompanied by slower growth rate and developmental delay and we observed an overall increased SRN39 biomass. The slow-down in growth and differences in transcriptome profiles of SRN39 were strongly associated with plant age. Additionally, enhanced lateral root growth was observed in SRN39 and three additional genotypes exuding primarily orobanchol. In summary, we demonstrate that the differences between SQR and SRN39 reach further than the changes in strigolactone profile in the root exudate and translate into alterations in growth and development.

09 BIOMASS FUELS↗

A differentiable model of the evolution of dark matter halo concentration

ABSTRACT We introduce a new model of the evolution of the concentration of dark matter haloes, c(t). For individual haloes, our model approximates c(t) as a power law with a time-dependent index, such that at early times, concentration has a nearly constant value of c ≈ 3–4, and as cosmic time progresses, c(t) smoothly increases. Using large samples of halo merger trees taken from the Bolshoi–Planck and MultiDark Planck 2 cosmological simulations, we demonstrate that our three-parameter model can approximate the evolution of the concentration of individual haloes with a typical accuracy of 0.1 dex for $t\gtrsim 2\, {\rm Gyr}$ for all Bolshoi–Planck and MultiDark Planck 2 haloes of present-day peak mass $M_{0}\gtrsim 10^{11.5}\, {\rm M}_{\odot }$. We additionally present a new model of the evolution of the concentration of halo populations, which we show faithfully reproduces both average concentration growth and the diversity of smooth trajectories of c(t), including capturing correlations with halo mass and halo assembly history. Our publicly available source code, diffprof, can be used to generate Monte Carlo realizations of the concentration histories of cosmologically representative halo populations. diffprof is differentiable due to its implementation in the jax autodiff library, which facilitates the incorporation of our model into existing analytical halo model frameworks.

Stevanovich, Dash↗

A differentiable perturbation-based weak lensing shear estimator

Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-per cent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (fpfs), using jax, which is publicly available as part of a new suite of analytic shear algorithms called anacal. This code can analytically calibrate the shear response of any non-linear observable constructed with the fpfs shapelets and detection modes utilizing autodifferentiation (ad), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the ad capability of jax, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new anacal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the 10 yr LSST. For the generic shear estimator, the magnitude of the multiplicative bias |m| is below 3 × 10 –3 (99.7 per cent confidence interval), and the effective galaxy number density is improved by 5 per cent. We also discuss some planned future additions to the anacal software suite to extend its applicability beyond the fpfs measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ProteoSign: an end-user online differential proteomics statistical analysis platform

Profiling of proteome dynamics is crucial for understanding cellular behavior in response to intrinsic and extrinsic stimuli and maintenance of homeostasis. Over the last 20 years, mass spectrometry (MS) has emerged as the most powerful tool for large-scale identification and characterization of proteins. Bottom-up proteomics, the most common MS-based proteomics approach, has always been challenging in terms of data management, processing, analysis and visualization, with modern instruments capable of producing several gigabytes of data out of a single experiment. Here, we present ProteoSign, a freely available web application, dedicated in allowing users to perform proteomics differential expression/abundance analysis in a userfriendly and self-explanatory way. Although several non-commercial standalone tools have been developed for post-quantification statistical analysis of proteomics data, most of them are not end-user appealing as they often require very stringent installation of programming environments, third-party software packages and sometimes further scripting or computer programming. To avoid this bottleneck, we have developed a user-friendly software platform accessible via a web interface in order to enable proteomics laboratories and core facilities to statistically analyse quantitative proteomics data sets in a resource-efficient manner. ProteoSign is available at http://bioinformatics.med.uoc.gr/ProteoSign and the source code at https://github.com/yorgodillo/ ProteoSign.

59 BASIC BIOLOGICAL SCIENCES↗

Chronic drought differentially alters the belowground microbiome of drought tolerant and drought susceptible genotypes of Populus trichocarpa

Populus trichocarpa is an ecologically important tree species and economically important biofeedstock. Belowground, P. trichocarpa interacts with diverse microorganisms in the rhizosphere and root endosphere. These plant-microbial interactions can bolster a variety of plant processes, ranging from pathogen suppression to drought tolerance, yet we know little about the impact of chronic drought stress on P. trichocarpa’s belowground microbiomes. To investigate the interactive effect of chronic drought on belowground microbial communities across genetically different P. trichocarpa hosts, we assessed archaeal/bacterial and fungal communities within the root endosphere, rhizosphere, and surrounding bulk soil of selected genotypes in a long-term drought experiment in Boardman, OR, USA. We sequenced the 16S rRNA and ITS2 gene region on samples collected from 16 distinct P. trichocarpa genotypes in plots with full or reduced irrigation. Eight of these genotypes have been previously identified as drought tolerant while the other eight genotypes were drought susceptible. While reduced irrigation influenced the composition of every archaeal/bacterial microbiome compartment, fungal communities were only affected in the rhizosphere and bulk soil compartments. Drought-tolerant bacteria, such as Actinobacteria, were differentially abundant in reduced irrigation across all belowground microbiomes. Host drought-tolerance influenced plant-associated microbiome compartments but had little impact on the bulk soil compartment. Drought-tolerant trees were enriched for potential growth-promoting microorganisms in the root endosphere and rhizosphere, including Sphingomonas bacteria and ectomycorrhizal fungi. Overall, associations of growth-promoting microbes in drought resistant P. trichocarpa genotypes can be leveraged to improve biofeedstock productivity in regions prone to periodic drought.

59 BASIC BIOLOGICAL SCIENCES↗

Solving generalized eigenvalue problems by ordinary differential equations on a quantum computer

Many eigenvalue problems arising in practice are often of the generalized form A x = λ B x . One particularly important case is symmetric, namely A , B are Hermitian and B is positive definite. The standard algorithm for solving this class of eigenvalue problems is to reduce them to Hermitian eigenvalue problems. For a quantum computer, quantum phase estimation is a useful technique to solve Hermitian eigenvalue problems. In this work, we propose a new quantum algorithm for symmetric generalized eigenvalue problems using ordinary differential equations. The algorithm has lower complexity than the standard one based on quantum phase estimation. Moreover, it works for a wider case than symmetric: B is invertible, B - 1 A is diagonalizable and all the eigenvalues are real.

Science & Technology - Other Topics↗

Measurement of the liquid argon scintillation pulse shape using differentiable simulation in the coherent CAPTAIN-Mills experiment

The coherent CAPTAIN-Mills (CCM) experiment is a liquid argon (LAr) light collection detector searching for MeV-scale neutrino and beyond Standard Model physics signatures. Two hundred eight-inch photomultiplier tubes instrument the seven-ton fiducial volume with 50% photocathode coverage to detect light produced by charged particles. CCM’s light-based approach reduces requirements of LAr purity, compared to other detection technologies, such that sub-MeV particles can be reliably detected without additional LAr filtration and with 𝒪⁡(1) parts per million of common contaminants. We present a measurement of LAr light production and propagation parameters, with uncertainties, obtained from a sample of MeV-scale electromagnetic events. The optimization of this high-dimensional parameter space was facilitated by a differentiable optical photon Monte Carlo simulation and detailed photomultiplier tube response characterization. This result accurately predicts the timing and spatial distribution of light due to scintillation and Cherenkov emission in the detector. This is the first description of photon propagation in LAr to include several effects, including anomalous dispersion of the index of refraction near the ultraviolet resonance, Mie scattering from impurities, and Cherenkov light production.

Cherenkov detectors↗