Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Distinguishing Mesoscale Polar Order (Unidirectional vs Bidirectional) of Cellulose Microfibrils in Plant Cell Walls Using Sum Frequency Generation Spectroscopy

Cellulose in plant cell walls are synthesized as crystalline microfibrils with diameters of 3–4 nm and lengths of around 1–10 μm. These microfibrils are known to be the backbone of cell walls, and their multiscale three-dimensional organization plays a critical role in cell wall functions including plant growth and recalcitrance to degradation. The mesoscale organization of microfibrils over a 1–100 nm range in cell walls is challenging to resolve because most characterization techniques investigating this length scale suffer from low spatial resolution, sample preparation artifacts, or inaccessibility of specific cell types. Here, we report a sum frequency generation (SFG) study determining the mesoscale polarity of cellulose microfibrils in intact plant cell walls. SFG is a nonlinear optical spectroscopy technique sensitive to the molecular-to-mesoscale order of noncentrosymmetric domains in amorphous matrices. However, the quantitative theoretical model to unravel the effect of polarity in packing of noncentrosymmetric domains on SFG spectral features has remained unresolved. In this work, we show how the phase synchronization principle of the SFG process is used to predict the relative intensities of vibrational modes with different polar angles from the noncentrosymmetric domain axis. Applying this model calculation for the first time and employing SFG microscopy, we found that cellulose microfibrils in certain xylem cell walls are deposited unidirectionally (or biased in one direction) instead of the bidirectional polarity which was believed to be dominant in plant cell walls from volume-averaged characterizations of macroscopic samples. Here, with this advancement in SFG analysis, one can now determine the relative polarity of noncentrosymmetric domains such as crystalline biopolymers interspersed in amorphous polymer matrices, which will open opportunities to study new questions that have not been conceived in the past.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spectroscopic Studies of Charge-Transfer Character and Photoresponses of F 4 TCNQ-Based Donor–Acceptor Complexes

F 4 TCNQ (2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane) is used widely as a hole-doping agent in photoresponsive organic semiconducting materials, yet relatively little is known about the photoresponses of the F 4 TCNQ· - anion generated via doping. Furthermore, there is still relatively little systematic exploration of how the properties of the local material or chemical environment impacts the driving force for generating these charge-transfer complexes. Here we present spectroscopic and photophysical studies of F 4 TCNQ in charge-transfer complexes (CTCs) with the electron donor N,N’-Diphenyl-N-N’-di-p-tolylbenzene-1,4-diamine (MPDA) both in dichloroethane solution and polystyrene matrices. Integer charge transfer (ICT) between donor and acceptor occurs readily in dichloroethane solvent to form F 4 TCNQ· - :MPDA + CTCs, due to a ~150 mV difference in MPDA + /MPDA and F 4 TCNQ/F 4 TCNQ· - reduction potentials. Ultrafast spectroscopic studies of the CTC as well as electrochemically generated F 4 TCNQ· - and MDPA + reveal that the photoresponses of these CTCs are dominated by that of the dopant anion, including rapid deactivation (800 fs) after excitation to the anion D 1 excited state, followed by slower (~10 ps) vibrational cooling in the anion D 0 state. Excitation to the higher-lying D 2 state results in a rapid relaxation to the D 1 state, in contrast to direct D 2 -D 0 relaxation previously observed for F 4 TCNQ· - in gas phase. CTCs embedded in polystyrene (PS) matrices are observed to lose their integer charge-transfer character upon evaporation of solvent, as evidenced by changes to electronic and vibrational absorption features associated with F 4 TCNQ· - . This change is attributed to the loss of solvent stabilization of the ion pair formed through the charge-transfer reaction. Ultrafast spectral measurements reveal that the photoresponses of the partial charge-transfer (PCT) species embedded in PS are still highly similar to those of the ICT species and unlike that of neutral F 4 TCNQ, implying the electronic properties of the PCT state are likewise dominated by properties of the reduced acceptor molecule. We conclude that excitation of ICT or PCT states introduces optical losses for photoresponses of doped organic semiconductor materials due to the large anion absorption cross section and its rapid, dissipative deactivation dynamics.

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

Reverse Polarizability of Rare Earth Ions (La 3+ , Gd 3+ , Lu 3+ , Y 3+ ) in Tellurite Glasses and Glass Ceramics for Optical Limiting

All-optical modulation using inherent third-order optical nonlinearity of a medium has garnered considerable interest in photonics and optoelectronics. Herein, nonlinear optical (NLO) properties of tellurite glasses and glass ceramics (GCs) containing four different rare earths (RE = La, Gd, Lu, and Y) have been deliberated in near-infrared regions under an ultrafast regime. The La-based glass exhibits ~10 times higher nonlinear refraction (n 2 ) and absorption (α 2 ) than reported NLO materials. The NLO susceptibility [χ (3) ] trend in the studied glasses is La > Gd > Lu > Y, matching with RE 3+ polarizability. Furthermore, Ln 2 Te 6 O 15 nanocrystallite-embedded transparent GCs exhibit a larger NLO coefficient due to the enhanced local field from oxygen vacancies in crystallites. Interestingly, the trend of χ (3) in GCs follows the sequence of Y > Lu > Gd > La, precisely opposite to the glasses. This observation challenges the general polarizability approach of RE 3+ ions, emphasizing that quadratic hyperpolarizability of RE 3+ is pivotal for NLO properties of GCs. Among the studied matrices, Y-containing GCs showed the lowest optical limiting (OL) threshold (5.4 mJ/cm 2 at 800 nm), much lower than those of the reported NLO materials, suggesting its potential as a femtosecond NIR-laser safety material. A combination of large α 2 and n 2 from the studied matrices indicates their advantage for harmonic generation, potentially aiding in the design of ultrafast signal processing devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polymer-Grafted Cellulose Nanofibrils with Enhanced Interfacial Compatibility for Stronger Poly(lactic acid) Composites

Cellulose nanofibrils (CNFs) are a promising reinforcement for biodegradable composite matrices such as poly(lactic acid) (PLA), but they require commercially scalable drying methods that preserve their fibrillar morphology along with improved interfacial interactions with polymer matrices. In this work, a water-based grafting-through polymerization scheme to modify CNFs improved spray drying behavior and reinforcement capacity in PLA composites. All polymer modifications yielded CNFs with a more fibrillar morphology after spray drying, increasing specific surface area by up to 490% compared to unmodified CNFs, values similar to conventional freeze drying. Polymer-grafted CNFs in PLA composites improved the tensile strength by 16% at 20 wt % loading and stiffness by 22% at a 10 wt % loading with two different graft-polymer chemistries compared to unmodified CNF composites. Surface energy heterogeneity measurements of the reinforcements and PLA matrix were employed to understand the improvements in composite properties. Polymer modifications lowered the total surface energy of the CNFs, and calculated ratios of work of adhesion to work of cohesion suggested improved interfacial compatibility for four of the modified CNFs with PLA. Rheological oscillatory shear studies of the composites correlated solid-like melt behavior, as demonstrated by storage moduli dominance, with higher tensile strength. Thermal analysis of the composites revealed that excessive plasticization by the poly(oligoethylene glycol methyl ether methacrylate)-grafted sample potentially offset mechanical property improvements imparted by the more fibrillar morphology. This work provides an opportunity for large-scale manufacturing of CNF/PLA composites via an entirely aqueous modification scheme and industrially relevant spray drying process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Relationship between microporous structure and light gas transport through glassy polymeric membranes revealed by molecular simulations

Microporous glassy polymers are attractive materials for gas separation membranes, due to their high permeability and tailorable selectivity, provided physical aging can be delayed. The archetypal microporous glassy polymer, PTMSP, can be blended with a hyper-crosslinked isatin–triptycene porous polymer network (PPN) to delay physical aging. However, while PPN is effective at reducing physical aging, it also affects the permeability of light gases through PTMSP. Molecular dynamics simulations were used here to shed fundamental light on the mechanisms responsible for these effects. Atomistic models are developed that satisfactorily reproduce experimental observations such as matrix density and cavity size distributions for neat PTMSP as well as for PTMSP–PPN blends. Analysis of the simulation results suggests that physical aging is delayed because the PPN inclusions slow down PTMSP relaxation while reducing the connectivity between free volume pockets. To understand how PPN inclusions affect light gas permeability, the atomistic models developed are used to probe CO 2 and CH 4 diffusion and sorption. These simulations are conducted for PTMSP matrices exhibiting varying density, towards reproducing experimental permeability data. Interrogating the simulation trajectories, it is discovered that while CH 4 travels preferentially through the free volume cavities, CO 2 preferentially interacts with the available surfaces, especially in the presence of PPN. These differences help interpret experimental observations. The transport of C 2 H 6 and H 2 S through PTMSP matrices was also investigated. These gases, and in particular H 2 S, were found to absorb within PTMSP, yielding very low diffusion coefficients. The ability to predict differences in diffusion pathways for various natural gas components unveils the possibility of engineering membranes to control permeability and selectivity towards large-scale applications.

Bao Le, Tran Thi [Univ. of Oklahoma, Norman, OK (U↗
DOE OSTI2021

Imaging pyrometry for most color cameras using a triple pass filter

A simple combination of the Planck blackbody emission law, optical filters, and digital image processing is demonstrated to enable most commercial color cameras (still and video) to be used as an imaging pyrometer for flames and explosions. The hardware and data processing described take advantage of the color filter array (CFA) that is deposited on the surface of the light sensor array present in most digital color cameras. In this work, a triple-pass optical filter incorporated into the camera lens allows light in three 10-nm wide bandpass regions to reach the CFA/light sensor array. These bandpass regions are centered over the maxima in the blue, green, and red transmission regions of the CFA, minimizing the spectral overlap of these regions normally present. A computer algorithm is used to retrieve the blue, green, and red image matrices from camera memory and correct for remaining spectral overlap. A second algorithm calibrates the corrected intensities to a gray body emitter of known temperature, producing a color intensity correction factor for the camera/filter system. The Wien approximation to the Planck blackbody emission law is used to construct temperature images from the three color (blue, green, red) matrices. A short pass filter set eliminates light of wavelengths longer than 750 nm, providing reasonable accuracy (±10%) for temperatures between 1200 and 6000 K. The effectiveness of this system is demonstrated by measuring the temperature of several systems for which the temperature is known.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

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

Matrix-formation dynamics dictate methyl nitrite conformer abundance

Methyl nitrite has two stable conformational isomers resulting from rotation about the primary C–O–N–O dihedral angle: cis-CH 3 ONO and trans-CH 3 ONO, with cis being more stable by ∼5 kJ/mol. Here, the barrier to rotational interconversion (∼45 kJ/mol) is too large for isomerization to occur under ambient conditions. This paper presents evidence of a change in conformer abundance when dilute CH 3 ONO is deposited onto a cold substrate; the relative population of the freshly deposited cis conformer is seen to increase compared to its gas-phase abundance, measured by in situ infrared spectroscopy. We observe abundance changes depending on the identity of the bath gas (N 2 , Ar, and Xe) and deposition angle. The observations indicate that the surface properties of the growing matrix influence conformer abundance—contrary to the widely held assumption that conformer abundance in matrices reflects gas-phase abundance. We posit that differences in the angle-dependent host-gas deposition dynamics affect the growing surfaces, causing changes in conformer abundances. Quantum chemistry calculations of the binding energies between CH 3 ONO and a single bath-gas component reveal that significant energetic stabilization is not observed in 1:1 complexes of N 2 :CH 3 ONO, Ar:CH 3 ONO, or Xe:CH 3 ONO. From our results, we conclude that the growing surface plays a significant role in trapping cis-CH 3 ONO more effectively than trans-CH 3 ONO, likely because cis-CH 3 ONO is more compact. Taken together, the observations highlight the necessity for careful characterization of conformers in matrix-isolated systems, emphasizing a need for further study into the deposition dynamics and surface structure of chemically inert matrices.

Hockey, Emily K.↗
DOE OSTI2020

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Particle hit clustering and identification using point set transformers in liquid argon time projection chambers

Liquid argon time projection chambers are often used in neutrino physics and dark-matter searches because of their high spatial resolution. The images generated by these detectors are extremely sparse, as the energy values detected by most of the detector are equal to 0, meaning that despite their high resolution, most of the detector is unused in a particular interaction. Instead of representing all of the empty detections, the interaction is usually stored as a sparse matrix, a list of detection locations paired with their energy values. Traditional machine learning methods that have been applied to particle reconstruction such as convolutional neural networks (CNNs), however, cannot operate over data stored in this way and therefore must have the matrix fully instantiated as a dense matrix. Operating on dense matrices requires a lot of memory and computation time, in contrast to directly operating on the sparse matrix. We propose a machine learning model using a point set neural network that operates over a sparse matrix, greatly improving both processing speed and accuracy over methods that instantiate the dense matrix, as well as over other methods that operate over sparse matrices. Compared to competing state-of-the-art methods, our method improves classification performance by 14%, segmentation performance by more than 22%, while taking 80% less time and using 66% less memory. Compared to state-of-the-art CNN methods, our method improves classification performance by more than 86%, segmentation performance by more than 71%, while reducing runtime by 91% and reducing memory usage by 61%.

calibration and fitting methods↗
DOE OSTI2020

A neural network for determination of latent dimensionality in Nonnegative Matrix Factorization

Non-negative Matrix Factorization (NMF) has proven to be a powerful unsupervised learning method for uncovering hidden features in complex and noisy datasets with applications in data mining, text recognition, dimension reduction, face recognition, anomaly detection, blind source separation, and many other fields. An important input for NMF is the latent dimensionality of the data, that is, the number of hidden features, K, present in the explored dataset. Unfortunately, and this quantity is rarely known a priori. The existing methods for determining latent dimensionality, such as Automatic Relevance Determination (ARD), are mostly heuristic and utilize different characteristics to estimate the number of hidden features. However, all of them require human presence to make a final determination of K. Here we utilize a supervised machine learning approach in combination with a recent method for model determination, called NMFk, to determine the number of hidden features automatically. NMFk performs a set of NMF simulations on an ensemble of matrices, obtained by bootstrapping the initial dataset, and estimates which K produces stable groups of latent features that reconstruct the initial dataset well. We then train a Multi-Layter Perceptron (MLP) classifier network to determine the correct number of latent features utilizing the statistics and characteristics of the NMF solution, obtained from NMFk. In order to train the MLP classifier, a training set of 58,660 matrices with predetermined latent features were factorized with NMFk. The MLP classifier in conjunction with NMFk maintains a greater than 95% success rate when applied to a held out test set. Additionally, when applied to two well-known benchmark datasets, the swimmer and MIT face data, NMFk/MLP correctly recovers the established number of hidden features. Finally, we compare the accuracy of our method to the ARD, AIC and Stability-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗
DOE OSTI2026

Linear complexity

We present factorization and solution phases for a new linear complexity direct solver designed for concurrent batch operations on fine-grained parallel architectures, for matrices amenable to hierarchical representation. We focus on the strong-admissibility-based $\mathscr{H}^{2}$ format, where strong recursive skeletonization factorization compresses remote interactions. We build upon previous implementations of $\mathscr{H}^{2}$ matrix construction for efficient factorization and solution algorithm design, which are illustrated graphically in stepwise detail. The algorithms are ‘blackbox’ in the sense that the only inputs are the matrix and right-hand side, without analytical or geometrical information about the origin of the system. We demonstrate linear complexity scaling in both time and memory on four representative families of dense matrices up to one million in size. Parallel scaling up to 16 threads is enabled by a multi-level matrix graph coloring and avoidance of dynamic memory allocations thanks to prefix-sum memory management. An experimental backward error analysis is included. We break down the timings of different phases, identify phases that are memory-bandwidth limited, and discuss alternatives for phases that may be sensitive to the trend to employ lower precisions for performance.

Boukaram, Wajih↗
DOE OSTI2025

Climate history modulates stress responses of common soil bacteria under experimental drought

Soil drying challenges microbial viability and survival, with bacteria employing various mechanisms to respond to shifts in osmolarity, including dormancy or metabolic upregulation of osmoprotectants. However, the extent to which these responses are shaped by an organism's phylogeny, or the climate history of a given environment is poorly understood. This study examines the responses of phylogenetically similar bacteria from semi-arid and humid tropical forest soils to osmotic and matric stress using synchrotron radiation-based Fourier Transform Infrared spectromicroscopy. This non-destructive approach depicts the biochemical phenotype for whole cells under control and stress conditions. We observed that, under osmotic stress, bacteria upregulated cell-signaling pathways, rapidly turned over lipid-storage compounds, and increased osmolyte production. In contrast, matric stress induced a more muted response, typically elevating the production of carbohydrate stress compounds, such as glycine betaine and trehalose. Whereas phylogenetically similar bacteria showed comparable biochemistry under control conditions, climate history played an important role in regulating responses to stress, whereby a stronger metabolic response was observed from semi-arid relative to tropical forest isolates. We conclude that bacterial stress response to drought can be more diverse than previously observed and regulated by both phylogeny and climate history.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Locally adaptive activation functions with slope recovery for deep and physics-informed neural networks

Here we propose two approaches of locally adaptive activation functions namely, layer-wise and neuron-wise locally adaptive activation functions, which improve the performance of deep and physics-informed neural networks. The local adaptation of activation function is achieved by introducing a scalable parameter in each layer (layer-wise) and for every neuron (neuron-wise) separately, and then optimizing it using a variant of stochastic gradient descent algorithm. In order to further increase the training speed, an activation slope-based slope recovery term is added in the loss function, which further accelerates convergence, thereby reducing the training cost. On the theoretical side, we prove that in the proposed method, the gradient descent algorithms are not attracted to sub-optimal critical points or local minima under practical conditions on the initialization and learning rate, and that the gradient dynamics of the proposed method is not achievable by base methods with any (adaptive) learning rates. We further show that the adaptive activation methods accelerate the convergence by implicitly multiplying conditioning matrices to the gradient of the base method without any explicit computation of the conditioning matrix and the matrix–vector product. The different adaptive activation functions are shown to induce different implicit conditioning matrices. Furthermore, the proposed methods with the slope recovery are shown to accelerate the training process.

97 MATHEMATICS AND COMPUTING↗