skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Nuclear induction line shape: Non-Markovian diffusion with boundaries
The dynamics of viscoelastic fluids are governed by a memory function, essential yet challenging to compute, especially when diffusion faces boundary restrictions. We propose a computational method that captures memory effects by analyzing the time-correlation function of the pressure tensor, a viscosity indicator, through the Stokes–Einstein equation’s analytic continuation into the Laplace domain. We integrate this equation with molecular dynamics simulations to derive necessary parameters. Our approach computes nuclear magnetic resonance (NMR) line shapes using a generalized diffusion coefficient, accounting for temperature and confinement geometry. This method directly links the memory function with thermal transport parameters, facilitating accurate NMR signal computation for non-Markovian fluids in confined geometries.  more » « less
Award ID(s):
2002313
PAR ID:
10481184
Author(s) / Creator(s):
;
Publisher / Repository:
American institute of physics (AIP)
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
160
Issue:
2
ISSN:
0021-9606
Page Range / eLocation ID:
1-11
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Inverse problems involving time-fractional differential equations have become increasingly important for modeling systems with memory-dependent dynamics, particularly in biotransport and viscoelastic materials. Despite their potential, these problems remain challenging due to issues of stability, non-uniqueness, and limited data availability. Recent advancements in Physics-Informed Neural Networks (PINNs) offer a data-efficient framework for solving such inverse problems, yet most implementations are restricted to integer-order derivatives. In this work, we develop a PINN-based framework tailored for inverse problems involving time-fractional derivatives. We consider two representative applications: anomalous diffusion and fractional viscoelasticity. Using both synthetic and experimental datasets, we infer key physical parameters including the generalized diffusion coefficient and the fractional derivative order in the diffusion model and the relaxation parameters in a fractional Maxwell model. Our approach incorporates a customized residual loss function scaled by the standard deviation of observed data to enhance robustness. Even under 25% Gaussian noise, our method recovers model parameters with relative errors below 10%. Additionally, the framework accurately predicts relaxation moduli in porcine tissue experiments, achieving similar error margins. These results demonstrate the framework’s effectiveness in learning fractional dynamics from noisy and sparse data, paving the way for broader applications in complex biological and mechanical systems. 
    more » « less
  2. Many  biological  macromolecules  rely  on  metal  ions  to  maintain  structural  integrity  and  control  their  regulatory  function. In biological fluids, detection and identification of metal ions requires sensitive analytical tools with clear readouts.  In  this  work,  we  sought  to  investigate  the  potential  of  solution  Nuclear  Magnetic  Resonance  (NMR)  spectroscopy to analyze metal ion solutions and mixtures. To enable 1H NMR detection, we prepared the complexes of  eight  metal  ions  with  the  chelating  agent,  1,2-bis(o-aminophenoxy)ethane-N,N,N′,N′-tetraacetic  acid  (BAPTA).  The 1H NMR spectra were collected for BAPTA samples as a function of metal ion concentrations. The analysis of NMR data revealed that all metal ions with a notable exception of Mg2+ bind BAPTA with high affinities and form complexes with 1:1 metal-to-chelator stoichiometry. Both methylene and aromatic regions of the BAPTA 1H NMR spectra  experience  significant  changes  upon  the  metal  ion  complex  formation.  We  identified  the  spectroscopic  signatures  of  trivalent  and  paramagnetic  ions  and  demonstrated  that  the  binary  Zn2+/Pb2+  metal  ion  mixture  can  be  successfully analyzed by NMR. We conclude that complexation with BAPTA followed by the 1H NMR analysis is a sensitive method to detect and identify both nutritive and xenobiotic metal ions. 
    more » « less
  3. This work presents a straightforward computational method to estimate the rotational diffusion coefficient (Dr) of cells and particles of various sizes using the continuum fluid mechanics theory. We calculate the torque (Γ) for cells and particles immersed in fluids to find the mobility coefficient μ and then obtain the Dr by substituting Γ in the Einstein relation. Geometries are constructed using triangular mesh, and the model is solved with computational fluid dynamics techniques. This method is less intensive and more efficient than the widely used models. We simulate eight different particle geometries and compare the results with previous literature. 
    more » « less
  4. Abstract AlphaFold2 has revolutionized protein structure prediction from amino‐acid sequence. In addition to protein structures, high‐resolution dynamics information about various protein regions is important for understanding protein function. Although AlphaFold2 has neither been designed nor trained to predict protein dynamics, it is shown here how the information returned by AlphaFold2 can be used to predict dynamic protein regions at the individual residue level. The approach, which is termed cdsAF2, uses the 3D protein structure returned by AlphaFold2 to predict backbone NMR NHS2order parameters using a local contact model that takes into account the contacts made by each peptide plane along the backbone with its environment. By combining for each residue AlphaFold2's pLDDT confidence score for the structure prediction accuracy with the predictedS2value using the local contact model, an estimator is obtained that semi‐quantitatively captures many of the dynamics features observed in experimental backbone NMR NHS2order parameter profiles. The method is demonstrated for a set nine proteins of different sizes and variable amounts of dynamics and disorder. 
    more » « less
  5. Abstract BackgroundThe biophysics of an organism span multiple scales from subcellular to organismal and include processes characterized by spatial properties, such as the diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales, through the generation of representative models. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent the spatial and stochastic features of a biological system, limiting their insights and applications. However, spatial models describing biological systems with spatial information are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them and highlights the need for simpler methods able to model the spatial features of biological systems. ResultsIn this work, we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models, the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters. We additionally investigate objects and processes implicitly represented by ODE model terms and parameters and improve the reproducibility of spatial, stochastic models. ConclusionWe developed and demonstrate a method for generating spatiotemporal, multicellular models from non-spatial population dynamics models of multicellular systems. We envision employing our method to generate new ODE model terms from spatiotemporal and multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes. 
    more » « less