Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Pan-spectral observing system simulation experiments of shortwave reflectance and long-wave radiance for climate model evaluation

Abstract. Top-of-atmosphere (TOA) spectrally resolved shortwave reflectances and long-wave radiances describe the response of the Earth's surface and atmosphere to feedback processes and human-induced forcings. In order to evaluate proposed long-duration spectral measurements, we have projected 21st Century changes from the Community Climate System Model (CCSM3.0) conducted for the Intergovernmental Panel on Climate Change (IPCC) A2 Emissions Scenario onto shortwave reflectance spectra from 300 to 2500 nm and long-wave radiance spectra from 2000 to 200 cm –1 at 8 nm and 1 cm –1 resolution, respectively. The radiative transfer calculations have been rigorously validated against published standards and produce complementary signals describing the climate system forcings and feedbacks. Additional demonstration experiments were performed with the Model for Interdisciplinary Research on Climate (MIROC5) and Hadley Centre Global Environment Model version 2 Earth System (HadGEM2-ES) models for the Representative Concentration Pathway 8.5 (RCP8.5) scenario. The calculations contain readily distinguishable signatures of low clouds, snow/ice, aerosols, temperature gradients, and water vapour distributions. The goal of this effort is to understand both how climate change alters reflected solar and emitted infrared spectra of the Earth and determine whether spectral measurements enhance our detection and attribution of climate change. This effort also presents a path forward to understand the characteristics of hyperspectral observational records needed to confront models and inline instrument simulation. Such simulation will enable a diverse set of comparisons between model results from coupled model intercomparisons and existing and proposed satellite instrument measurement systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Impacts of future climate change on urban flood volumes in Hohhot in northern China: benefits of climate change mitigation and adaptations

As China becomes increasingly urbanised, flooding has become a regular occurrence in its major cities. Assessing the effects of future climate change on urban flood volumes is crucial to informing better management of such disasters given the severity of the devastating impacts of flooding (e.g. the 2016 flooding events across China). Although recent studies have investigated the impacts of future climate change on urban flooding, the effects of both climate change mitigation and adaptation have rarely been accounted for together in a consistent framework. In this study, we assess the benefits of mitigating climate change by reducing greenhouse gas (GHG) emissions and locally adapting to climate change by modifying drainage systems to reduce urban flooding under various climate change scenarios through a case study conducted in northern China. The urban drainage model – Storm Water Management Model – was used to simulate urban flood volumes using current and two adapted drainage systems (i.e. pipe enlargement and low-impact development, LID), driven by bias-corrected meteorological forcing from five general circulation models in the Coupled Model Intercomparison Project Phase 5 archive. Results indicate that urban flood volume is projected to increase by 52 % over 2020–2040 compared to the volume in 1971–2000 under the business-as-usual scenario (i.e. Representative Concentration Pathway (RCP) 8.5). The magnitudes of urban flood volumes are found to increase nonlinearly with changes in precipitation intensity. On average, the projected flood volume under RCP 2.6 is 13 % less than that under RCP 8.5, demonstrating the benefits of global-scale climate change mitigation efforts in reducing local urban flood volumes. Comparison of reduced flood volumes between climate change mitigation and local adaptation (by improving drainage systems) scenarios suggests that local adaptation is more effective than climate change mitigation in reducing future flood volumes. This has broad implications for the research community relative to drainage system design and modelling in a changing environment. Furthermore, this study highlights the importance of accounting for local adaptation when coping with future urban floods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

DECOVALEX-2015: an international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

Here, an important part of the performance and safety assessment of disposal systems for radioactive waste and spent nuclear fuel in deep geological formations is to evaluate the impact of the coupled effects of mechanical deformation, fluid and gas flow through the repository, and thermal loading from the decaying waste. It was recognized early on in safety and performance assessments conducted for geologic disposal sites that to be able to conduct such an evaluation, there was a need to rigorously enhance the theoretical background and to develop models capable of simulating coupled thermo-hydro-mechanical (THM) processes. More recently, chemical processes have also been added to enable the study of fully coupled THMC processes in geosystems.

58 GEOSCIENCES↗
DOE OSTI2025

AMOC Variability in Climate Models and Its Dependence on the Mean State

Abstract Understanding internal variability of the climate system is critical when isolating internal and anthropogenically forced signals. Here, we investigate the modes of Atlantic Meridional Overturning Circulation ( AMOC ) variability using perturbation experiments with the Institut Pierre‐Simon Laplace's (IPSL) coupled model and compare them to Coupled Model Intercomparison Project Phase 6 (CMIP6) pre‐industrial control simulations. We identify two characteristic modes of variability—decadal‐to‐multidecadal ( DMD var ) and centennial ( CEN var ). The former is driven largely by temperature anomalies in the subpolar North Atlantic, while the latter is driven by salinity in the western subpolar North Atlantic. The amplitude of each mode scales linearly with the mean AMOC strength in the IPSL experiments. The DMD var amplitude correlates well with the AMOC mean strength across CMIP6 models, while the CEN var mode does not. These findings suggest that the strength of DMD var depends robustly on the North Atlantic mean state, while the CEN var mode may be model‐dependent.

Ferster, Brady S. [LOCEAN‐IPSL Sorbonne Université↗
DOE OSTI2021

Subtropical Eastern North Pacific SST bias in earth system models

Karthik Balaguru (PNNL), Luke Van Roekel (LLNL), L. Ruby Leung (PNNL) and Milena Veneziani (LLNL) Abstract: This study systematically evaluates the warm SST bias in the Subtropical Eastern North Pacific (SENP), a problem plaguing most CMIP6 models, using the Energy Exascale Earth System Model version 1 (E3SM). Compared to the model at its standard resolution (1o atmosphere, 30-60 km ocean), the SST bias improves in several regions at high resolution (0.25o atmosphere, 18-6 km ocean) except in subtropical eastern boundary regions, such as the SENP. A mixed layer heat budget analysis reveals that errors in the surface radiative fluxes and horizontal heat advection are primarily responsible for the bias. Analysis of HighResMIP models indicates that the lack of improvement in the SENP SST bias with resolution is not native to E3SM alone. A larger CMIP6 inter-model correlation of SST bias with biases in surface winds than radiation highlights the important role of large-scale circulation biases in the SENP SST bias in coupled models.

Balaguru, Karthik↗
DOE OSTI2021

Assessment of the burning-plasma operational space in ITER by using a control-oriented core-SOL-divertor model

In future tokamaks, the control of burning plasmas will require careful regulation of the plasma density and temperature. Along with the design of effective burn-control systems, understanding how the fusion power varies in the density-temperature space is vital for the operation of fusion power plants. Here in this work, the steady-state operational space of ITER is studied using a control-oriented core-plasma model coupled to a two-point model of the scrape-off-layer (SOL) and divertor regions. The two models are coupled through the exchange of input-output parameters. The deuterium and tritium recycling from the wall are output parameters of the SOL-divertor model that are used as input parameters in the core-plasma density balance. Furthermore, the separatrix temperature, which is an output parameter of the SOL-divertor model, is incorporated into the radial core-plasma temperature profiles. Therefore, the temperature-dependent power balance of the plasma core is intimately linked to the SOL-divertor model. Both the power entering the SOL from the core, as determined by the core-plasma power balance, and the separatrix density, as dictated by the core-plasma density balance, are input parameters to the SOL-divertor model. They are control knobs in the SOL-divertor model that can be regulated using the core-plasma actuators: auxiliary power and pellet injection. There are various operational limitations, such as the saturation of the aforementioned actuators, that will prevent ITER from accessing certain high-fusion plasma regimes. The achievable tritium concentration in the fueling lines and the maximum sustainable heat load on the divertor will impose further restrictions. By accounting for these limitations, the ITER operational space is computed based on the coupled core-SOL-divertor model and visualized using Plasma Operation Contour (POPCON) plots that map performance metrics, such as the fusion to auxiliary power ratio, over the density-temperature space. Comparisons are drawn between plasmas with different recycling, confinement, and SOL-divertor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

MIRAGE: A Framework for Data-Driven Collaborative High-Resolution Simulation

Information about how human populations shift in response to various stimuli is limited because no single model is capable of addressing these stimuli simultaneously, and integration of the best existing models has been challenging because of the vast disparity among constituent model purposes, architectures, scales, and execution environments. To demonstrate a potential model coupling for approaching this problem, three major model components are integrated into a fully coupled system that executes a worldwide infection-infected routine where a human population requires a food source for sustenance and an infected population can spread an infection when it is in contact with the remaining healthy population. To enable high-resolution data-driven model federation and an ability to capture dynamics and behaviors of billions of humans, a high-performance computing agent-based framework has been created and is demonstrated in this chapter.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Superconducting and structural properties of {delta}-MoC{sub 0.681} cubic molybdenum carbide phase

The superconducting and lattice properties of {delta}-MoC{sub 0.681} were studied by electromagnetic measurements, synchrotron X-ray diffraction, neutron diffraction, and electron diffraction. The superconducting properties (T{sub c}=12 K) of {delta}-MoC{sub 0.681} were well characterized by a weak coupling model. The carbon vacancies present in the host cubic structure were found to be robust, although the material was synthesized from stoichiometric carbon and Mo powder under a high-pressure of 6 GPa. A thermodynamically-stable structure with ordered vacancies did not account for the robust features of {delta}-MoC{sub 0.681} since the vacancies are unlikely to be ordered in long range in the host structure. A model based on inherent phonon instability theoretically predicted for a stoichiometric MoC phase might be responsible for the robust features of {delta}-MoC{sub 0.681}. - Graphical Abstract: The cubic molybdenum carbide shows an excellent superconductivity with robust carbon vacancies. Inherent phonon instability theoretically predicted for a stoichiometric MoC phase might be responsible for the vacancies rather than a thermodynamically-stable structure with vacancies ordering. Highlights: Black-Right-Pointing-Pointer The 12 K superconductivity is well characterized by a weakly coupling model. Black-Right-Pointing-Pointer Carbon vacancies are robust and disordered in the cubic host structure. Black-Right-Pointing-Pointer Inherent phonon instability might be responsible for the robust carbon vacancies in {delta}-MoC{sub 0.681}.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Extreme Deformations in Lithium Ion Batteries

A simultaneously coupled modeling approach to study the electrochemical and thermal behavior of lithium-ion batteries under large mechanical deformation has been developed. The thermo-electrochemical pseudo-2D (P2D) battery model is coupled with a mechanical material model. Mechanical, thermal, and electrochemical models are implemented as user-defined sub-routines in the commercial multi-physics code LS-DYNA. The mechanical strain experienced by anode, cathode and separator results in thickness and porosity changes in each layer which in turn influences electrochemical behavior. The evolution of concentration profiles and cell potential are studied under different mechanical loading conditions. Internal short-circuits caused by mechanical deformation and corresponding physical behaviors are also elucidated. We discuss the competing effects of improved transport at higher temperature due to the internal short-circuit versus a drop in the effective ionic conductivity and electrolyte diffusivity due to mechanical deformation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Quantifying impacts of heat waves on power grid operation

Climate change is projected to cause an increase in the severity and frequency of extreme weather events such as heat waves and droughts. Such changes present planning and operating challenges and risks to many economic sectors. In the electricity sector, statistics of extreme events in the past have been used to help plan for future peak loads, determine associated infrastructure requirements, and evaluate operational risks, but industry-standard planning tools have yet to be coupled with or informed by temperature models to explore the impacts of the "new normal" on planning studies. For example, high ambient temperatures during heat waves reduce the output capacity and efficiency of gas fired combustion turbines just when they are needed most to meet peak demands. This paper describes the development and application of a production cost and unit commitment model coupled to high resolution, hourly temperature data and a temperature dependent load model. The coupled system has the ability to represent the impacts of hourly temperatures on load conditions and available capacity and efficiency of combustion turbines, and therefore capture the potential impacts on system reliability and production cost. Ongoing work expands this capability to address the impacts of water availability and temperature on power grid operation.

Climate change↗
DOE OSTI2024

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

CAMPUS: Fully Coupled Multiphysics Modeling and Simulation of LWR and FBR Nuclear Fuel Behavior

Nuclear fuels operate in an environment that involves complex multiphysics phenomena, and the multiphysics behavior is often tightly coupled. The prediction of fuel behavior under irradiation involves very complex material behavior that is not seen in any other engineering discipline. Understanding and predicting the evolving material properties and phenomenal behaviors of nuclear fuels are major challenges for fuel design and behavior prediction. The multiphysics behavior is often a tightly coupled thermal, mechanical and chemical phenomenon associated with a large number of control parameters and uncertainties, which requires an advanced fuel performance modeling and simulation capability to improve the existing reactors performance, to address environmental and reactor safety concerns. Nuclear fuel performance modeling and simulation capability is necessary as a result of the high cost and difficulty in experimental measurements. So, a fuel performance code called CityU Advanced Multiphysics Nuclear Fuels Performance with User-defined Simulations (CAMPUS) was developed to predict both light water reactor and fast reactor fuel behavior. An effective approach was developed using self-defined multiphysics models based on the framework of COMSOL Multiphysics to manage the nonlinearities associated with light water reactor and fast reactor fuels performance analysis. The modelling and simulations approach is of particular importance since it demonstrates highly interrelated and nonlinear fuel performance models can be developed and implemented into a multiphysics platform, and thus more efforts can be put into physics model development and improvement. The modelling results can be used to identify the important phenomena and parameters, to direct and design experiments and to assist in post analysis of expensive nuclear experiments for light water reactor and fast reactor nuclear fuel performance analysis. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2015

Rate turnover in mechano-catalytic coupling: A model and its microscopic origin

A novel aspect in the area of mechano-chemistry concerns the effect of external forces on enzyme activity, i.e., the existence of mechano-catalytic coupling. Recent experiments on enzyme-catalyzed disulphide bond reduction in proteins under the effect of a force applied on the termini of the protein substrate reveal an unexpected biphasic force dependence for the bond cleavage rate. Here, using atomistic molecular dynamics simulations combined with Smoluchowski theory, we propose a model for this behavior. For a broad range of forces and systems, the model reproduces the experimentally observed rates by solving a reaction-diffusion equation for a “protein coordinate” diffusing in a force-dependent effective potential. The atomistic simulations are used to compute, from first principles, the parameters of the model via a quasiharmonic analysis. Additionally, the simulations are also used to provide details about the microscopic degrees of freedom that are important for the underlying mechano-catalysis.

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

Twistor Actions for Integrable Systems

Many integrable systems can be reformulated as holomorphic vector bundles on twistor space. This is a powerful organizing principle in the theory of integrable systems. One shortcoming is that it is formulated at the level of the equations of motion. From this perspective, it is mysterious that integrable systems have Lagrangians. In this paper, we study a Chern-Simons action on twistor space and use it to derive the Lagrangians of some integrable sigma models. Our focus is on examples that come from dimensionally reduced gravity and supergravity. The dimensional reduction of general relativity to two spacetime dimensions is an integrable coset sigma model coupled to a dilaton and 2d gravity. The dimensional reduction of supergravity to two spacetime dimensions is an integrable coset sigma model coupled to matter fermions, a dilaton, and 2d supergravity. We derive Lax operators and Lagrangians for these 2d integrable systems using the Chern-Simons theory on twistor space. In the supergravity example, we use an extended setup in which twistor Chern-Simons theory is coupled to a pair of matter fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗