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: Modeling and analysis of ensemble average solvation energy and solute–solvent interfacial fluctuations
Abstract Variational implicit solvation models (VISMs) have gained extensive popularity in the molecular-level solvation analysis of biological systems due to their cost-effectiveness and satisfactory accuracy. Central in the construction of VISM is an interface separating the solute and the solvent. However, traditional sharp-interface VISMs fall short in adequately representing the inherent randomness of the solute–solvent interface, a consequence of thermodynamic fluctuations within the solute–solvent system. Given that experimentally observable quantities are ensemble averaged, the computation of the ensemble average solvation energy (EASE)–the averaged solvation energy across all thermodynamic microscopic states–emerges as a key metric for reflecting thermodynamic fluctuations during solvation processes. This study introduces a novel approach to calculating the EASE. We devise two diffuse-interface VISMs: one within the classic Poisson–Boltzmann (PB) framework and another within the framework of size-modified PB theory, accounting for the finite-size effects. The construction of these models relies on a new diffuse interface definition u ( x ) u\left(x), which represents the probability of a point x xfound in the solute phase among all microstates. Drawing upon principles of statistical mechanics and geometric measure theory, we rigorously demonstrate that the proposed models effectively capture EASE during the solvation process. Moreover, preliminary analyses indicate that the size-modified EASE functional surpasses its counterpart based on the classic PB theory across various analytic aspects. Our work is the first step toward calculating EASE through the utilization of diffuse-interface VISM.  more » « less
Award ID(s):
2306991 2110914 2306992
PAR ID:
10610755
Author(s) / Creator(s):
; ;
Publisher / Repository:
De Gruyter
Date Published:
Journal Name:
Computational and Mathematical Biophysics
Volume:
12
Issue:
1
ISSN:
2544-7297
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We extend the Calderón–Zygmund theory for nonlocal equations tostrongly coupled system of linear nonlocal equations A s u = f {\mathcal{L}^{s}_{A}u=f}, where the operator A s {\mathcal{L}^{s}_{A}}is formally given by A s u = n A ( x , y ) | x - y | n + 2 s ( x - y ) ( x - y ) | x - y | 2 ( u ( x ) - u ( y ) ) 𝑑 y . \mathcal{L}^{s}_{A}u=\int_{\mathbb{R}^{n}}\frac{A(x,y)}{|x-y|^{n+2s}}\frac{(x-%y)\otimes(x-y)}{|x-y|^{2}}(u(x)-u(y))\,dy. For 0 < s < 1 {0<1}and A : n × n {A:\mathbb{R}^{n}\times\mathbb{R}^{n}\to\mathbb{R}}taken to be symmetric and serving asa variable coefficient for the operator, the system under consideration is the fractional version of the classical Navier–Lamé linearized elasticity system. The study of the coupled system of nonlocal equations is motivated by its appearance in nonlocal mechanics, primarily in peridynamics. Our regularity result states that if A ( , y ) {A(\,\cdot\,,y)}is uniformly Holder continuous and inf x n A ( x , x ) > 0 {\inf_{x\in\mathbb{R}^{n}}A(x,x)>0}, then for f L loc p {f\in L^{p}_{\rm loc}}, for p 2 {p\geq 2}, the solution vector u H loc 2 s - δ , p {u\in H^{2s-\delta,p}_{\rm loc}}for some δ ( 0 , s ) {\delta\in(0,s)}. 
    more » « less
  2. Abstract Let 𝑋 be a Kähler manifold with semiample canonical bundle K X K_{X}.It is proved in [W. Jian, Y. Shi and J. Song, A remark on constant scalar curvature Kähler metrics on minimal models,Proc. Amer. Math. Soc.147(2019), 8, 3507–3513] that, for any Kähler class 𝛾, there exists δ > 0 \delta>0such that, for all t ( 0 , δ ) t\in(0,\delta), there exists a unique cscK metric g t g_{t}in K X + t γ K_{X}+t\gamma.In this paper, we prove that { ( X , g t ) } t ( 0 , δ ) \{(X,g_{t})\}_{t\in(0,\delta)}have uniformly bounded Kähler potentials, volume forms and diameters.As a consequence, these metric spaces are pre-compact in the Gromov–Hausdorff sense. 
    more » « less
  3. Abstract We introduce a distributional Jacobian determinant det D V β ( D v ) \det DV_{\beta}(Dv)in dimension two for the nonlinear complex gradient V β ( D v ) = | D v | β ( v x 1 , v x 2 ) V_{\beta}(Dv)=\lvert Dv\rvert^{\beta}(v_{x_{1}},-v_{x_{2}})for any β > 1 \beta>-1, whenever v W loc 1 , 2 v\in W^{1\smash{,}2}_{\mathrm{loc}}and β | D v | 1 + β W loc 1 , 2 \beta\lvert Dv\rvert^{1+\beta}\in W^{1\smash{,}2}_{\mathrm{loc}}.This is new when β 0 \beta\neq 0.Given any planar ∞-harmonic function 𝑢, we show that such distributional Jacobian determinant det D V β ( D u ) \det DV_{\beta}(Du)is a nonnegative Radon measure with some quantitative local lower and upper bounds.We also give the following two applications. Applying this result with β = 0 \beta=0, we develop an approach to build up a Liouville theorem, which improves that of Savin.Precisely, if 𝑢 is an ∞-harmonic function in the whole R 2 \mathbb{R}^{2}with lim inf R inf c R 1 R B ( 0 , R ) | u ( x ) c | d x < , \liminf_{R\to\infty}\inf_{c\in\mathbb{R}}\frac{1}{R}\barint_{B(0,R)}\lvert u(x)-c\rvert\,dx<\infty,then u = b + a x u=b+a\cdot xfor some b R b\in\mathbb{R}and a R 2 a\in\mathbb{R}^{2}.Denoting by u p u_{p}the 𝑝-harmonic function having the same nonconstant boundary condition as 𝑢, we show that det D V β ( D u p ) det D V β ( D u ) \det DV_{\beta}(Du_{p})\to\det DV_{\beta}(Du)as p p\to\inftyin the weak-⋆ sense in the space of Radon measure.Recall that V β ( D u p ) V_{\beta}(Du_{p})is always quasiregular mappings, but V β ( D u ) V_{\beta}(Du)is not in general. 
    more » « less
  4. Abstract Let K / Q p K/\mathbf{Q}_{p}be unramified.Inside the Emerton–Gee stack X 2 \mathcal{X}_{2}, one can consider the locus of two-dimensional mod 𝑝 representations of Gal ( K ̄ / K ) \mathrm{Gal}(\overline{K}/K)having a crystalline lift with specified Hodge–Tate weights.We study the case where the Hodge–Tate weights are irregular, which is an analogue for Galois representations of the partial weight one condition for Hilbert modular forms.We prove that if the gap between each pair of weights is bounded by 𝑝 (the irregular analogue of a Serre weight), then this locus is irreducible.We also establish various inclusion relations between these loci. 
    more » « less
  5. Abstract LetXbe acompact orientable non-Haken 3-manifold modeled on the Thurston geometry Nil {\operatorname{Nil}}. We show that the diffeomorphism group Diff ( X ) {\operatorname{Diff}(X)}deformation retracts to the isometry group Isom ( X ) {\operatorname{Isom}(X)}. Combining this with earlier work by many authors, this completes the determination the homotopy type of Diff ( X ) {\operatorname{Diff}(X)}for any compact, orientable, prime 3-manifoldX. 
    more » « less