🔬 condensed matter

Pressure effects in the properties of simple monohydric alcohols. Lessons from molecular dynamics simulations of united atom type UAM-EW model

This study utilizes isobaric-isothermal molecular dynamics simulations with a united atom force field to investigate how pressure (up to 3 kbar) affects the macroscopic properties and microscopic structure of methanol, ethanol, and 1-propanol, validating the model against experimental data.

Original authors: M. Aguilar, L. Pusztai, O. Pizio

Published 2026-03-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Aguilar, L. Pusztai, O. Pizio

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Squeezing the Squeeze

Imagine you have a crowded dance floor filled with people holding hands in specific patterns. This is like a liquid alcohol (like the stuff in your hand sanitizer or beer). Usually, these people (molecules) are dancing at a comfortable pace, holding hands loosely, and moving around freely.

This paper asks a simple question: What happens to this dance floor if we shrink the room and push everyone closer together?

The scientists used a super-powerful computer to simulate what happens when they squeeze three common alcohols—Methanol, Ethanol, and 1-Propanol—from normal room pressure all the way up to 3,000 times that pressure (3 kbar). That's like going from the surface of the Earth to the bottom of the deepest ocean trench, but for liquids.

The Tools: The "United Atom" Lego Set

To run this simulation, the researchers didn't build a model of every single tiny electron in the alcohol molecules. That would take too long. Instead, they used a "United Atom" model.

The Analogy: Imagine building a house with Lego bricks. A "United Atom" model treats a whole group of atoms (like a carbon chain) as a single, big Lego brick. It's a simplified version, but it's fast and surprisingly accurate.
The Goal: They wanted to see if this simplified "Lego" model could predict how the alcohol behaves under extreme pressure, and if those predictions matched real-world experiments.

What They Found: The Dance Floor Changes

Here are the four main things they looked at, explained simply:

1. Density: The Crowd Gets Tighter

What it is: How many people fit in the room.
The Result: As they squeezed the room (increased pressure), the alcohol got denser. The molecules packed closer together.
The Verdict: The computer model was very good at predicting this. For Methanol, it was spot-on. For Ethanol and Propanol, it was a little off (predicting they got a bit too dense too quickly), but overall, the model worked well.

2. Compressibility: How Easy is it to Squeeze?

What it is: How much the liquid shrinks when you push on it. Think of a sponge (easy to squeeze) vs. a rock (hard to squeeze).
The Result: As the pressure went up, the alcohol became harder to squeeze (less compressible). The molecules were already so close that there was no more room to push them together.
The Verdict: The model got the shape of the curve right, even if the exact numbers were slightly different from real life.

3. The Dielectric Constant: The "Stickiness" of Electricity

What it is: This sounds fancy, but it's basically a measure of how well the liquid can hold onto an electric charge. Alcohol molecules are like tiny magnets (dipoles).
The Result: When you squeeze the room, the "magnets" are forced closer together. They start to align better and hold onto each other more strongly. So, the liquid becomes better at conducting electricity (the dielectric constant goes up).
The Verdict: The model predicted this increase perfectly. The more you squeeze, the more "electrically sticky" the alcohol gets.

4. Self-Diffusion: The Dance Moves Slow Down

What it is: How fast the molecules can move around and swap places.
The Result: In a crowded room, you can't run. As the pressure increased, the molecules got stuck. They moved much slower.
The Verdict: The model correctly predicted that the alcohol would get "sluggish" under pressure. Methanol slowed down the most, while the bigger alcohol molecules (Ethanol and Propanol) slowed down a bit less.

The Secret Sauce: The Hydrogen Bond Network

The most interesting part of the paper is about the microscopic structure.

The Analogy: Imagine the alcohol molecules are people holding hands (Hydrogen bonds).
The Discovery: When the room shrank, the people didn't let go of each other's hands. In fact, the "hand-holding" (hydrogen bonding) stayed almost exactly the same!
The Twist: While the hands stayed the same, the tails of the molecules (the non-polar parts) got squished and rearranged. The molecules had to twist and turn to fit into the smaller space, creating a more ordered, rigid structure, even though the core "hand-holding" didn't change much.

Why Does This Matter?

You might ask, "Who cares about squeezing alcohol?"

Better Models: This study proves that their simplified computer model (the Lego set) is a reliable tool. Scientists can now use it to predict how alcohols behave in places we can't easily test, like inside high-pressure industrial machines or deep underground.
Future Experiments: The computer found things that real scientists haven't measured yet (like the exact structure of Ethanol under high pressure). This gives experimentalists a "map" of what to look for in their labs.
Understanding Mixtures: Since alcohol and water mix so well, understanding how alcohol behaves under pressure helps us understand how water-alcohol mixtures behave in extreme environments, which is crucial for chemistry and biology.

The Bottom Line

The researchers successfully used a computer simulation to "squeeze" three types of alcohol. They found that while the molecules get packed tighter and move slower, their core connections (hydrogen bonds) remain surprisingly stable. The computer model worked so well that it can now be trusted to guide future experiments and help us understand the hidden world of liquids under extreme pressure.

1. Problem Statement

The study addresses the lack of comprehensive understanding regarding the pressure dependence of thermodynamic, dynamic, and dielectric properties of simple monohydric alcohols (methanol, ethanol, and 1-propanol) at the molecular level. While previous research has focused on water-alcohol mixtures or ambient pressure conditions, there is a scarcity of data and theoretical models describing the behavior of neat alcohols under high pressure (up to 3 kbar). Furthermore, existing force fields often fail to simultaneously reproduce multiple target properties (density, dielectric constant, surface tension) accurately across varying pressure ranges. The authors aim to validate a recently proposed united-atom force field (UAM-EW) against experimental data and elucidate the microscopic structural evolution of these alcohols under compression.

2. Methodology

System and Models: The study utilizes the UAM-EW united-atom, non-polarizable force field for methanol (MeOH), ethanol (EtOH), and 1-propanol (PrOH). In this model, methyl ( $CH_3$ ) and methylene ( $CH_2$ ) groups are treated as single interaction sites, while oxygen (O) and hydrogen (H) are explicit. The potential combines Lennard-Jones (LJ) and Coulombic terms.
Simulation Protocol:
- Ensemble: Isothermal-isobaric (NPT) ensemble.
- Conditions: Temperature fixed at 298.15 K; Pressure varied from 1 bar to 3000 bar (3 kbar).
- Software: GROMACS 5.1.2.
- System Size: 3000 molecules in a cubic box with periodic boundary conditions.
- Parameters: Time step of 0.002 ps; V-rescale thermostat ( $\tau_T = 0.5$ ps) and Parrinello-Rahman barostat ( $\tau_P = 2.0$ ps).
- Interactions: Non-bonded cutoff at 1.1 nm; Particle Mesh Ewald (PME) for long-range electrostatics; LINCS algorithm for bond constraints.
Analysis:
- Thermodynamics: Density ( $\rho$ ) and isothermal compressibility ( $\kappa_T$ ).
- Dielectric: Static dielectric constant ( $\varepsilon$ ) calculated via fluctuations of the total dipole moment.
- Dynamics: Self-diffusion coefficients ( $D$ ) derived from the mean squared displacement (Einstein relation).
- Microstructure: Site-site pair distribution functions (PDFs, $g(r)$ ), coordination numbers, hydrogen bond counts ( $n_{HB}$ ), and the distance-dependent Kirkwood factor ( $G_K(r)$ ).

3. Key Contributions

Validation of UAM-EW Model: The paper provides a rigorous validation of the UAM-EW force field under high-pressure conditions, demonstrating its capability to predict properties for three different alcohols simultaneously.
High-Pressure Structural Insights: It offers novel insights into the microscopic structural changes of neat alcohols under pressure, a regime where experimental data (particularly for EtOH and PrOH) is currently unavailable.
Decoupling of Polar and Non-Polar Effects: The study distinguishes between the behavior of the polar hydroxyl groups (which remain relatively stable) and the non-polar alkyl chains (which undergo significant ordering) under compression.
Kirkwood Factor Analysis: The authors introduce the analysis of the finite-system Kirkwood factor ( $G_K$ ) to understand dipole alignment and cooperative hydrogen bonding effects under pressure.

4. Key Results

A. Thermodynamic and Dynamic Properties

Density: The UAM-EW model shows excellent agreement with experimental density data for methanol. For ethanol, the model slightly overestimates density growth at high pressures, while for 1-propanol, it underestimates density at low pressures but improves at intermediate pressures (0.8–1.4 kbar). Overall, the model is deemed satisfactory across the 1–3 kbar range.
Isothermal Compressibility ( $\kappa_T$ ): While absolute values show some deviation from experiments, the trend of $\kappa_T$ decreasing with pressure is accurately reproduced by the simulations for all three alcohols.
Dielectric Constant ( $\varepsilon$ ): $\varepsilon$ increases with pressure for all three alcohols due to reduced intermolecular distances and stronger dipole correlations. The model captures this trend well, with deviations of only a few percent. Methanol shows the most pronounced increase.
Self-Diffusion Coefficient ( $D$ ): $D$ decreases with increasing pressure for all alcohols, consistent with increased density. The model underestimates $D$ for methanol across the entire range but shows better agreement for ethanol and 1-propanol.

B. Microscopic Structure

Hydrogen Bonding Network:
- The average number of hydrogen bonds per molecule ( $n_{HB}$ ) remains nearly constant (e.g., MeOH: 1.923 at 1 bar $\to$ 1.935 at 3 kbar).
- The first coordination shell of oxygen atoms ( $g_{O-O}$ and $g_{O-H}$ ) is largely preserved, indicating that the local H-bond network is robust against pressure.
- Pressure primarily affects the cooperative structure of the network rather than the number of bonds.
Non-Polar Tail Ordering:
- Significant structural changes occur in the non-polar regions. The pair distribution functions for carbon sites ( $g_{C-C}$ ) show increased oscillation magnitude and the emergence of distinct coordination shells (up to four for MeOH) as pressure rises.
- This indicates a transition toward a more ordered arrangement of the alkyl "tails" under compression.
Kirkwood Factor ( $G_K$ ):
- $G_K(r)$ exhibits oscillations corresponding to coordination shells, indicating preferential parallel alignment of dipole moments.
- The effect of pressure on $G_K$ is most significant for methanol. The increase in $G_K$ at high pressure is attributed to the evolution of cooperative hydrogen bonding features rather than changes in the pair-level bond count.

5. Significance and Future Directions

Scientific Impact: This work bridges the gap between theoretical modeling and experimental data for alcohols at high pressures. It confirms that united-atom models can effectively capture complex pressure-induced structural transitions, specifically the ordering of non-polar chains while maintaining the integrity of the hydrogen bond network.
Guidance for Experiments: The simulation results provide a predictive framework for future high-pressure experiments, particularly for ethanol and 1-propanol where experimental data is currently lacking.
Model Development: The authors suggest that while the UAM-EW model is robust, future work should focus on:
- Calculating adiabatic compressibility and speed of sound.
- Simulating viscosity and dielectric relaxation phenomena.
- Validating partial pair distribution functions against X-ray and neutron diffraction data (currently unavailable for high-pressure EtOH/PrOH).
- Refining force fields to better capture the rigidity of alkyl chains under extreme compression.

In conclusion, the paper successfully demonstrates that the UAM-EW united-atom model is a reliable tool for investigating the pressure-dependent behavior of simple monohydric alcohols, revealing that pressure primarily induces ordering in the non-polar domains while preserving the local hydrogen-bonding topology.