PubMed Health⌕ Search

Biomedical subjects

F Bresme

Publications and source records attributed to F Bresme.

6 recordsLinked to original sources
PubMed2006

Algorithms for the computation of solutions of the Ornstein-Zernike equation.

We introduce a robust and efficient methodology to solve the Ornstein-Zernike integral equation using the pseudoarc length (PAL) continuation method that reformulates the integral equation in an equivalent but nonstandard form. This enables the computation of solutions in regions where the compressibility experiences large changes or where the existence of multiple solutions and so-called branch points prevents Newton's method from converging. We illustrate the use of the algorithm with a difficult problem that arises in the numerical solution of integral equations, namely the evaluation of the so-called no-solution line of the Ornstein-Zernike hypernetted chain (HNC) integral equation for the Lennard-Jones potential. We are able to use the PAL algorithm to solve the integral equation along this line and to connect physical and nonphysical solution branches (both isotherms and isochores) where appropriate. We also show that PAL continuation can compute solutions within the no-solution region that cannot be computed when Newton and Picard methods are applied directly to the integral equation. While many solutions that we find are new, some correspond to states with negative compressibility and consequently are not physical.

Journal Article↗
PubMed2006

Interactions of polarizable media in water: a molecular dynamics approach.

We investigate the interactions of polarizable solutes in water as a function of the solute permittivity. A generic and computationally efficient simulation methodology for the investigation of systems involving dielectric discontinuities is introduced. We report results for interactions between two polarizable cylindrical solutes of nanometer dimensions, which demonstrate that the interactions between the solutes strongly depend on the solute permittivity epsilon. For low permittivity, epsilon approximately 1-2, the interactions are dominated by surface tension forces whose origin lies in the formation of a vapor cavity between the two hydrophobic solutes. This effect leads to a drying transition, where the intersolute force changes discontinuously at a specific solute-solute separation. We find that a moderate permittivity, epsilon approximately 20, enhances the solvation of the polarizable objects inhibiting this drying transition. In the limit of moderately high permittivity, the interactions are dominated by solvation forces. These forces are much larger than those calculated using macroscopic models of dielectrics, which consider water as a continuum dielectric medium. Our results emphasize the importance of including the solvent explicitly to investigate dielectric discontinuities and interactions between polarizable media in water.

Biopolymers↗
PubMed2000

Order-disorder transition in the solid phase of a charged hard sphere model.

We investigate the solid phases of the restricted primitive model (RPM). Monte Carlo simulations show the existence of an order-disorder transition from a substitutionally disordered face centered cubic lattice (fcc) to a new ordered fcc structure which is proposed as the ground state of the RPM at the close packing density. Our results suggest that the new phase might turn out in a new triple point in the RPM phase diagram involving three solid phases: CsCl, fcc ordered and fcc disordered structures. The order-disorder transition is also studied using the cell theory. The theory shows good agreement with the simulation results and suggests that the transition is weakly first order.

Journal Article↗