Charlie Ray Wand
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
Charlie Ray Wand has a total of 28 co-authors. Here are the top co-authors by number of shared publications:
- Flor R. Siperstein
The University of Manchester
4 shared publications - Paola Carbone
Department of Chemistry, University of Manchester, Oxford Rd, Manchester, UK
3 shared publications - Daan Frenkel
Yusuf Hamied Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW
2 shared publications - Kieran Harris
The University of Manchester
2 shared publications - Maria Panoukidou
School of Physics and Astronomy
2 shared publications Peter Visser
AkzoNobel
2 shared publications
Latest Publications
- Open Access
- 13 January 2026
- Published by Taylor & Francis in Molecular Simulation
- Vol. 52 (2) , 120-133
A comparison of molecular dynamics chemical reaction heuristics applied to crosslink formation
- 13 January 2026
- Published by Taylor & Francis in Molecular Simulation
- Vol. 52 (2) , 120-133
Open AccessA comparison of molecular dynamics chemical reaction heuristics applied to crosslink formation
- Open Access
- 17 April 2024
- Published by Royal Society of Chemistry (RSC) in RSC Applied Interfaces
- Vol. 1 (4) , 812-820
Segregation in epoxy/amine systems on iron oxide surfaces
Preferential adsorption of amine crosslinkers leads to composition heterogeneity in coatings which is strongly influenced by the surface structure.
- 17 April 2024
- Published by Royal Society of Chemistry (RSC) in RSC Applied Interfaces
- Vol. 1 (4) , 812-820
Open AccessSegregation in epoxy/amine systems on iron oxide surfaces
Preferential adsorption of amine crosslinkers leads to composition heterogeneity in coatings which is strongly influenced by the surface structure.
- Open Access
- 1 October 2022
- Published by Wiley in Chemistry – A European Journal
- Vol. 28 (68)
Water‐Mediated Epoxy/Surface Adhesion: Understanding the Interphase Region
Epoxy resins coatings are commonly found in corrosion protection coatings but the presence of water can affect their adhesion to the substrate, often weakening the adhesion of the coating to the solid, reducing its efficiency. Nevertheless, small amounts of water can enhance the epoxy/substrate interactions. In this work, the interphase region of an epoxy precursor and metal oxide substrates is investigated using molecular simulations and it is found that water accumulates between the epoxy layer and the solid substrate. At high water concentrations (9 wt %) the interaction between the epoxy precursor and the solid surface is weakened regardless of the nature of the solid, but at low water concentrations the nature of the solid surface becomes important. For hematite, the presence of water decreases the strength of adhesion but for goethite the presence of a small amount of water (3 wt %) enhances the adhesion to the surface resulting in a densification at the interface.
- 1 October 2022
- Published by Wiley in Chemistry – A European Journal
- Vol. 28 (68)
Open AccessWater‐Mediated Epoxy/Surface Adhesion: Understanding the Interphase Region
Epoxy resins coatings are commonly found in corrosion protection coatings but the presence of water can affect their adhesion to the substrate, often weakening the adhesion of the coating to the solid, reducing its efficiency. Nevertheless, small amounts of water can enhance the epoxy/substrate interactions. In this work, the interphase region of an epoxy precursor and metal oxide substrates is investigated using molecular simulations and it is found that water accumulates between the epoxy layer and the solid substrate. At high water concentrations (9 wt %) the interaction between the epoxy precursor and the solid surface is weakened regardless of the nature of the solid, but at low water concentrations the nature of the solid surface becomes important. For hematite, the presence of water decreases the strength of adhesion but for goethite the presence of a small amount of water (3 wt %) enhances the adhesion to the surface resulting in a densification at the interface.
- Open Access
- 1 January 2022
- Published by Elsevier in SSRN Electronic Journal
Molecular Control of Epoxy-Amine / Iron Oxide Interphase Properties
Understanding the polymer / metal oxide interphase is central to improving the performance lifetime of corrosion resistant coatings, where network polymers commonly form via step growth mechanisms in the presence of pigments. Here we employ a holistic analytical approach encompassing ATR-FTIR, DSC and molecular dynamics simulations to consider how crosslinker structure affects adsorption and incorporation into the network, using a stoichiometric mixture of DGEBA with m-xylylenediamine (MXDA) cured in the presence of hematite (Fe 2 O 3 ) and goethite (FeOOH) powders. We find that the rigid MXDA molecule has two distinct binding modes on both hematite and goethite, and that synergistic hydrogen bonding modes observed on goethite limit interconversion between the two. Moreover, we find that binding persists in fully cured composite samples, determining the levels of residual amine. In contrast to tetra-amine systems, however, we find that the Tg of composite specimens is independent of added hematite and goethite volumes. Molecular dynamics simulations demonstrate this is due to electrostatic binding between the cationic Fe sites and electronegative heteroatoms in MXDA. This renders both amine functionalities unavailable for incorporation into the network and hence does not affect polymer dynamics.
- 1 January 2022
- Published by Elsevier in SSRN Electronic Journal
Open AccessMolecular Control of Epoxy-Amine / Iron Oxide Interphase Properties
Understanding the polymer / metal oxide interphase is central to improving the performance lifetime of corrosion resistant coatings, where network polymers commonly form via step growth mechanisms in the presence of pigments. Here we employ a holistic analytical approach encompassing ATR-FTIR, DSC and molecular dynamics simulations to consider how crosslinker structure affects adsorption and incorporation into the network, using a stoichiometric mixture of DGEBA with m-xylylenediamine (MXDA) cured in the presence of hematite (Fe 2 O 3 ) and goethite (FeOOH) powders. We find that the rigid MXDA molecule has two distinct binding modes on both hematite and goethite, and that synergistic hydrogen bonding modes observed on goethite limit interconversion between the two. Moreover, we find that binding persists in fully cured composite samples, determining the levels of residual amine. In contrast to tetra-amine systems, however, we find that the Tg of composite specimens is independent of added hematite and goethite volumes. Molecular dynamics simulations demonstrate this is due to electrostatic binding between the cationic Fe sites and electronegative heteroatoms in MXDA. This renders both amine functionalities unavailable for incorporation into the network and hence does not affect polymer dynamics.
- 13 October 2021
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 37 (42) , 12409-12418
Adsorption of Epoxy Oligomers on Iron Oxide Surfaces: The Importance of Surface Treatment and the Role of Entropy
Epoxy-based coatings are widely used in a range of industries as protective coatings. The performance of the final solid–polymer system is dependent on the physicochemical properties of the interface and the interaction between the polymer and the solid substrate. In this study, we perform atomistic molecular dynamics simulations to investigate the binding of a common component in epoxy resins, diglycidyl ether of bisphenol A (DGEBA), on two iron oxide surfaces, hematite (0001) and magnetite (100), and investigate the effect of surface hydroxylation on the binding energy. We show that adsorption of DGEBA on hematite is more favorable than on magnetite and that the adsorbed molecules are highly localized on the pristine hematite surface but mobile on highly hydroxylated hematite surfaces and magnetite surfaces irregardless of surface hydroxylation fraction. A high degree of hydroxylation significantly reduces the binding energy of DGEBA on hematite but not on magnetite. The free-energy calculations confirm the trends observed upon hydroxylation, but the magnitude of the potential of mean force is lower than the binding energy due to the entropic contributions. Therefore, it can be suggested that DGEBA will adsorb more strongly on a surface containing a higher content of hematite than magnetite and that the presence of hydroxyl groups will weaken this adsorption. The presence of hydroxyl groups increases mobility of the chains, which can affect the coating rigidity.
- 13 October 2021
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 37 (42) , 12409-12418
Adsorption of Epoxy Oligomers on Iron Oxide Surfaces: The Importance of Surface Treatment and the Role of Entropy
Epoxy-based coatings are widely used in a range of industries as protective coatings. The performance of the final solid–polymer system is dependent on the physicochemical properties of the interface and the interaction between the polymer and the solid substrate. In this study, we perform atomistic molecular dynamics simulations to investigate the binding of a common component in epoxy resins, diglycidyl ether of bisphenol A (DGEBA), on two iron oxide surfaces, hematite (0001) and magnetite (100), and investigate the effect of surface hydroxylation on the binding energy. We show that adsorption of DGEBA on hematite is more favorable than on magnetite and that the adsorbed molecules are highly localized on the pristine hematite surface but mobile on highly hydroxylated hematite surfaces and magnetite surfaces irregardless of surface hydroxylation fraction. A high degree of hydroxylation significantly reduces the binding energy of DGEBA on hematite but not on magnetite. The free-energy calculations confirm the trends observed upon hydroxylation, but the magnitude of the potential of mean force is lower than the binding energy due to the entropic contributions. Therefore, it can be suggested that DGEBA will adsorb more strongly on a surface containing a higher content of hematite than magnetite and that the presence of hydroxyl groups will weaken this adsorption. The presence of hydroxyl groups increases mobility of the chains, which can affect the coating rigidity.
- Open Access
- 16 August 2021
- Published by Royal Society of Chemistry (RSC) in Soft Matter
- Vol. 17 (36) , 8343-8353
Comparison of equilibrium techniques for the viscosity calculation from DPD simulations
Dissipative Particle Dynamics (DPD) is a powerful mesoscopic modelling technique that is routinely used to predict complex fluid morphology and structural properties. While its ability to quickly scan the conformational space is well known, it is unclear if DPD can correctly calculate the viscosity of complex fluids. In this work, we estimate the viscosity of several unentangled polymer solutions using both the Einstein and Green–Kubo formulas. For this purpose, an Einstein relation is derived analogous to the revised Green–Kubo formula suggested by Jung and Schmid, J. Chem. Phys., 2016, 144, 204104. We show that the DPD simulations reproduce the dynamical behaviour predicted by the theory irrespectively of the values of the conservative and friction parameters used and estimate a Schmidt number compatible to that of a fluid system. Moreover, we observe that the Einstein method requires shorter trajectories to achieve the same statistical accuracy as the Green–Kubo formula. This work shows that DPD can confidently be used to calculate the viscosity of complex fluids and that the statistical accuracy of short trajectories can be improved by using our revised Einstein formula.
- 16 August 2021
- Published by Royal Society of Chemistry (RSC) in Soft Matter
- Vol. 17 (36) , 8343-8353
Open AccessComparison of equilibrium techniques for the viscosity calculation from DPD simulations
Dissipative Particle Dynamics (DPD) is a powerful mesoscopic modelling technique that is routinely used to predict complex fluid morphology and structural properties. While its ability to quickly scan the conformational space is well known, it is unclear if DPD can correctly calculate the viscosity of complex fluids. In this work, we estimate the viscosity of several unentangled polymer solutions using both the Einstein and Green–Kubo formulas. For this purpose, an Einstein relation is derived analogous to the revised Green–Kubo formula suggested by Jung and Schmid, J. Chem. Phys., 2016, 144, 204104. We show that the DPD simulations reproduce the dynamical behaviour predicted by the theory irrespectively of the values of the conservative and friction parameters used and estimate a Schmidt number compatible to that of a fluid system. Moreover, we observe that the Einstein method requires shorter trajectories to achieve the same statistical accuracy as the Green–Kubo formula. This work shows that DPD can confidently be used to calculate the viscosity of complex fluids and that the statistical accuracy of short trajectories can be improved by using our revised Einstein formula.
- 28 September 2020
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 36 (41) , 12288-12298
The Relationship between Wormlike Micelle Scission Free Energy and Micellar Composition: The Case of Sodium Lauryl Ether Sulfate and Cocamidopropyl Betaine
The scission energy is the difference in free energy between two hemispherical caps and the cylindrical region of a wormlike micelle. This energy difference determines the logarithm of the average micelle length, which affects several macroscopic properties such as the viscosity of viscoelastic fluids. Here we use a recently published method by Wang et al (Langmuir 2018 34 1564-1573) to directly calculate the scission energy of micelles composed of monodisperse sodium laurylethersulphate (SLESnEO), an anionic surfactant. Utilising dissipative particle dynamics (DPD), we perform a systematic study varying the number of ethoxyl groups (n) and salt concentration. The scission energy increases with increasing salt concentration, indicating that the formation of longer micelles is favoured. We attribute this to the increased charge screening that reduces the repulsion between head groups. However, the scission energy decreases with increasing number of ethoxyl groups as the flexibility of the head group increases and the sodium ion becomes less tightly bound to the head group. We then extend the analysis to look at the effect of a common co-surfactant, cocamidopropyl betaine (CAPB) and find that its addition stabilises wormlike micelles at a lower salt concentration.
- 28 September 2020
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 36 (41) , 12288-12298
The Relationship between Wormlike Micelle Scission Free Energy and Micellar Composition: The Case of Sodium Lauryl Ether Sulfate and Cocamidopropyl Betaine
The scission energy is the difference in free energy between two hemispherical caps and the cylindrical region of a wormlike micelle. This energy difference determines the logarithm of the average micelle length, which affects several macroscopic properties such as the viscosity of viscoelastic fluids. Here we use a recently published method by Wang et al (Langmuir 2018 34 1564-1573) to directly calculate the scission energy of micelles composed of monodisperse sodium laurylethersulphate (SLESnEO), an anionic surfactant. Utilising dissipative particle dynamics (DPD), we perform a systematic study varying the number of ethoxyl groups (n) and salt concentration. The scission energy increases with increasing salt concentration, indicating that the formation of longer micelles is favoured. We attribute this to the increased charge screening that reduces the repulsion between head groups. However, the scission energy decreases with increasing number of ethoxyl groups as the flexibility of the head group increases and the sodium ion becomes less tightly bound to the head group. We then extend the analysis to look at the effect of a common co-surfactant, cocamidopropyl betaine (CAPB) and find that its addition stabilises wormlike micelles at a lower salt concentration.
- 1 December 2019
- Published by Elsevier in Journal of Colloid and Interface Science
- Vol. 557, 34-44
Constructing the phase diagram of sodium laurylethoxysulfate using dissipative particle dynamics
- 1 December 2019
- Published by Elsevier in Journal of Colloid and Interface Science
- Vol. 557, 34-44
Constructing the phase diagram of sodium laurylethoxysulfate using dissipative particle dynamics
- 25 November 2019
- Published by American Physical Society (APS) in Physical Review E
- Vol. 100 (5) , 052702
Chiral nematic liquid crystals in torus-shaped and cylindrical cavities
We present a Monte Carlo simulation study of chiral nematic liquid crystals confined in torus-shaped and cylindrical cavities. For an achiral nematic with planar degenerate anchoring confined to a toroidal or cylindrical cavity, the ground state is defect free, with an untwisted director field. As chirality is introduced, the ground state remains defect free but the director field becomes twisted within the cavity. For homeotropic anchoring, the ground state for an achiral nematic within a toroidal cavity consists of two disclination rings, one large and one small, that follow the major circumference of the torus. As chirality is introduced and increased, this ground state becomes unstable with respect to twisted configurations. The closed nature of the toroidal cavity requires that only a half integer number of twists can be formed and this leads to the ground state being either a single disclination line that encircles the torus twice or a pair of intertwined disclination rings forming stable, knotted defect structures.
- 25 November 2019
- Published by American Physical Society (APS) in Physical Review E
- Vol. 100 (5) , 052702
Chiral nematic liquid crystals in torus-shaped and cylindrical cavities
We present a Monte Carlo simulation study of chiral nematic liquid crystals confined in torus-shaped and cylindrical cavities. For an achiral nematic with planar degenerate anchoring confined to a toroidal or cylindrical cavity, the ground state is defect free, with an untwisted director field. As chirality is introduced, the ground state remains defect free but the director field becomes twisted within the cavity. For homeotropic anchoring, the ground state for an achiral nematic within a toroidal cavity consists of two disclination rings, one large and one small, that follow the major circumference of the torus. As chirality is introduced and increased, this ground state becomes unstable with respect to twisted configurations. The closed nature of the toroidal cavity requires that only a half integer number of twists can be formed and this leads to the ground state being either a single disclination line that encircles the torus twice or a pair of intertwined disclination rings forming stable, knotted defect structures.
- Open Access
- 1 April 2019
- Published by Elsevier in The Journal of Chemical Thermodynamics
- Vol. 131, 620-629
Solubilities of pyrene in organic solvents: Comparison between chemical potential calculations using a cavity-based method and direct coexistence simulations
- 1 April 2019
- Published by Elsevier in The Journal of Chemical Thermodynamics
- Vol. 131, 620-629
Open AccessSolubilities of pyrene in organic solvents: Comparison between chemical potential calculations using a cavity-based method and direct coexistence simulations