Alireza Foroozani Behbahani
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Alireza Foroozani Behbahani has a total of 24 co-authors. Here are the top co-authors by number of shared publications:
- Vagelis Harmandaris
University of Crete
9 shared publications - Manolis Doxastakis
University of Tennessee
5 shared publications - Craig Burkhart
Goodyear Tire & Rubber Company
3 shared publications - Friederike Schmid
Johannes Gutenberg University Mainz
3 shared publications Patrycja Polińska
Avenue Gordon Smith
3 shared publications- Wei Li
Institute for Chemical Reaction Design and Discovery, Hokkaido University, Sapporo, Japan
3 shared publications
Latest Publications
- Open Access
- 8 October 2026
- Published by American Chemical Society (ACS) in Macromolecules
From Single-Chain Dynamics to Structure Formation: Dynamic Self-Consistent Field Theory and Molecular Dynamics of (Co)polymer Melts across Entanglement Regimes
Dynamic self-consistent field theory (DSCFT) provides an efficient continuum framework for studying structure formation in inhomogeneous polymer systems, but its predictive accuracy depends on the choice of nonlocal mobilities. Here, we construct mobility functions for moderately and strongly entangled homopolymer and diblock copolymer systems from the relaxation dynamics of single-chain structure factors, based on molecular dynamics (MD) simulations of the Kremer–Grest model and analytical reptation theory. Single-chain mobilities are combined such that the resulting DSCFT accounts for the dependence of fluxes on local chain densities. The theory is then applied to the spinodal decomposition of symmetric homopolymer blends and diblock copolymer melts following a quench into the (micro)phase-separation regime. Predictions of DSCFT are systematically compared with MD simulations. Mobility functions derived from single-chain dynamics are found to reproduce the kinetics of structure formation more accurately than conventional Debye-type mobilities. We additionally investigate the influence of adding stochastic currents (noise) that are correlated according to the fluctuation–dissipation relation. At low noise levels, they enable the generation of equilibrium initial states and facilitate defect annealing. At high noise levels, however, nonlinear effects lead to discrepancies between DSCFT and MD simulations.
- 8 October 2026
- Published by American Chemical Society (ACS) in Macromolecules
Open AccessFrom Single-Chain Dynamics to Structure Formation: Dynamic Self-Consistent Field Theory and Molecular Dynamics of (Co)polymer Melts across Entanglement Regimes
Dynamic self-consistent field theory (DSCFT) provides an efficient continuum framework for studying structure formation in inhomogeneous polymer systems, but its predictive accuracy depends on the choice of nonlocal mobilities. Here, we construct mobility functions for moderately and strongly entangled homopolymer and diblock copolymer systems from the relaxation dynamics of single-chain structure factors, based on molecular dynamics (MD) simulations of the Kremer–Grest model and analytical reptation theory. Single-chain mobilities are combined such that the resulting DSCFT accounts for the dependence of fluxes on local chain densities. The theory is then applied to the spinodal decomposition of symmetric homopolymer blends and diblock copolymer melts following a quench into the (micro)phase-separation regime. Predictions of DSCFT are systematically compared with MD simulations. Mobility functions derived from single-chain dynamics are found to reproduce the kinetics of structure formation more accurately than conventional Debye-type mobilities. We additionally investigate the influence of adding stochastic currents (noise) that are correlated according to the fluctuation–dissipation relation. At low noise levels, they enable the generation of equilibrium initial states and facilitate defect annealing. At high noise levels, however, nonlinear effects lead to discrepancies between DSCFT and MD simulations.
- Open Access
- 6 March 2026
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 59 (6) , 3455-3464
Stress Relaxation in Monodisperse Entangled Polymer Melts: Correlation between Viscoelastic Response and Single-Chain Relaxation via Molecular Dynamics Simulations
- 6 March 2026
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 59 (6) , 3455-3464
Open AccessStress Relaxation in Monodisperse Entangled Polymer Melts: Correlation between Viscoelastic Response and Single-Chain Relaxation via Molecular Dynamics Simulations
- 5 March 2026
- Published by IOP Publishing in Physica Scripta
- Vol. 101 (9) , 095911
Percolation conductivity in composites with a binary blend matrix and a non-uniform distribution of conductive sites: a lattice-based approach
Despite the widespread use of conductive composites with a binary blend matrix in electrical applications, computational studies examining the effects of conductive particle localization within individual phases and at their interface on conductance and the percolation threshold remain limited. To address this gap, we employ a simple cubic lattice model to simulate the localization of conductive sites within individual phases and at the interface between them. Our results show that exclusive localization of conductive sites within a single phase lowers the percolation threshold and enhances conductance, especially near the percolation threshold, whereas localization at the interface further reduces the threshold. When all conductive sites are exclusively localized in one phase, the percolation threshold can be expressed in two ways: first, as the fraction of conductive sites calculated relative to the total number of sites, denoted as Π; and second, as the fraction calculated relative to the available sites in that phase, denoted as Π f . In this situation, the second phase effectively acts as an excluded volume, reducing the number of available sites. At a constant concentration of the excluded-volume phase, increasing the size of this phase causes the percolation threshold to first decrease and then converge to a constant value. The constant value corresponds to a state in which the structure of conductive sites within the available volume is similar to that of a homogeneous system. In this case, the percolation threshold expressed relative to the available sites is given by Π f = Π homogen , which equals 0.31 in the case of a simple cubic lattice; additionally, the percolation threshold relative to all sites is given by Π = (1 − φ A )Π homogen , where φ A is the fraction of the excluded-volume phase. We also studied situations in which particles are distributed between the two phases with different probabilities. In such cases, the percolation threshold lies between the two extreme cases: a fully homogeneous system and a fully segregated system in which one phase acts as an excluded volume.
- 5 March 2026
- Published by IOP Publishing in Physica Scripta
- Vol. 101 (9) , 095911
Percolation conductivity in composites with a binary blend matrix and a non-uniform distribution of conductive sites: a lattice-based approach
Despite the widespread use of conductive composites with a binary blend matrix in electrical applications, computational studies examining the effects of conductive particle localization within individual phases and at their interface on conductance and the percolation threshold remain limited. To address this gap, we employ a simple cubic lattice model to simulate the localization of conductive sites within individual phases and at the interface between them. Our results show that exclusive localization of conductive sites within a single phase lowers the percolation threshold and enhances conductance, especially near the percolation threshold, whereas localization at the interface further reduces the threshold. When all conductive sites are exclusively localized in one phase, the percolation threshold can be expressed in two ways: first, as the fraction of conductive sites calculated relative to the total number of sites, denoted as Π; and second, as the fraction calculated relative to the available sites in that phase, denoted as Π f . In this situation, the second phase effectively acts as an excluded volume, reducing the number of available sites. At a constant concentration of the excluded-volume phase, increasing the size of this phase causes the percolation threshold to first decrease and then converge to a constant value. The constant value corresponds to a state in which the structure of conductive sites within the available volume is similar to that of a homogeneous system. In this case, the percolation threshold expressed relative to the available sites is given by Π f = Π homogen , which equals 0.31 in the case of a simple cubic lattice; additionally, the percolation threshold relative to all sites is given by Π = (1 − φ A )Π homogen , where φ A is the fraction of the excluded-volume phase. We also studied situations in which particles are distributed between the two phases with different probabilities. In such cases, the percolation threshold lies between the two extreme cases: a fully homogeneous system and a fully segregated system in which one phase acts as an excluded volume.
- Open Access
- 17 December 2024
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 58 (1) , 767-786
Relaxation Dynamics of Entangled Linear Polymer Melts via Molecular Dynamics Simulations
- 17 December 2024
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 58 (1) , 767-786
Open AccessRelaxation Dynamics of Entangled Linear Polymer Melts via Molecular Dynamics Simulations
- Open Access
Relaxation Dynamics of Entangled Linear Polymer Melts via Molecular Dynamics Simulations
We present an extensive analysis of the relaxation dynamics of entangled linear polymer melts via long-time molecular dynamics simulations of a generic bead-spring model. We study the mean-squared displacements, the autocorrelation function of the end-to-end vector, $P(t)$, the single-chain dynamic structure factor, $S(q,t)$, and the linear viscoelastic properties, especially the shear stress relaxation modulus, $G(t)$. The simulation data are compared with the theoretically expected scaling laws for different time regimes of entangled melts, and with analytical expressions that account for different relaxation mechanisms in the tube model, namely, reptation, contour length fluctuation (CLF), and constraint release (CR). CLF involves a $t^{1/4}$ scaling regime in the time-dependence of $(1-P(t))$. With increasing chain length, a gradual development of this scaling regime is observed. In the absence of CR, the tube model further predicts that at long times, the chain dynamics is governed by one central quantity, the ``surviving tube fraction'' $μ(t)$. As a result, one expects $S(q,t) \propto G(t) \propto P(t)$ in that time regime. We test this prediction by comparing $S(q,t)$ and $G(t)$ with $P(t)$. For both quantities, proportionality with $P(t)$ is not observed, indicating that CR has an important effect on the relaxation of these two quantities. Instead, to a very good approximation, we find $G(t)\propto P(t)^{2}$ at late times, which is consistent with the dynamic tube dilation or double reptation approximations for the CR process. In addition, we calculate non-local mobility functions, which can be used in dynamic density functional theories for entangled inhomogeneous polymer blends, and discuss the effect of entanglements on the shape of these functions.Open AccessRelaxation Dynamics of Entangled Linear Polymer Melts via Molecular Dynamics Simulations
We present an extensive analysis of the relaxation dynamics of entangled linear polymer melts via long-time molecular dynamics simulations of a generic bead-spring model. We study the mean-squared displacements, the autocorrelation function of the end-to-end vector, $P(t)$, the single-chain dynamic structure factor, $S(q,t)$, and the linear viscoelastic properties, especially the shear stress relaxation modulus, $G(t)$. The simulation data are compared with the theoretically expected scaling laws for different time regimes of entangled melts, and with analytical expressions that account for different relaxation mechanisms in the tube model, namely, reptation, contour length fluctuation (CLF), and constraint release (CR). CLF involves a $t^{1/4}$ scaling regime in the time-dependence of $(1-P(t))$. With increasing chain length, a gradual development of this scaling regime is observed. In the absence of CR, the tube model further predicts that at long times, the chain dynamics is governed by one central quantity, the ``surviving tube fraction'' $μ(t)$. As a result, one expects $S(q,t) \propto G(t) \propto P(t)$ in that time regime. We test this prediction by comparing $S(q,t)$ and $G(t)$ with $P(t)$. For both quantities, proportionality with $P(t)$ is not observed, indicating that CR has an important effect on the relaxation of these two quantities. Instead, to a very good approximation, we find $G(t)\propto P(t)^{2}$ at late times, which is consistent with the dynamic tube dilation or double reptation approximations for the CR process. In addition, we calculate non-local mobility functions, which can be used in dynamic density functional theories for entangled inhomogeneous polymer blends, and discuss the effect of entanglements on the shape of these functions. - 3 October 2024
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 161 (13)
Relaxation dynamics of a liquid in the vicinity of an attractive surface: The process of escaping from the surface
We analyze the displacements of the particles of a glass-forming molecular liquid perpendicular to a confining solid surface using extensive molecular dynamics simulations with atomistic models. In the vicinity of an attractive surface, the liquid molecules are trapped. Transient localization of liquid molecules near the surface introduces a relaxation process related to the escape of molecules from the surface into the dynamics of the interfacial liquid layer. To describe this process, we analyze several dynamical observables of the confined liquid. The self-intermediate scattering function and the mean-squared displacement of the particles located in the interfacial layer are dominated by the process of escaping from the surface. This relaxation process is also associated with a strong heterogeneity in the mobility of the interfacial particles. The studied model liquid is hydrogenated methyl methacrylate. For the confining wall, we consider different models, namely a periodic single layer of graphene and a frozen amorphous configuration of the bulk liquid (frozen wall). Near graphene, where the liquid molecules form a layered structure and adopt parallel-to-surface orientation, a clear separation between small-scale movements of the molecules near the surface and the process of escaping from the surface is observed. This is reflected in the three-step relaxation of the interfacial layer. However, near the frozen wall, where the liquid molecules do not have a preferential alignment, a clear three-step relaxation is not seen, even though the dynamical quantities are controlled by the process of escaping from the surface.
- 3 October 2024
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 161 (13)
Relaxation dynamics of a liquid in the vicinity of an attractive surface: The process of escaping from the surface
We analyze the displacements of the particles of a glass-forming molecular liquid perpendicular to a confining solid surface using extensive molecular dynamics simulations with atomistic models. In the vicinity of an attractive surface, the liquid molecules are trapped. Transient localization of liquid molecules near the surface introduces a relaxation process related to the escape of molecules from the surface into the dynamics of the interfacial liquid layer. To describe this process, we analyze several dynamical observables of the confined liquid. The self-intermediate scattering function and the mean-squared displacement of the particles located in the interfacial layer are dominated by the process of escaping from the surface. This relaxation process is also associated with a strong heterogeneity in the mobility of the interfacial particles. The studied model liquid is hydrogenated methyl methacrylate. For the confining wall, we consider different models, namely a periodic single layer of graphene and a frozen amorphous configuration of the bulk liquid (frozen wall). Near graphene, where the liquid molecules form a layered structure and adopt parallel-to-surface orientation, a clear separation between small-scale movements of the molecules near the surface and the process of escaping from the surface is observed. This is reflected in the three-step relaxation of the interfacial layer. However, near the frozen wall, where the liquid molecules do not have a preferential alignment, a clear three-step relaxation is not seen, even though the dynamical quantities are controlled by the process of escaping from the surface.
- Open Access
- 13 May 2024
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 6 (2) , 023161
Local viscoelastic properties and shear stress propagation in bulk and confined polymer melts and low-molecular weight liquids
Through the analysis of the spatial correlations of local stress, we detect the propagation of long-ranged liquid-elasticity-mediated shear stress waves in polymeric and low-molecular weight liquids. The propagation of shear waves is effectively planar; i.e., propagates in the plane. The autocorrelation functions of the local stress of a region are affected by both the relaxation of stress in that region and the propagation of stress from the region to the rest of the sample. However, due to the planar propagation of shear waves, the transfer of from those slices of the simulation box that are periodic in the plane is negligible. This allows direct probing of the position-dependent local stress relaxation modulus of liquid in the vicinity of a confining surface.- 13 May 2024
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 6 (2) , 023161
Open AccessLocal viscoelastic properties and shear stress propagation in bulk and confined polymer melts and low-molecular weight liquids
Through the analysis of the spatial correlations of local stress, we detect the propagation of long-ranged liquid-elasticity-mediated shear stress waves in polymeric and low-molecular weight liquids. The propagation of shear waves is effectively planar; i.e., propagates in the plane. The autocorrelation functions of the local stress of a region are affected by both the relaxation of stress in that region and the propagation of stress from the region to the rest of the sample. However, due to the planar propagation of shear waves, the transfer of from those slices of the simulation box that are periodic in the plane is negligible. This allows direct probing of the position-dependent local stress relaxation modulus of liquid in the vicinity of a confining surface. - Open Access
Correction: Foroozani Behbahani, A.; Harmandaris, V. Gradient of Segmental Dynamics in Stereoregular Poly(methyl methacrylate) Melts Confined between Pristine or Oxidized Graphene Sheets. Polymers 2021, 13, 830
The authors wish to make the following corrections to this paper:
Open AccessCorrection: Foroozani Behbahani, A.; Harmandaris, V. Gradient of Segmental Dynamics in Stereoregular Poly(methyl methacrylate) Melts Confined between Pristine or Oxidized Graphene Sheets. Polymers 2021, 13, 830
The authors wish to make the following corrections to this paper:
- 17 September 2021
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 54 (18) , 8693-8713
Dynamics of Long Entangled Polyisoprene MeltsviaMultiscale Modeling
- 17 September 2021
- Published by American Chemical Society (ACS) in Macromolecules
- Vol. 54 (18) , 8693-8713
Dynamics of Long Entangled Polyisoprene MeltsviaMultiscale Modeling