SEARCH · Engineering Papers

Results for “weak scaling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Variability of Eastern North Atlantic Summertime Marine Boundary Layer Clouds and Aerosols Across Different Synoptic Regimes Identified With Multiple Conditions

Abstract This study estimates the meteorological covariations of aerosol and marine boundary layer (MBL) cloud properties in the eastern North Atlantic (ENA) region, characterized by diverse synoptic conditions. Using a deep‐learning‐based clustering model with mid‐level and surface daily meteorological data, we identify seven distinct synoptic regimes during the summer from 2016 to 2021. Our analysis, incorporating reanalysis data and satellite retrievals, shows that surface aerosols and MBL clouds exhibit clear regime‐dependent characteristics, whereas lower tropospheric aerosols do not. This discrepancy likely arises from synoptic regimes determined by daily large‐scale conditions, which may overlook air mass histories that predominantly dictate lower tropospheric aerosol conditions. Focusing on three regimes dominated by northerly winds, we analyze the Atmospheric Radiation Measurement Program (ARM) ENA observations on Graciosa Island in the Azores. In the subtropical anticyclone regime, fewer cumulus clouds and more single‐layer stratocumulus clouds with light drizzle are observed, along with the highest cloud droplet number concentration (Nd), surface cloud condensation nuclei (CCN) and surface aerosol levels. The post‐trough regime features more broken or multi‐layer stratocumulus clouds with slightly higher surface rain rate, and lower Nd and surface CCN levels. The weak trough regime is characterized by the deepest MBL clouds, primarily cumulus and broken stratocumulus clouds, with the strongest surface rain rate and the lowest Nd, surface CCN and surface aerosol levels, indicating strong wet scavenging. These findings highlight the importance of considering the covariation of cloud and aerosol properties driven by large‐scale regimes when assessing aerosol indirect effects using observations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Eulerian thermo-mechanical simulations of heterogeneous solid propellants using an approximate projection method

We present a mathematical model that describes thermo-mechanical deformations and thermal gradients on the unsteady burning of a heterogeneous solid propellant. A scaling study shows that the deformations in solid at combustion timescales can be treated as quasi-static. The resulting thermo-mechanical formulation is formulated on a Cartesian grid and makes use of a weak form of Chorin-type projection method to deal with large difference in shear modulus of constituent materials. Further, a one-dimensional verification study is carried out by comparing numerical simulations with those of an analytical model. In addition, convergence studies for a two-dimensional propellant sandwich configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a two-dimensional random propellant pack and the time-averaged burn rate and strains are reported.

42 ENGINEERING↗

DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Massively parallel and universal approximation of nonlinear functions using diffractive processors

Nonlinear computation is essential for a wide range of information processing tasks, yet implementing nonlinear functions using optical systems remains a challenge due to the weak and power-intensive nature of optical nonlinearities. Overcoming this limitation without relying on nonlinear optical materials could unlock unprecedented opportunities for ultrafast and parallel optical computing systems. Here, we demonstrate that large-scale nonlinear computation can be performed using linear optics through optimized diffractive processors composed of passive phase-only surfaces. In this framework, the input variables of nonlinear functions are encoded into the phase of an optical wavefront—e.g., via a spatial light modulator (SLM)—and transformed by an optimized diffractive structure with spatially varying point-spread functions to yield output intensities that approximate a large set of unique nonlinear functions–all in parallel. We provide proof establishing that this architecture serves as a universal function approximator for an arbitrary set of bandlimited nonlinear functions, also covering wavelength-multiplexed nonlinear functions as well as multi-variate and complex-valued functions that are all-optically cascadable. Our analysis also indicates the successful approximation of typical nonlinear activation functions commonly used in neural networks, including the sigmoid, tanh, ReLU (rectified linear unit), and softplus. We numerically demonstrate the parallel computation of one million distinct nonlinear functions, accurately executed at wavelength-scale spatial density at the output of a diffractive optical processor. Furthermore, we experimentally validated this framework using in situ optical learning and approximated 35 unique nonlinear functions in a single shot using a compact setup consisting of an SLM and an image sensor. These results establish diffractive optical processors as a scalable platform for massively parallel universal nonlinear function approximation, paving the way for new capabilities in analog optical computing based on linear materials.

Rahman, Md Sadman Sakib [University of California,↗

DOE OSTI2022

A Massive, Position-Resolving, High-Energy-Resolution Detector for Non-Accelerator Cosmic and Intensity Frontier Particle Physics

We proposed to take the first steps in the development of a detector that promises energy resolution of tens of eV FWHM combined with robust nuclear-recoil discrimination and high fidelity, mm-precision position reconstruction for applications in non-accelerator particle physics at the Cosmic and Intensity Frontiers. The detector would obtain these excellent resolutions by sensing athermal phonons produced by particle interactions in crystalline, dielectric targets using a sensitive, highly multiplexable superconducting phonon sensor, the kinetic inductance detector (KID). This detector would be applicable to: the search for low-mass particle dark matter candidates with masses below 5~GeV via direct detection of scattering of dark matter particles with terrestrial nuclei; detection of coherent elastic neutrino-nucleus scattering to test for new physics such as a non-standard value of the weak nuclear charge, non-standard neutrino interactions (perhaps driven by a neutrino magnetic moment), or the existence of sterile neutrinos; and, searches for neutrinoless double-beta decay. During the funding period, we demonstrated scaling up of the detector concept from a 22-mm by 22-mm by 1-mm, 1-g prototype with 0.9~keV FWHM energy resolution to a 75-mm diameter by 1-mm, 9-g prototype while improving the inferred energy resolution to 0.7~keV~FWHM. In the process, we solved many problems associated with scaling device fabrication to large wafers and vastly reduced the fraction of the surface occupied by inactive but energy-absorbing metal. We also demonstrated a new technique that substantially simplifies the process of characterizing a new detector. These results provide a good foundation for future work scaling up the design to 4-mm thickness substrates and improving the resolution to reach 0.035~keV~FWHM, yielding a detector with compelling potential for the above applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Nematic Fluctuations in the Non-Superconducting Iron Pnictide BaFe 1.9-x Ni 0.1 Cr x As 2

The main driven force of the electronic nematic phase in iron-based superconductors is still under debate. Here, we report a comprehensive study on the nematic fluctuations in a non-superconducting iron pnictide system BaFe 1.9-x Ni 0.1 Cr x As 2 by electronic transport, angle-resolved photoemission spectroscopy (ARPES), and inelastic neutron scattering (INS) measurements. Previous neutron diffraction and transport measurements suggested that the collinear antiferromagnetism persists to x = 0.8, with similar Néel temperature TN and structural transition temperature Ts around 32 K, but the charge carriers change from electron type to hole type around x = 0.5. In this study, we have found that the in-plane resistivity anisotropy also highly depends on the Cr dopings and the type of charge carriers. While ARPES measurements suggest possibly weak orbital anisotropy onset near Ts for both x = 0.05 and x = 0.5 compounds, INS experiments reveal clearly different onset temperatures of low-energy spin excitation anisotropy, which is likely related to the energy scale of spin nematicity. These results suggest that the interplay between the local spins on Fe atoms and the itinerant electrons on Fermi surfaces is crucial to the nematic fluctuations of iron pnictides, where the orbital degree of freedom may behave differently from the spin degree of freedom, and the transport properties are intimately related to the spin dynamics.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Environmental sensitivities of shallow-cumulus dilution – Part 1: Selected thermodynamic conditions

Cumulus entrainment, and its consequent dilution of buoyant cloud cores, strongly regulates the life cycle of shallow cumuli yet remains poorly understood. Herein, new insights into this problem are obtained through large-eddy simulations that systematically investigate the sensitivity of shallow-cumulus dilution to cloud-layer relative humidity (RH), cloud- and subcloud-layer depths, and continentality (i.e., the land–ocean contrast). The simulated cloud-core dilution is found to be strongly sensitive to continentality, with fractional dilution rates twice as large over the ocean as over land. Using a similarity theory based on the turbulent-kinetic-energy (TKE) budget, the reduced cloud-core dilution over land is attributed to larger cloud-base mass flux (mb), driven by stronger surface heating and subcloud turbulence. As mb increases, the fractional dilution rate must decrease to maintain energetic equilibrium. A positive sensitivity is also found to cloud-layer RH, with the core dilution increasing by 25 %–50 % for a 10 % enhancement in RH. This sensitivity is interpreted using the buoyancy-sorting hypothesis, in that mixtures of cloud and environmental air are more likely to become negatively buoyant and detrain (rather than diluting the cloud core) in drier cloud layers. By contrast, the sensitivities of (marine) shallow-cumulus dilution to cloud- and subcloud-layer depths are weak, with a 3 % decrease for a doubling for the former and a 4 % reduction in dilution for a 50 % deeper subcloud layer. These surprisingly weak sensitivities are readily explained by offsetting effects in the TKE similarity theory. Altogether, these experimental findings provide useful, though still incomplete, guidance for flow-dependent shallow-cumulus entrainment formulations in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Detection of the significant impact of source clustering on higher order statistics with DES Year 3 weak gravitational lensing data

We measure the impact of source galaxy clustering on higher order summary statistics of weak gravitational lensing data. By comparing simulated data with galaxies that either trace or do not trace the underlying density field, we show that this effect can exceed measurement uncertainties for common higher order statistics for certain analysis choices. We evaluate the impact on different weak lensing observables, finding that third moments and wavelet phase harmonics are more affected than peak count statistics. Using Dark Energy Survey (DES) Year 3 (Y3) data, we construct null tests for the source-clustering-free case, finding a p-value of p = 4 × 10 −3 (2.6σ) using third-order map moments and p = 3 × 10 −11 (6.5σ) using wavelet phase harmonics. The impact of source clustering on cosmological inference can be either included in the model or minimized through ad hoc procedures (e.g. scale cuts). We verify that the procedures adopted in existing DES Y3 cosmological analyses were sufficient to render this effect negligible. Failing to account for source clustering can significantly impact cosmological inference from higher order gravitational lensing statistics, e.g. higher order N-point functions, wavelet-moment observables, and deep learning or field-level summary statistics of weak lensing maps.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Redshift inference from the combination of galaxy colours and clustering in a hierarchical Bayesian model – Application to realistic N -body simulations

ABSTRACT Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broad-band imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. A hierarchical Bayesian model has recently been developed to estimate those from the robust combination of prior information, photometry of single galaxies, and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-z’s tightens the redshift posteriors and overcomes biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar{z} \approx 0.05,$ the final biases in the posterior are $\Delta \bar{z} \lesssim 0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Relaxation of the degenerate one-dimensional Fermi gas

In this work, we study how a system of one-dimensional spin-1/2 fermions at temperatures well below the Fermi energy approaches thermal equilibrium. The interactions between fermions are assumed to be weak and are accounted for within the perturbation theory. In the absence of an external magnetic field, spin degeneracy strongly affects relaxation of the Fermi gas. For sufficiently short-range interactions, the rate of relaxation scales linearly with temperature. Focusing on the case of the system near equilibrium, we linearize the collision integral and find the exact solution of the resulting relaxation problem. Furthermore, we discuss the application of our results to the evaluation of the transport coefficients of the one-dimensional Fermi gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Interaction-driven quantum anomalous Hall insulator in a Dirac semimetal

Here, the interaction-driven quantum anomalous Hall (QAH) insulator has been sought for a long time in a Dirac semimetal with linear band touching points at the Fermi level. By combining exact diagonalization, density matrix renormalization group, and analytical methods, we study a spinless fermion system on the checkerboard lattice with twofold rotational symmetry, which realizes two Dirac band touching points in the absence of interaction. At weak coupling, the Dirac semimetal is stable. At a finite density-density repulsive interaction, we analyze possible symmetry broken states, and find that an QAH state is stabilized when the interaction strength exceeds the energy scale controlling the separation between the Dirac points. Through numerical simulations, we verify the existence of the QAH phase with spontaneous time-reversal symmetry breaking and quantized Chern number C = 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

An acoustic investigation of the near-surface turbulence on Mars

The Perseverance rover is carrying out an original acoustic experiment on Mars: the SuperCam microphone records the spherical acoustic waves generated by laser sparks at distances from 2 m to more than 8 m. These N-shaped acoustic waves scatter from the multiple local heterogeneities of the turbulent atmosphere. Therefore, large and random fluctuations of sound travel time and intensity develop as the waves cross the medium. The variances of the travel times and the scintillation index (normalized variance of the sound intensity) are studied within the mathematical formalism of the propagation of spherical acoustic waves through thermal turbulence to infer statistical properties of the Mars atmospheric temperature fluctuation field. The comparison with the theory is made by simplifying assumptions that do not include wind fluctuations and diffraction effects. Two Earth years (about one Martian year) of observations acquired during the maximum convective period (10:00–14:00 Mars local time) show a good agreement between the dataset and the formalism: the travel time variance diverges from the linear Chernov solution exactly where the density of occurrence of the first caustic reaches its maximum. Moreover, on average, waves travel faster than the mean speed of sound due to a fast path effect, which is also observed on Earth. To account for the distribution of turbulent eddies, several power spectra are tested and the best match to observation is obtained with a generalized von Karman spectrum with a shallower slope than the Kolmogorov cascade, φ(k)∝(1 + k 2 L 2 ) -4/3 . It is associated with an outer scale of turbulence, L, of 11 cm at 2 m above the surface and a standard deviation of 6 K over 9 s for the temperature. These near-surface atmospheric properties are consistent with a weak to moderate wave scattering regime around noon with little saturation. Overall, this study presents an innovative and promising methodology to probe the near-surface atmospheric turbulence on Mars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Host Dark Matter Halos of SDSS Red and Blue Quasars: No Significant Difference in Large-scale Environment

The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models that attribute the color to the random orientation of the accretion disk along the line of sight, and evolutionary models that invoke connections between quasar systems and their environments. We test these schemas by probing the dark matter halo environments of optically selected quasars as a function of g – i optical color by measuring the two-point correlation functions of ~0.34 million eBOSS quasars as well as the gravitational deflection of cosmic microwave background photons around ~0.66 million XDQSO photometric quasar candidates. We do not detect a trend of halo bias with optical color through either analysis, finding that optically selected quasars at 0.8 < z < 2.2 occupy halos of characteristic mass M h ~ 3 × 10 12 h –1 M ⊙ regardless of their color. This result implies that a quasar's large-scale halo environment is not strongly connected to its observed optical color. We also confirm the findings of fundamental differences in the radio properties of red and blue quasars by stacking 1.4 GHz FIRST images at their positions, suggesting the observed differences cannot be attributed to orientation. Instead, the differences between red and blue quasars likely arise on nuclear-galactic scales, perhaps owing to reddening by a nuclear dusty wind. Finally, we show that optically selected quasars' halo environments are also independent of their r – W2 optical–infrared colors, while previous work has suggested that mid-infrared-selected obscured quasars occupy more massive halos. We discuss the implications of this result for models of quasar and galaxy coevolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Dark Energy Survey Year 3 Results: Multi-Probe Modeling Strategy and Validation

This paper details the modeling pipeline and validates the baseline analysis choices of the DES Year 3 joint analysis of galaxy clustering and weak lensing (a so-called "3$\times$2pt" analysis). These analysis choices include the specific combination of cosmological probes, priors on cosmological and systematics parameters, model parameterizations for systematic effects and related approximations, and angular scales where the model assumptions are validated. We run a large number of simulated likelihood analyses using synthetic data vectors to test the robustness of our baseline analysis. We demonstrate that the DES Year 3 modeling pipeline, including the calibrated scale cuts, is sufficiently accurate relative to the constraining power of the DES Year 3 analyses. Our systematics mitigation strategy accounts for astrophysical systematics, such as galaxy bias, intrinsic alignments, source and lens magnification, baryonic effects, and source clustering, as well as for uncertainties in modeling the matter power spectrum, reduced shear, and estimator effects. We further demonstrate excellent agreement between two independently-developed modeling pipelines, and thus rule out any residual uncertainties due to the numerical implementation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Direct Measurement of Galaxy Assembly Bias using DESI DR1 Data

We report the first direct measurement of galaxy assembly bias, a critical systematic in cosmology, from the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey. We introduce a novel, cosmology-independent method to measure the halo occupation distribution (HOD) by combining a state-of-the-art group catalog with weak gravitational lensing. For groups binned by total luminosity, we determine the galaxy occupation number $N_{\rm gal}$ from group-galaxy cross-correlations, while weak lensing constrains the average halo mass $M_h$. Applying this to a volume-limited sample at $z{\in}[0.05,0.2]$, we measure the dependence of HOD, $N_{\rm gal}(M_h)$, on large-scale overdensity $δ_{g}$. Focusing on the satellite galaxies, we find an assembly bias parameter of $Q_{\rm sat}{=}0.05{\pm}0.14$, a result consistent with zero and in tension with many empirical galaxy formation models. Our method provides a robust approach for characterizing galaxy assembly bias to achieve precision cosmology with DESI and future Stage-V surveys.

Shao, Zhiwei [Shanghai Jiaotong U.; Shanghai Jiao ↗

DOE OSTI2022

Probing the Neutrino-Mass Scale with the KATRIN Experiment

The absolute mass scale of neutrinos is an intriguing open question in contemporary physics. The as-yet-unknown mass of the lightest and, at the same time, most abundant massive elementary particle species bears fundamental relevance to theoretical particle physics, astrophysics, and cosmology. The most model-independent experimental approach consists of precision measurements of the kinematics of weak decays, notably tritium β decay. With the KATRIN experiment, this direct neutrino-mass measurement has entered the sub-eV domain, recently pushing the upper limit on the electron-based neutrino mass down to 0.8 eV (90% CL) on the basis of first-year data out of ongoing, multiyear operations. Here, we review the experimental apparatus of KATRIN, the progress of data taking, and initial results. While KATRIN is heading toward the target sensitivity of 0.2 eV, other scientific goals are pursued. We discuss the search for light sterile neutrinos and an outlook on future keV-scale sterile-neutrino searches as well as further physics opportunities beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗