Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

EPSI 2013-2014. Annual Progress Report

Initially a separate orthogonal expansion for the low-k components and FEM for high-k components was considered, or a multi-scale FEM method, but a more attractive hybrid FDM/FEM model for parallel computations has been devised with minimal impact to formerly developed code. This system involves using Lagrange polynomials along magnetic field lines to establish a FDM discretization of the toroidal contributions to the perpendicular Laplacian, which is the operator of importance in gyrokinetic solution of the vector potential. The method was found analytically to be stable, due to some dissipation from varying field-line lengths. The accuracy of this method has been evaluated using the method of manufactured solutions (with Bei Wang of Princeton University), whereby we find that for a Poisson problem the maximum error roughly diminishes linearly as a function of the grid resolution (uniformly refined), whereas the L2 norm as square root of the resolution improvement. Most of the error appears to be contributed near boundaries and where the FEM grid has anisotropy, so we are trying to improve this analysis using improved tessellations obtained from Mark Shepard, and publish the results as soon as the problems with the grid have been eliminated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structure and Potential–Dependent Selectivity in Redox–Metallopolymers: Electrochemically Mediated Multicomponent Metal Separations

Electro-responsive functional materials can play a critical role in selective metal recovery and recycling due to the need for molecular differentiation between transition metals in complex mixtures. Redox-active metallopolymers are a promising platform for electrochemical separations, offering versatile structural tuning and fast electron transfer. First, through a judicious selection of polymer structure between a main-chain metallopolymer (polyferrocenylsilane) and a pendant-group metallopolymer (polyvinylferrocene), charge-transfer interactions and binding strength toward competing metal ions are tuned, which as a result, dictate selectivity. For example, almost an order of magnitude increase in separation factor between chromate and meta-vanadate can be achieved, depending on polymer structure. Second, these metallopolymer electrodes exhibit potential-dependent selectivity that can even flip ion preference, based solely on electrical means—indicating a control parameter that is orthogonal to structural modifications. Finally, this work presents a framework for evaluating electrochemical separations in multicomponent ion mixtures and elucidates the underlying charge-transfer mechanisms resulting in molecular selectivity through a combination of spectroscopy and electronic structure calculations. Furthermore, the findings demonstrate the applicability of redox-metallopolymers in tailored electrochemical separations for environmental remediation, value-added metal recovery, waste recycling, and even mining processing.

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

Effects of Frequency and Motion Paradigm on Perception of Tilt and Translation During Periodic Linear Acceleration

Previous studies have demonstrated an effect of frequency on the gain of tilt and translation perception. Results from different motion paradigms are often combined to extend the stimulus frequency range. For example, Off-Vertical Axis Rotation (OVAR) and Variable Radius Centrifugation (VRC) are useful to test low frequencies of linear acceleration at amplitudes that would require impractical sled lengths. The purpose of this study was to compare roll-tilt and lateral translation motion perception in 12 healthy subjects across four paradigms: OVAR, VRC, sled translation and rotation about an earth-horizontal axis. Subjects were oscillated in darkness at six frequencies from 0.01875 to 0.6 Hz (peak acceleration equivalent to 10 deg, less for sled motion below 0.15 Hz). Subjects verbally described the amplitude of perceived tilt and translation, and used a joystick to indicate the direction of motion. Consistent with previous reports, tilt perception gain decreased as a function of stimulus frequency in the motion paradigms without concordant canal tilt cues (OVAR, VRC and Sled). Translation perception gain was negligible at low stimulus frequencies and increased at higher frequencies. There were no significant differences between the phase of tilt and translation, nor did the phase significantly vary across stimulus frequency. There were differences in perception gain across the different paradigms. Paradigms that included actual tilt stimuli had the larger tilt gains, and paradigms that included actual translation stimuli had larger translation gains. In addition, the frequency at which there was a crossover of tilt and translation gains appeared to vary across motion paradigm between 0.15 and 0.3 Hz. Since the linear acceleration in the head lateral plane was equivalent across paradigms, differences in gain may be attributable to the presence of linear accelerations in orthogonal directions and/or cognitive aspects based on the expected motion paths.

Beaton, K. H.↗
DOE OSTI2017

Solution Deposition of Phenylphosphinic Acid Leads to Highly Ordered, Covalently Bound Monolayers on TiO 2 (110) Without Annealing

Solution-deposited phosphonic acids, O=P(OH) 2 R, have been used to impart new functionality to a variety of metal oxide surfaces for applications ranging from organic field-effect transistors to biocompatible coatings on implants. Interestingly, the as-deposited monolayers are easily rinsed off, becoming robust and strongly adherent only after a long, low-temperature thermal anneal (e.g., 18 h at 120 °C). The need for this thermal treatment has raised questions about the nature of the bonding of the as-deposited monolayer. Is it merely physisorbed, requiring heat treatment for covalent bonding? To understand the first stages of monolayer formation, we have studied the reactivity and molecular bonding geometry of a prototypical, solution-deposited phosphinic acid, O=PH(OH)R, on the prototypical metal oxide surface rutile (110). We show that solution deposition produces near ideal, dense phenylphosphinate monolayers covalently bound in a bridged bidentate geometry. Three nearly orthogonal molecular vibrations—the P–H stretch vibration and the symmetric and antisymmetric OPO stretch vibrations—provided an unambiguous signature of the three-dimensional structure and binding of the adsorbed monolayer; scanning tunneling microscopy provided information on long-range order and intermolecular conformation; and X-ray photoemission spectroscopy provided coverage quantification. Despite their covalent bidentate attachment and significantly higher binding energy than the corresponding carboxylic acid, a H 2 O rinse removed most of the phosphinate monolayer, demonstrating that hydrolytic stability does not result from covalent attachment alone. The H 2 O rinse also oxidized ~25% of the phosphinate monolayer to the corresponding phosphonate, producing a species that was somewhat more resistant to H 2 O rinsing.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Bifurcation to square-wave switching in orthogonally delay-coupled semiconductor lasers: Theory and experiment

We analyze the dynamics of two semiconductor lasers with so-called orthogonal time-delayed mutual coupling: the dominant TE (x) modes of each laser are rotated by 90 deg. (therefore, TM polarization or y) before being coupled to the other laser. Although this laser system allows for steady-state emission in either one or in both polarization modes, it may also exhibit stable time-periodic dynamics including square waveforms. A theoretical mapping of the switching dynamics unveils the region in parameter space where one expects to observe long-term time-periodic mode switching. Detailed numerical simulations illustrate the role played by the coupling strength, the mode frequency detuning, or the mode gain to loss difference. We complement our theoretical study with several experiments and measurements. We present time series and intensity spectra associated with the characteristics of the square waves and other waveforms observed as a function of the strength of the delay coupling. The experimental observations are in very good agreement with the analysis and the numerical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The appearance, motion, and disappearance of three-dimensional magnetic null points

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the magnetic field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2007

Shell Models of RMHD Turbulence and the Heating of Solar Coronal Loops

A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B|| and stratification, nicknamed 'Shell-Atm,' is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfve'n waves propagating along B||. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and haracteristic properties of heating events are distributed as power laws. Crosscorrelations show a delay of dissipation compared with energy content.

corona↗
DOE OSTI2016

The Path to High Q-Factors in Superconducting Accelerating Cavities: Flux Expulsion and Surface Resistance Optimization

Accelerating cavities are devices resonating in the radio-frequency (RF) range used to accelerate charged particles in accelerators. Superconducting accelerating cavities are made out of niobium and operate at the liquid helium temperature. Even if superconducting, these resonating structures have some RF driven surface resistance that causes power dissipation. In order to decrease as much as possible the power losses, the cavity quality factor must be increased by decreasing the surface resistance. In this dissertation, the RF surface resistance is analyzed for a large variety of cavities made with different state-of-the-art surface treatments, with the goal of finding the surface treatment capable to return the highest Q-factor values in a cryomodule-like environment. This study analyzes not only the superconducting properties described by the BCS surface resistance, which is the contribution that takes into account dissipation due to quasi-particle excitations, but also the increasing of the surface resistance due to trapped flux. When cavities are cooled down below their critical temperature inside a cryomodule, there is always some remnant magnetic field that may be trapped increasing the global RF surface resistance. This thesis also analyzes how the fraction of external magnetic field, which is actually trapped in the cavity during the cooldown, can be minimized. This study is performed on an elliptical single-cell horizontally cooled cavity, resembling the geometry of cavities cooled in accelerator cryomodules. The horizontal cooldown study reveals that, as in case of the vertical cooldown, when the cooling is performed fast, large thermal gradients are created along the cavity helping magnetic flux expulsion. However, for this geometry the complete magnetic flux expulsion from the cavity equator is more difficult to achieve. This becomes even more challenging in presence of orthogonal magnetic field, that is easily trapped on top of the cavity equator causing temperature rising. The physics behind the magnetic flux expulsion is also analyzed, showing that during a fast cooldown the magnetic field structures, called vortices, tend to move in the same direction of the thermal gradient, from the Meissner state region to the mixed state region, minimizing the Gibbs free energy. On the other hand, during a slow cool down, not only the vortices movement is limited by the absence of thermal gradients, but, also, at the end of the superconducting transition, the magnetic field concentrates along randomly distributed normal-conducting region from which it cannot be expelled anymore. The systematic study of the surface resistance components performed for the different surface treatments, reveals that the BCS surface resistance and the trapped flux surface resistance have opposite trends as a function of the surface impurity content, defined by the mean free path. At medium field value, the BCS surface resistance is minimized for nitrogen-doped cavities and significantly larger for standard niobium cavities. On the other hand, Nitrogen-doped cavities show larger dissipation due to trapped flux. This is consequence of the bell-shaped trend of the trapped flux sensitivity as a function of the mean free path. Such experimental findings allow also a better understanding of the RF dissipation due to trapped flux. The best compromise between all the surface resistance components, taking into account the possibility of trapping some external magnetic field, is given by light nitrogen-doping treatments. However, the beneficial effects of the nitrogen-doping is completely lost when large amount of magnetic field is trapped during the cooldown, underlying the importance of both cooldown and magnetic field shielding optimization in high quality factors cryomodules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

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

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

36 MATERIALS SCIENCE↗
NASA NTRS2000

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.↗
NASA NTRS2007

Fitting a Two-Component Scattering Model to Polarimetric SAR Data from Forests

Two simple scattering mechanisms are fitted to polarimetric synthetic aperture radar (SAR) observations of forests. The mechanisms are canopy scatter from a reciprocal medium with azimuthal symmetry and a ground scatter term that can represent double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants or Bragg scatter from a moderately rough surface, which is seen through a layer of vertically oriented scatterers. The model is shown to represent the behavior of polarimetric backscatter from a tropical forest and two temperate forest sites by applying it to data from the National Aeronautic and Space Agency/Jet Propulsion Laboratory's Airborne SAR (AIRSAR) system. Scattering contributions from the two basic scattering mechanisms are estimated for clusters of pixels in polarimetric SAR images. The solution involves the estimation of four parameters from four separate equations. This model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem. The model is used to develop an understanding of the ground-trunk double-bounce scattering that is present in the data, which is seen to vary considerably as a function of incidence angle. Two parameters in the model fit appear to exhibit sensitivity to vegetation canopy structure, which is worth further exploration. Results from the model fit for the ground scattering term are compared with estimates from a forward model and shown to be in good agreement. The behavior of the scattering from the ground-trunk interaction is consistent with the presence of a pseudo-Brewster angle effect for the air-trunk scattering interface. If the Brewster angle is known, it is possible to directly estimate the real part of the dielectric constant of the trunks, a key variable in forward modeling of backscatter from forests. It is also shown how, with a priori knowledge of the forest height, an estimate for the attenuation coefficient of the canopy can be obtained directly from the multi-incidence-angle polarimetric observations. This attenuation coefficient is another key variable in forward models and is generally related to the canopy density.

temperate forests↗
NASA NTRS2017

JIGSAW-GEO (1.0): Locally Orthogonal Staggered Unstructured Grid Generation for General Circulation Modelling on the Sphere

An algorithm for the generation of non-uniform, locally orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate very high-quality staggered VoronoiDelaunay meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of high-quality unstructured spheroidal Delaunay triangulations is introduced. A locally orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Frontal-Delaunay refinement technique allows for the generation of very high-quality unstructured triangulations, satisfying a priori bounds on element size and shape. Grid quality is further improved through the application of hill-climbing-type optimisation techniques. Overall, the algorithm is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. A selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling are presented. These grids are shown to satisfy the geometric constraints associated with contemporary unstructured C-grid-type finite-volume models, including the Model for Prediction Across Scales (MPAS-O). The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution-type studies is discussed in detail.

geophysical fluid dynamics↗
DOE OSTI2014

PAPR-Constrained Pareto-Optimal Waveform Design for OFDM-STAP Radar

We propose a peak-to-average power ratio (PAPR) constrained Pareto-optimal waveform design approach for an orthogonal frequency division multiplexing (OFDM) radar signal to detect a target using the space-time adaptive processing (STAP) technique. The use of an OFDM signal does not only increase the frequency diversity of our system, but also enables us to adaptively design the OFDM coefficients in order to further improve the system performance. First, we develop a parametric OFDM-STAP measurement model by considering the effects of signaldependent clutter and colored noise. Then, we observe that the resulting STAP-performance can be improved by maximizing the output signal-to-interference-plus-noise ratio (SINR) with respect to the signal parameters. However, in practical scenarios, the computation of output SINR depends on the estimated values of the spatial and temporal frequencies and target scattering responses. Therefore, we formulate a PAPR-constrained multi-objective optimization (MOO) problem to design the OFDM spectral parameters by simultaneously optimizing four objective functions: maximizing the output SINR, minimizing two separate Cramer-Rao bounds (CRBs) on the normalized spatial and temporal frequencies, and minimizing the trace of CRB matrix on the target scattering coefficients estimations. We present several numerical examples to demonstrate the achieved performance improvement due to the adaptive waveform design.

Sen, Satyabrata↗
NASA NTRS1989

Analysis of radar images of the active volcanic zone at Krafla, Iceland: The effects of look azimuth biasing

The geomorphic expression of Mid-Ocean-Ridge (MOR) volcanism in a subaerial setting occurs uniquely on Earth in Iceland, and the most recent MOR eruptive activity has been concentrated in the Northeastern Volcanic Zone in an area known as Krafla. Within the Krafla region are many of the key morphologic elements of MOR-related basaltic volcanism, as well as volcanic explosion craters, subglacial lava shields, tectonic fissure swarms known as gjar, and basaltic-andesite flows with well developed ogives (pressure-ridges). The objective was to quantify the degree to which the basic volcanic and structural features can be mapped from directional SAR imagery as a function of the look azimuth. To accomplish this, the current expression of volcanic and tectonic constructs was independently mapped within the Krafla region on the E, W, and N-looking SAR images, as well as from SPOT Panchromatic imagery acquired in 1987. The initial observations of the E, W, and N images indicates that fresh a'a lava surfaces are extremely radar bright (rough at 3 cm to meter scales) independent of look direction; this suggests that these flows do not have strong flow direction related structures at meter and cm scales, which is consistent with typical Icelandic a'a lava surfaces in general. The basic impression from a preliminary analysis of the effects of look azimuth biasing on interpretation of the geology of an active MOR volcanic zone is that up to 30 percent of the diagnostic features can be missed at any given look direction, but that having two orthogonal look direction images is probably sufficient to prevent gross misinterpretation.

Garvin, J. B.↗
NASA NTRSDate not supplied

Mapping Aerosol Lidar Ratios Over Ocean for CALIPSO Using Modis AOD-Constrained Retrievals and A Global Aerosol Model

The CALIPSO aerosol algorithms currently assign one lidar ratio value for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms by developing regional and seasonal lidar ratio climatologies. In this study, aerosol lidar ratios are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) datasets. The analysis is subsampled for only those profiles that are cloud-free and contain one aerosol type (based on CALIOP feature classification). In addition, the CALIOP profiles are collocated with aerosol volume fractions simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. In this talk, the twelve-year (2006-2017) mean spatial distributions of inferred aerosol lidar ratios for CALIOP-classified marine aerosols and the corresponding modeled sea salt volume fractions (SSVF) will be shown. Model-assisted climatological lidar ratio maps on a seasonal scale will also be presented, developed from the empirical relationship found between the modeled SSVF and lidar ratios. A comparison with past studies will be provided as well as results of a sensitivity study regarding the variability of retrieved lidar ratios as a function of horizontal averaging resolution. While the majority of this talk will focus on lidar ratios for CALIOP-classified marine aerosols, preliminary results will be shown for other aerosol types over ocean, such as dust and elevated smoke.

Travis Toth↗
DOE OSTI2023

Beam dynamics corrections in the measurement of the anomalous precession frequency at the Muon $g-2$ experiment at Fermilab

The Muon $g-2$ experiment at Fermilab (E989) aims to measure the anomalous magnetic moment of the muon with an accuracy of 140 ppb (parts per billions). This accuracy, obtained by adding in quadrature a statistical and a systematic contribution of comparable value (100 ppb), will allow to reduce the experimental uncertainty (from the previous E821 experiment at BNL \cite{bnl}) of a factor of 4, and represents one of the most precise tests of the Standard Model (SM) theory of elementary particles. The first result on the Run-1 dataset \cite{prl} was released on April 7, 2021, showing a very good agreement with the previous result from BNL experiment, with a slightly better uncertainty. The corresponding experimental average increases the significance of the discrepancy between the measured and Standard Model prediction of 4.2$\sigma$ \cite{white_paper}.\\ \noindent The measured quantity is the muon magnetic anomaly a$_\mu$=$\frac{g_\mu-2}{2}$ where g$_\mu$ is the gyromagnetic factor of the muon. Dirac's equation predicts g$_\mu$ = 2, while radiative corrections, dominated by the QED contribution due to an exchange of a virtual photon, causes a per-mille correction on this quantity. By including all the SM contributions, a$_\mu$ is known at 370 ppb. The E989 experiment measures $a_\mu$ injecting positive muons with momentum of 3.1 GeV/$c$ polarized longitudinally in bunches (called $fills$) with an average rate of 12 Hz, in a storage ring of 14 meters diameter. Due to the parity violation in the weak muon decay, high energy positrons produced are emitted preferably in the muon's spin direction. By counting the number of positrons with energy greater then 1.7 GeV in function of the time, the frequency precession of the muon spin is measured, that together with the measurement of the magnetic field, allows to extract $a_\mu$. The positrons are detected with 24 electromagnetic calorimeters, that measure the energy and the arrival time of the positrons, each made of 54 crystals of lead fluoride (PbF$_2$) read by silicon photomultipliers (SiPM). Together with calorimeters, two tracking detectors are used to make non-destructive measurements of the muon beam characteristic by reconstructing the muon decay position extrapolating backward the decay positrons. The knowledge of the beam motion inside the ring plays a fundamental role in the analysis of $a_{\mu}$, where the measured anomalous precession frequency $\omega_a$ must be corrected for four main beam dynamics effects. A first correction is associated with the presence of an electric field responsible of the vertical focusing on the storage ring, where vertical direction is orthogonal to the orbit and horizontal direction is along the storage ring radius. Due to the oscillations in the vertical direction, the so called vertical betatron oscillations, a second correction is necessary to account for an average angle associated with the muon motion off the ideal orbit. A third correction is caused by lost muons in the ring which have a different spin phase at the injection respect to the decay ones. Finally, due to the correlation between the vertical and horizontal beam motion of the muons and the acceptance of the calorimeter, a correction (called ``phase-acceptance") arises. In Run1, due to the presence of two damaged resistors in one electrostatic quadrupole, this correction had a prominent role. The replacement of the damaged resistors before Run2 decreased this effect by one order of magnitude, and a further reduction in Run-3 was provided by an improved orbit.\\ \noindent The work of this Thesis focuses on the beam dynamics corrections on $\omega_a$. Due to the prominent role in Run-1 a special attention has been put to the phase-acceptance correction\footnote{The effect of this correction on the E821 BNL measurement of the $g-2$ has been evaluated to be within the quoted systematic error.}. Chapter \ref{ch:anomalous_magnetic_moment} introduces the anomalous magnetic moment of the muon. Chapter \ref{ch:early_experiment} describes the history of the Muon $g-2$ experiments. Chapter \ref{ch:theory} discusses the Standard Model prediction and possible new physics scenario. Chapter \ref{ch:muon_g2_experiment} describes the E989 experiment, whose experimental technique and the apparatus are discussed focusing on the improvements needed to reach the final goal on $a_\mu$ measurement. The original work of the Thesis is discussed in the last 5 chapters. Chapter \ref{ch:experiment_technique} presents the analysis technique to extract $\omega_a$, Chapter \ref{ch:beam_dynamics} describes the beam dynamics in the E989 experiment, Chapter \ref{ch:beam_dynamics_correction} discusses the beam dynamics correction to $\omega_a$, and Chapter \ref{sec:pa_corr} focuses on the phase acceptance correction. Finally Chapter \ref{ch:run23_analysis_improvements} contains the Run-2/3 improvements which are expected to allow for an increased precision ons the beam dynamics co...

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Identifying Important Ions and Positions in Mass Spectrometry Imaging Data Using CUR Matrix Decompositions

Mass spectrometry imaging enables label-free, high-resolution spatial mapping of the chemical composition of complex, biological samples. Typical experiments require selecting ions and/or positions from the images: ions for fragmentation studies to identify keystone compounds and positions for follow up validation measurements using microdissection or other orthogonal techniques. Unfortunately, with modern imaging machines, these must be selected from an overwhelming amount of raw data. Existing techniques to reduce the volume of data, the most popular of which are principle component analysis and non-negative matrix factorization, have the disadvantage that they return difficult-to-interpret linear combinations of actual data elements. In this work, we show that CX and CUR matrix decompositions can be used directly to address this selection need. CX and CUR matrix decompositions use empirical statistical leverage scores of the input data to provide provably good low-rank approximations of the measured data that are expressed in terms of actual ions and actual positions, as opposed to difficult-to-interpret eigenions and eigenpositions. We show that this leads to effective prioritization of information for both ions and positions. In particular, important ions can be found either by using the leverage scores as a ranking function and using a deterministic greedy selection algorithm or by using the leverage scores as an importance sampling distribution and using a random sampling algorithm; however, selection of important positions from the original matrix performed significantly better when they were chosen with the random sampling algorithm. Also, we show that 20 ions or 40 locations can be used to reconstruct the original matrix to a tolerance of 17% error for a widely studied image of brain lipids; and we provide a scalable implementation of this method that is applicable for analysis of the raw data where there are often more than a million rows and/or columns, which is larger than SVD-based low-rank approximation methods can handle. Lastly, these results introduce the concept of CX/CUR matrix factorizations to mass spectrometry imaging, describing their utility and illustrating principled algorithmic approaches to deal with the overwhelming amount of data generated by modern mass spectrometry imaging.

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

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex

How early sensory experience during “critical periods” of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent [S. Chenget al.,Cell185, 311–327.e24 (2022)]. Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial–ventricular axis in V1. By applying multitasking theory, we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types. By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell-type-specific gene expression programs. By contrast, the second program impacts L2/3 cell-type identity, regulating a subset of cell-type-specific genes and shifting the distribution of cells within the L2/3 cell-type manifold. Through an integrated analysis of spatial transcriptomics with single-nucleus RNA-seq data, we describe how vision patterns cortical L2/3 cell types during the critical period.

Science & Technology - Other Topics↗