Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic coupling and spin-wave dispersions in a magnetic hybrid system made of an artificial spin-ice structure and an extended NiFe underlayer

We present a combined experimental and numerical study of the spin-wave dispersion in a NiFe artificial spin-ice (ASI) system consisting of an array of stadium-shaped nanoislands deposited on the top of a continuous NiFe film with non-magnetic spacer layers of varying thickness. The spin-wave dispersion, measured by wavevector resolved Brillouin light scattering spectroscopy in the Damon–Eshbach configuration, consists of a rich number of modes, with either stationary or propagating character. We find that the lowest frequency mode displays a bandwidth of ~0.5 GHz, which is independent of the presence of the film underneath. On the contrary, the Brillouin light scattering intensity of some of the detected modes strongly depends on the presence of the extended thin-film underlayer. Micromagnetic simulations unveil the details of the dynamic coupling between the ASI lattice and film underlayer. Interestingly, the ASI lattice facilitates dynamics of the film either specific wavelengths or intensity modulation peculiar to the modes of the ASI elements imprinted in the film. Our results demonstrate that propagating spin waves can be modulated at the nanometer length scale by harnessing the dynamic mode coupling in the vertical, i.e., the out-of-plane direction of suitably designed magnonic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ALP-Assisted Electroweak Phase Transition and Baryogenesis

Axion-like particles (ALPs) can be naturally lighter than the electroweak scale. We consider an ALP that couples to the Standard Model Higgs to achieve the strong first-order electroweak phase transition. We discuss the two-field dynamics of the phase transition and the associated computation in detail and identify the viable parameter space. The ALP mass can be from the MeV to GeV scale. Baryon asymmetry can be explained by local baryogenesis without violating the electron electric dipole moment bound. The viable parameter space can be probed through Higgs exotic decay, rare kaon decay, the electron electric dipole moment, and the effective number of neutrinos in the cosmic microwave background. The gravitational-wave signal is too weak to be detected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Prediction of Long-Term Geochemical Change in Bentonite Based on the Interpretative THMC Model of the FEBEX In Situ Test

Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.

Zheng, Liange↗
DOE OSTI2024

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal confinement and transport in spherical tokamaks: a review

Here, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger $E_r \times B$ shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high β T . Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower β T , both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher β T , microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of crystallization fouling on a horizontal-tube falling-film evaporator for thermal desalination

Scale formation within horizontal-tube falling-film evaporators is an important issue for thermal desalination, because of the deterioration of heat transfer performance and added maintenance costs. Predicting the fouling process is crucial for desalination plant design, operation, and maintenance. In this paper, a model is developed to predict the crystallization fouling process during seawater falling-film flow over a horizontal tube bundle. The scaling is estimated based on deposition theory, with couple heat transfer including the in-tube steam condensation, conduction through the tube wall and the scale layer, and falling-film evaporation. The spatial and temporal variations of temperature, heat transfer coefficient, and scale thickness are predicted, and the effects of variations in process parameters of steam, feed seawater, and tube properties on the scale thickness, evaporation rate, and heat transfer coefficient for falling-film evaporation are analyzed. Comparisons to existing experimental data show the scaling layer thickness falls in the general range of reported measurements, with a slight overestimation due to neglecting the delay of scaling onset under real conditions. The scale layer thickness increases dramatically as the in-tube steam pressure is increased or the seawater flow rate is decreased. In general, the scaling layer becomes thicker on the lower tubes in the tube bundle, due to the heat and mass transfer of the falling film. The effect of tube material on scaling appears mainly dependent on thermal conductivity; thus, using polymer rather than a stainless steel tube decreases scaling and evaporation rate 96% and 88%, respectively. Finally, this work has potential to guide thermal desalination plant design, material selection, and operating parameter optimization.

42 ENGINEERING↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2025

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Gravitational wave probes of massive gauge bosons at the cosmological collider

We extend the reach of the "cosmological collider" for massive gauge boson production during inflation from the CMB scales to the interferometer scales. Considering a Chern-Simons coupling between the gauge bosons and the pseudoscalar inflaton, one of the transverse gauge modes is efficiently produced and its inverse decay leaves an imprint in the primordial scalar and tensor perturbations. We study the correlation functions of these perturbations and derive the updated constraints on the parameter space from CMB observables. We then extrapolate the tensor power spectrum to smaller scales consistently taking into account the impact of the gauge field on inflationary dynamics. Our results show that the presence of massive gauge fields during inflation can be detected from characteristic gravitational wave signals encompassing the whole range of current and planned interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Giant phonon anomalies in the proximate Kitaev quantum spin liquid α-RuCl 3

The Kitaev quantum spin liquid epitomizes an entangled topological state, for which two flavors of fractionalized low-energy excitations are predicted: the itinerant Majorana fermion and the Z 2 gauge flux. Here, it was proposed recently that fingerprints of fractional excitations are encoded in the phonon spectra of Kitaev quantum spin liquids through a novel fractional-excitation-phonon coupling. Here, we detect anomalous phonon effects in α-RuCl 3 using inelastic X-ray scattering with meV resolution. At high temperature, we discover interlaced optical phonons intercepting a transverse acoustic phonon between 3 and 7 meV. Upon decreasing temperature, the optical phonons display a large intensity enhancement near the Kitaev energy, J K ~8 meV, that coincides with a giant acoustic phonon softening near the Z 2 gauge flux energy scale. These phonon anomalies signify the coupling of phonon and Kitaev magnetic excitations in α-RuCl 3 and demonstrates a proof-of-principle method to detect anomalous excitations in topological quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Future-proofing photovoltaics module reliability through a unifying predictive modeling framework

Solar energy, especially photovoltaics (PV), plays a significant role in the global energy transition required for decarbonization. Technological advancements in increasing efficiencies coupled with cost reductions through economies of scale made PV a cost-competitive energy source set to exceed the milestone of 1 TW globally deployed capacity in 2022. However, ongoing downward price pressure coupled with increasing service life expectations creates tremendous challenges for reliability engineers in ensuring the safe and reliable operation of PV modules and systems over the anticipated service life. Today's reliability research efforts aim for module lifetimes of up to 50 years. Still, our current tools fall short in accurately assessing degradation mechanisms and failure modes over such extended periods. We argue that the well-established PV reliability learning cycle needs to be accelerated to keep up with the rapid technological advancements and high expectations that are put on PV and its role in the global energy transition. In this article, we explore the evolution of the PV reliability learning cycle and highlight the significance that predictive modeling capabilities will have on future PV module reliability. We propose creating a unifying modeling framework - which once established will enable the holistic assessment of PV module reliability, accelerate the PV reliability learning cycle, and bring us closer to quantitative service life predictions of PV modules.

14 SOLAR ENERGY↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
DOE OSTI2021

Omics-to-Reactive-Transport (ORT): A workflow linking genome-scale metabolic models with reactive transport codes

Motivation: Nutrient and contaminant behavior in the subsurface are governed by multiple coupled hydrobiogeochemical processes which occur across different temporal and spatial scales. Accurate description of macroscopic system behavior requires accounting for the effects of microscopic and especially microbial processes. Microbial processes mediate precipitation and dissolution and change aqueous geochemistry, all of which impacts macroscopic system behavior. As `omics data describing microbial processes is increasingly affordable and available, novel methods for using this data quickly and effectively for improved ecosystem models are needed. Results: We propose a workflow (`Omics to Reactive Transport – ORT) for utilizing metagenomic and environmental data to describe the effect of microbiological processes in macroscopic reactive transport models. This workflow utilizes and couples two open-source software packages: KBase (a software platform for systems biology) and PFLOTRAN (a reactive transport modeling code). We describe the architecture of ORT and demonstrate an implementation using metagenomic and geochemical data from a river system. Our demonstration uses microbiological drivers of nitrification and denitrification to predict nitrogen cycling patterns which agree with those provided with generalized stoichiometries. While our example uses data from a single measurement, our workflow can be applied to spatiotemporal metagenomic datasets to allow for iterative coupling between KBASE and PFLOTRAN. Live, interactive models, which incorporate the results from this narrative into a PFLOTRAN simulation, are available (without login) at https://pflotranmodeling.paf.subsurfaceinsights.com/pflotran-simple-model/.

Rubinstein, Rebecca L↗
DOE OSTI2023

Gravity-improved metastability bounds for the Type-I seesaw mechanism

Right-handed neutrinos (RHN) destabilize the electroweak vacuum by increasing its decay rate. In the SM, the latter is dominated by physics at the RG scale at which λ reaches its minimum, ${\mu}_{\ast}^{\textrm{SM}}$ ~10 17 GeV. For large neutrino Yukawa coupling Y ν , RHNs can push μ* beyond the Planck scale, implying that gravitational effects need to be taken into account. In this work, we perform the first comprehensive study of electroweak vacuum metastability in the type-I seesaw mechanism including these effects. Our analysis covers both low- and high-scale seesaw models, with two as well as three RHNs and for multiple values of the Higgs’ non-minimal coupling to gravity. We find that gravitational effects can significantly stabilize the vacuum, leading to weaker metastability bounds. We show that metastability sets the strongest bounds for low-scale seesaws with M N > 1 TeV. For high-scale seesaws, we find upper bounds on the allowed masses for the RHNs, which are relevant for high-scale leptogenesis. We also point out that Tr(${Y}_{\nu}^{\dagger}$Y ν ), which is commonly used to express these metastability bounds, cannot be used for all of parameter space. Instead, we argue that bounds can always be expressed reliably through Tr(${Y}_{\nu}^{\dagger}$Y ν ${Y}_{\nu}^{\dagger}$Y ν ). Lastly, we use this insight to develop a new technique for an easier RG analysis applicable to scenarios with degenerate RHN masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗