Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics simulation”

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 OSTI2021

Solvation Structure and Dynamics of Mg(TFSI) 2 Aqueous Electrolyte

Using ab initio molecular dynamics (AIMD) simulations, classical molecular dynamics (CMD) simulations, small-angle X-ray scattering (SAXS), and pulsed-field gradient nuclear magnetic resonance (PFG-NMR), the solvation structure and ion dynamics of magnesium bis(trifluoromethanesulfonyl)imide (Mg(TFSI) 2 ) aqueous electrolyte at 1, 2, and 3 m concentrations are investigated. From AIMD and CMD simulations, the first solvation shell of an Mg 2+ ion is found to be composed of six water molecules in an octahedral configuration and the solvation shell is rather rigid. The TFSI - ions prefer to stay in the second solvation shell and beyond. Meanwhile, the comparable diffusion coefficients of positive and negative ions in Mg(TFSI) 2 aqueous electrolytes have been observed, which is mainly due to the formation of the stable [Mg(H 2 O) 6 ] 2+ complex, and, as a result, the increased effective Mg ion size. Finally, the calculated correlated transference numbers are lower than the uncorrelated ones even at the low concentration of 2 and 3 m, suggesting the enhanced correlations between ions in the multivalent electrolytes. This work provides a molecular-level understanding of how the solvation structure and multivalency of the ion affect the dynamics and transport properties of the multivalent electrolyte, providing insight for rational designs of electrolytes for improved ion transport properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Superhydrophilicity of $α$-alumina surfaces results from tight binding of interfacial waters to specific aluminols

Understanding the microscopic driving force of water wetting is challenging and important for design of materials. The relations between structure, dynamics and hydrogen bonds of interfacial water can be investigated using molecular dynamics simulations. Here, contact angles at the alumina (0001) and ($11\bar{2}0$) surfaces are studied using both classical molecular dynamics simulations and experiments. To test the superhydrophilicity, the free energy cost of removing waters near the interfaces are calculated using the density fluctuations method. The strength of hydrogen bonds is determined by their lifetime and geometry. Both surfaces are superhydrophilic and the (0001) surface is more hydrophilic. Interactions between surfaces and interfacial waters promote a templating effect whereby the latter are aligned in a pattern that follows the underlying lattice of the surfaces. Translational and rotational dynamics of interfacial water molecules are slower than in bulk water. Hydrogen bonds between water and both surfaces are asymmetric, water-to-aluminol ones are stronger than aluminol-to-water ones. Molecular dynamics simulations eliminate the impacts of surface contamination when measuring contact angles and the results reveal the microscopic origin of the macroscopic superhydrophilicity of alumina surfaces: strong water-to-aluminol hydrogen bonds.

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

Depletion-Driven Assembly of Polymer-Coated Nanocrystals

Depletion-driven assembly has been widely studied for micron-sized colloids, but questions remain at the nanoscale where the governing physics are impacted by stabilizing surface ligands or wrapping polymers, whose dimensions are similar to those of the colloidal core and the depletant. Here, we probe how wrapping colloidal tin-doped indium oxide nanocrystals with polymers affects their depletion-induced interactions and assembly in solutions of polyethylene glycol. Copolymers of polyacrylic acid grafted with polyethylene oxide provide nanocrystal wrappings with different effective polymer graft densities and molecular weights. (Ultra) small angle X-ray scattering, coarse-grained molecular dynamics simulation, and molecular thermodynamic theory were combined to analyze how depletant size and polymer wrapping characteristics affect depletion interactions, structure, and phase behavior. The results show how depletant molecular weight, as well as surface density and molecular weight of polymer grafts, set thresholds for assembly. These signatures are unique to depletion-driven assembly of nanoscale colloids and mirror phase behaviors of grafted nanoparticle–polymer composites. Optical and rheological responses of depletion-driven assemblies of nanocrystals with different polymer shell architectures were probed to learn how their structural differences impact properties. Finally, we discuss how these wrinkles for depletion-driven assembly at the nanoscale may provide fresh opportunities for designing responsive depletion interactions and dynamically reconfigurable materials.

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

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

Conformation and dynamics of ring polymers under symmetric thin film confinement

Understanding the structure and dynamics of polymers under confinement has been of widespread interest, and one class of polymers that have received comparatively little attention under confinement is that of ring polymers. The properties of non-concatenated ring polymers can also be important in biological fields because ring polymers have been proven to be a good model to study DNA organization in the cell nucleus. From our previous study, linear polymers in a cylindrically confined polymer melt were found to segregate from each other as a result of the strong correlation hole effect that is enhanced by the confining surfaces. By comparison, our subsequent study of linear polymers in confined thin films at similar levels of confinements found only the onset of segregation. Here, we use molecular dynamics simulation to investigate the chain conformations and dynamics of ring polymers under planar (1D) confinement as a function of film thickness. Our results show that conformations of ring polymers are similar to the linear polymers under planar confinement, except that ring polymers are less compressed in the direction normal to the walls. While we find that the correlation hole effect is enhanced under confinement, it is not as pronounced as the linear polymers under 2D confinement. Finally, we show that chain dynamics far above Tg are primarily affected by the friction from walls based on the monomeric friction coefficient we get from the Rouse mode analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
DOE OSTI2020

In Silico Insights into Protein–Protein Interaction Disruptive Mutations in the PCSK9-LDLR Complex

Gain-of-function mutations in PCSK9 (proprotein convertase subtilisin/kexin type 9) lead to reduced uptake of LDL (low density lipoprotein) cholesterol and, therefore, increased plasma LDL levels. However, the mechanism by which these mutants reduce LDL reuptake is not fully understood. Here, we have used molecular dynamics simulations, MM/PBSA (Molecular Mechanics/Poisson–Boltzmann Surface Area) binding affinity calculations, and residue interaction networks, to investigate the protein–protein interaction (PPI) disruptive effects of two of PCSK9's gain-of-function mutations, Ser127Arg and Asp374Tyr on the PCSK9 and LDL receptor complex. In addition to these PPI disruptive mutants, a third, non-interface mutation (Arg496Trp) is included as a positive control. Our results indicate that Ser127Arg and Asp374Tyr confer significantly improved binding affinity, as well as different binding modes, when compared to the wild-type. These PPI disruptive mutations lie between the EGF(A) (epidermal growth factor precursor homology domain A) of the LDL receptor and the catalytic domain of PCSK9 (Asp374Tyr) and between the prodomain of PCSK9 and the β-propeller of the LDL receptor (Ser127Arg). The interactions involved in these two interfaces result in an LDL receptor that is sterically inhibited from entering its closed conformation. This could potentially implicate the prodomain as a target for small molecule inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗
DOE OSTI2025

Assessing the Effect of Explicit Polarizability on Models of Carbon Dioxide Solvation in Ionic Liquids

Ionic liquids are an important possible carbon capture material because of their anomalously high sorption selectivity for carbon dioxide over other gases common in air. Many research groups have investigated the molecular origins of this property and provided important insights, including using 1D and 2D-IR spectroscopy. Molecular dynamics simulations have been indispensable to the interpretation of these experiments. In prior molecular dynamics simulation work, charge-scaled force fields have typically been used to provide a mean-field treatment of effects vital to ionic liquid systems such as charge transfer and polarization. Here, we compare models of carbon dioxide solvated in ionic liquids with explicit polarization to models of the same with implicit polarizability through charge-scaling. We calculate structural, dynamical, and spectroscopic properties, and make comparisons to the same items measured in experiment. In this study, we focus on two ionic liquids: 1-butyl-3- methylimidazolium (BMIM + ) paired with bis(trifluoromethane sulfonyl imide) (Tf 2 N − ) and 1-butyl-3-methylimidazolium (BMIM+) paired with hexafluorophosphate (PF 6 − ). We find that many structural, dynamical, and spectroscopic properties are changed when polarization is modeled explicitly. We also find that explicit polarizability softens local ion cages around the carbon dioxide and that the long-time diffusion of the carbon dioxide is gated by the reorganization of the ionic liquid molecules. Comparisons to experiment show modest improvement of many observables compared with experiment for the explicitly polarizable model over the charge-scaled model. Overall, our results show that charge-scaled force fields are likely sufficient to compute spectroscopic properties of carbon dioxide in ionic liquids and suggest some interpretive rules for understanding their structural and dynamical properties. Those using charge-scaled force fields should generally assume that the ion cages around solutes such as carbon dioxide are too stiff and cation-rich in their models and adjust their interpretations and predictions accordingly.

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

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Molecular design and fabrication of PIM-1/polyphosphazene blend membranes with high performance for CO 2 /N 2 separation

New polymeric blend membranes for CO 2 separation were synthesized based on insights from molecular dynamics simulation. A molecular-level structure-property relationship in polymers of intrinsic microporosity (PIM) based blend membranes, was investigated in detail computationally. Calculated local density profiles and energy of interaction of the blend membranes, composed of PIM-1 and various polyphosphazenes, showed that using the polyphosphazene with a higher concentration of ether side chains can improve the compatibility with PIM-1. Furthermore, based on the findings of computational studies, blend membranes were experimentally fabricated from PIM-1 and polyphosphazenes with various polyether side chain concentrations. Polyether concentration in polyphosphazenes was correlated with the film properties and gas transport performance of the blend membranes. Blend membranes showed very high CO 2 permeability (3100-5300 barrer) and improved CO 2 /N 2 selectivity (24-28), outperforming all other PIM-based blend membranes reported to date. Moreover, the CO 2 permeability performance of the blend membranes was tested 566 hours under real post-combustion flue gas from a coal-fired power plant, including CO 2 , N 2 , H 2 O, O 2 , SO x and NO x .

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

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Driving and characterizing nucleation of urea and glycine polymorphs in water

Crystal nucleation is relevant across the domains of fundamental and applied sciences. However, in many cases, its mechanism remains unclear due to a lack of temporal or spatial resolution. To gain insights into the molecular details of nucleation, some form of molecular dynamics simulations is typically performed; these simulations, in turn, are limited by their ability to run long enough to sample the nucleation event thoroughly. To overcome the timescale limits in typical molecular dynamics simulations in a manner free of prior human bias, here, we employ the machine learning-augmented molecular dynamics framework “reweighted autoencoded variational Bayes for enhanced sampling (RAVE).” We study two molecular systems—urea and glycine—in explicit all-atom water, due to their enrichment in polymorphic structures and common utility in commercial applications. From our simulations, we observe multiple back-and-forth nucleation events of different polymorphs from homogeneous solution; from these trajectories, we calculate the relative ranking of finite-sized polymorph crystals embedded in solution, in terms of the free-energy difference between the finite-sized crystal polymorph and the original solution state. We further observe that the obtained reaction coordinates and transitions are highly nonclassical.

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

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

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

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

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

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗