Daria Ruth Galimberti
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Daria Ruth Galimberti has a total of 88 co-authors. Here are the top co-authors by number of shared publications:
- Marie-Pierre Gaigeot
Universite Paris-Saclay, Univ Evry, CY Cergy Paris Universite, CNRS, LAMBE UMR8587, 91025 Evry-Courcouronnes, France
11 shared publications - Alberto Milani
Department of Energy, Micro and Nanostructured Materials Laboratory - NanoLab, Energy, Politecnico di Milano, Via Ponzio 34/3, Milano 20133, Italy
10 shared publications - Chiara Castiglioni
Dipartimento di Chimica, Materiali e Ingegneria Chimica Giulio Natta, Politecnico di Milano, 20133 Milano, Italy
8 shared publications - Simone Pezzotti
Laboratoire CPCV, Département de Chimie, École Normale Supérieure
7 shared publications - Herma M. Cuppen
Institute for Molecules and Materials, Radboud University, Nijmegen, the Netherlands
5 shared publications - Matteo Tommasini
Politecnico di Milano
5 shared publications
Latest Publications
- Open Access
- 23 October 2026
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 22 (21) , 11471-11483
SPectroscopy Informed Configurational Enhanced Sampling (SPICES)
We present a new enhanced-sampling scheme for obtaining the underlying conformer starting from a spectrum: the SPectroscopy Informed Configurational Enhanced Sampling (SPICES) method. SPICES is driven by a history-dependent biased molecular dynamics, in which the biases are governed by spectral information rather than by energetic considerations. The performance of SPICES is hence rather insensitive to the target conformer’s relative energy compared to the global minimum structure. Furthermore, it does not require pretraining or prior knowledge of the system. We have critically assessed the performance of SPICES and demonstrated that it can identify the correct conformer of our test system, diphenylalanine (FF), among the 400 stable conformations in 22 spectral calculations (median value). We also show that SPICES is robust to noise in the experimental spectrum and to small discrepancies between theory and experiment.
- 23 October 2026
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 22 (21) , 11471-11483
Open AccessSPectroscopy Informed Configurational Enhanced Sampling (SPICES)
We present a new enhanced-sampling scheme for obtaining the underlying conformer starting from a spectrum: the SPectroscopy Informed Configurational Enhanced Sampling (SPICES) method. SPICES is driven by a history-dependent biased molecular dynamics, in which the biases are governed by spectral information rather than by energetic considerations. The performance of SPICES is hence rather insensitive to the target conformer’s relative energy compared to the global minimum structure. Furthermore, it does not require pretraining or prior knowledge of the system. We have critically assessed the performance of SPICES and demonstrated that it can identify the correct conformer of our test system, diphenylalanine (FF), among the 400 stable conformations in 22 spectral calculations (median value). We also show that SPICES is robust to noise in the experimental spectrum and to small discrepancies between theory and experiment.
- 5 July 2026
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry A
- Vol. 130 (28) , 5564-5572
The Infrared Spectra of Neutral Dimethyl-Sulfide, -Disulfide and -Sulfoxide Biomarkers in Molecular Beams
The suggested detection of sulfur-based biomarkers in the atmosphere of exoplanet K2–18b has triggered intensive discussions about the use of observational data from the James Webb Space Telescope (JWST) to search for life on exoplanets. This methodology, however, requires precise infrared spectroscopic data of biomarker molecules, ideally without the influence of pressure broadening or contaminants. In this work, we use the light of the free-electron laser FELIX to record the infrared spectrum of three sulfur-based biomarkers, dimethyl sulfide (DMS), dimethyl disulfide (DMDS), and dimethyl sulfoxide (DMSO), in the isolated conditions of a molecular beam. The use of FELIX allows the sampling of a wide spectral range, spanning from 650 to 3300 cm–1, where JWST is most sensitive. In addition, separate spectra of the (DMS)2, (DMDS)2 and (DMSO)2 dimers are recorded, allowing for a clear assessment of the effect of complexation on the infrared spectra of the molecules. Density functional theory calculations and Born–Oppenheimer molecular dynamics simulations are performed to interpret the measured infrared spectra and to estimate absolute cross sections.
- 5 July 2026
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry A
- Vol. 130 (28) , 5564-5572
The Infrared Spectra of Neutral Dimethyl-Sulfide, -Disulfide and -Sulfoxide Biomarkers in Molecular Beams
The suggested detection of sulfur-based biomarkers in the atmosphere of exoplanet K2–18b has triggered intensive discussions about the use of observational data from the James Webb Space Telescope (JWST) to search for life on exoplanets. This methodology, however, requires precise infrared spectroscopic data of biomarker molecules, ideally without the influence of pressure broadening or contaminants. In this work, we use the light of the free-electron laser FELIX to record the infrared spectrum of three sulfur-based biomarkers, dimethyl sulfide (DMS), dimethyl disulfide (DMDS), and dimethyl sulfoxide (DMSO), in the isolated conditions of a molecular beam. The use of FELIX allows the sampling of a wide spectral range, spanning from 650 to 3300 cm–1, where JWST is most sensitive. In addition, separate spectra of the (DMS)2, (DMDS)2 and (DMSO)2 dimers are recorded, allowing for a clear assessment of the effect of complexation on the infrared spectra of the molecules. Density functional theory calculations and Born–Oppenheimer molecular dynamics simulations are performed to interpret the measured infrared spectra and to estimate absolute cross sections.
- 14 May 2026
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 148 (20) , 21069-21079
Characterizing Gateway Modes for Solid–Solid Phase Transitions in Organic Crystals: The Thermosalient 4-DBpFO
We disclose the nature of the 4-DBpFO thermosalient (or "jumping") crystal phase transition 40 cm-1 (i.e., 1.2 THz) "gateway" mode using computed low-frequency Raman spectra (LFRS) benchmarked to experimental spectra and show that it is a coherent/in-phase collective motion of the inner-outer rings. We explore and reveal the roles of anharmonicity, thermal disorder, and volume fluctuations in determining the THz spectral shape. The anharmonic couplings and cell fluctuations appear to be more than just innocent spectators, and indeed, the intermolecular correlations for this mode involve multiple cells. To obtain this information, we use our Activity Weighted Velocities (AWV) method to integrate classical trajectories with the accuracy of hybrid Density Functional Theory (DFT) Raman activities. We demonstrate that this combination effectively computes accurate LFRS and offers a way to predict the Raman activity of thermally disordered systems and entropically stabilized phases, where current methods (harmonic spectra and DFT-Molecular Dynamics) fail.
- 14 May 2026
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 148 (20) , 21069-21079
Characterizing Gateway Modes for Solid–Solid Phase Transitions in Organic Crystals: The Thermosalient 4-DBpFO
We disclose the nature of the 4-DBpFO thermosalient (or "jumping") crystal phase transition 40 cm-1 (i.e., 1.2 THz) "gateway" mode using computed low-frequency Raman spectra (LFRS) benchmarked to experimental spectra and show that it is a coherent/in-phase collective motion of the inner-outer rings. We explore and reveal the roles of anharmonicity, thermal disorder, and volume fluctuations in determining the THz spectral shape. The anharmonic couplings and cell fluctuations appear to be more than just innocent spectators, and indeed, the intermolecular correlations for this mode involve multiple cells. To obtain this information, we use our Activity Weighted Velocities (AWV) method to integrate classical trajectories with the accuracy of hybrid Density Functional Theory (DFT) Raman activities. We demonstrate that this combination effectively computes accurate LFRS and offers a way to predict the Raman activity of thermally disordered systems and entropically stabilized phases, where current methods (harmonic spectra and DFT-Molecular Dynamics) fail.
- Open Access
- 7 April 2026
- Published by Scilight Press Pty Ltd in Photochemistry and Spectroscopy
- Vol. 2 (2) , 1
Bridging Molecular and Periodic Models: A Theoretical Study of the Spectroscopic Properties of Polythiophene †
Polymeric organic semiconductors play a pivotal role in the development of new organic thermoelectric materials, sensing, and organic optoelectronics. Understanding the relation between their structure and properties is essential for designing the next generation of materials with improved features. Theoretical modelling based on first-principles calculations is a powerful aid for interpreting experimental structural and spectroscopic data of π-conjugated polymers, both in their pristine and doped states. Two approaches are currently available, i.e., oligomeric models and one-dimensional (1D) periodic crystal, each of one with its advantages and disadvantages. The combination of both allows to overcome the limitations of each. However, to do this, it is fundamental to systematically address the differences between these two approaches (molecular models vs. crystals). Here, we present a first systematic comparison between oligomeric and (1D) crystal models for describing the structural and spectroscopic properties of polyconjugated polymers. Using polythiophene (PT) and poly(3-hexylthiophene) (P3HT) as prototypical benchmark systems, we examine the relationship between finite-size oligomers and the corresponding infinite polymer treated as a one-dimensional (1D) crystal under periodic boundary conditions. Density functional theory (DFT) calculations are performed to examine the convergence of geometrical parameters and vibrational features (IR and Raman) as the oligomer length increases. The 1D crystal model is then evaluated as the limiting case of the oligomeric series. Its advantages in simplifying vibrational mode analysis are illustrated. This study establishes a methodological framework connecting molecular and periodic models, with direct relevance to the interpretation of experimental vibrational spectra of conducting polymers.
- 7 April 2026
- Published by Scilight Press Pty Ltd in Photochemistry and Spectroscopy
- Vol. 2 (2) , 1
Open AccessBridging Molecular and Periodic Models: A Theoretical Study of the Spectroscopic Properties of Polythiophene †
Polymeric organic semiconductors play a pivotal role in the development of new organic thermoelectric materials, sensing, and organic optoelectronics. Understanding the relation between their structure and properties is essential for designing the next generation of materials with improved features. Theoretical modelling based on first-principles calculations is a powerful aid for interpreting experimental structural and spectroscopic data of π-conjugated polymers, both in their pristine and doped states. Two approaches are currently available, i.e., oligomeric models and one-dimensional (1D) periodic crystal, each of one with its advantages and disadvantages. The combination of both allows to overcome the limitations of each. However, to do this, it is fundamental to systematically address the differences between these two approaches (molecular models vs. crystals). Here, we present a first systematic comparison between oligomeric and (1D) crystal models for describing the structural and spectroscopic properties of polyconjugated polymers. Using polythiophene (PT) and poly(3-hexylthiophene) (P3HT) as prototypical benchmark systems, we examine the relationship between finite-size oligomers and the corresponding infinite polymer treated as a one-dimensional (1D) crystal under periodic boundary conditions. Density functional theory (DFT) calculations are performed to examine the convergence of geometrical parameters and vibrational features (IR and Raman) as the oligomer length increases. The 1D crystal model is then evaluated as the limiting case of the oligomeric series. Its advantages in simplifying vibrational mode analysis are illustrated. This study establishes a methodological framework connecting molecular and periodic models, with direct relevance to the interpretation of experimental vibrational spectra of conducting polymers.
- 25 February 2026
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 22 (5) , 2490-2500
Anharmonic Adsorption Free Energies Beyond GGA for Monomer Ethanol in H-ZSM-5
Zeolites play a pivotal role as catalysts in biomass conversion, and understanding how alcohols interact with these materials is essential for designing next-generation catalysts that reduce dependence on nonrenewable fuels and enable a more sustainable energy landscape. Achieving this understanding requires accurate free-energy calculations. In this study, we present the anharmonic enthalpy, entropy, and vibrational free energy associated with ethanol adsorption on H-ZSM-5, computed at a hybrid DFT level of theory. We critically assess the need to go beyond both GGA functionals and the harmonic approximation to achieve chemically accurate vibrational free energies. The vibrational contributions to the free energy are obtained using a newly implemented strategy that combines the DOS-P method with a QM/QM molecular-dynamics scheme. This approach yields significant improvement in accuracy compared to the commonly used PBE + D2 functional. By combining the anharmonic vibrational free energies computed at the B3LYP + D2 level with electronic energies refined at the CCSD(T) level, we obtain an anharmonic adsorption enthalpy, ΔH, of -101.9 kJ mol-1, the entropic term, -TΔS, equal to 61.0 kJ mol-1, and a free energy of adsorption, ΔG, of -40.8 kJ mol-1 at 313 K. We performed a critical comparison between these values and the experimental data. Notably, the entropic term agrees well with the experimental value of 62.0 kJ mol-1. We also examine the significance of anharmonicity, finding a correction of approximately 1 kJ mol-1 for both the vibrational entropy and enthalpy, and about 3.7 kJ mol-1 for the ZPE.
- 25 February 2026
- Published by American Chemical Society (ACS) in Journal of Chemical Theory and Computation
- Vol. 22 (5) , 2490-2500
Anharmonic Adsorption Free Energies Beyond GGA for Monomer Ethanol in H-ZSM-5
Zeolites play a pivotal role as catalysts in biomass conversion, and understanding how alcohols interact with these materials is essential for designing next-generation catalysts that reduce dependence on nonrenewable fuels and enable a more sustainable energy landscape. Achieving this understanding requires accurate free-energy calculations. In this study, we present the anharmonic enthalpy, entropy, and vibrational free energy associated with ethanol adsorption on H-ZSM-5, computed at a hybrid DFT level of theory. We critically assess the need to go beyond both GGA functionals and the harmonic approximation to achieve chemically accurate vibrational free energies. The vibrational contributions to the free energy are obtained using a newly implemented strategy that combines the DOS-P method with a QM/QM molecular-dynamics scheme. This approach yields significant improvement in accuracy compared to the commonly used PBE + D2 functional. By combining the anharmonic vibrational free energies computed at the B3LYP + D2 level with electronic energies refined at the CCSD(T) level, we obtain an anharmonic adsorption enthalpy, ΔH, of -101.9 kJ mol-1, the entropic term, -TΔS, equal to 61.0 kJ mol-1, and a free energy of adsorption, ΔG, of -40.8 kJ mol-1 at 313 K. We performed a critical comparison between these values and the experimental data. Notably, the entropic term agrees well with the experimental value of 62.0 kJ mol-1. We also examine the significance of anharmonicity, finding a correction of approximately 1 kJ mol-1 for both the vibrational entropy and enthalpy, and about 3.7 kJ mol-1 for the ZPE.
- Open Access
- 28 October 2025
- Published by Elsevier in Vibrational Spectroscopy
- Vol. 141
Influence of the computational methods in the conformational assignment of experimental infrared spectra
- 28 October 2025
- Published by Elsevier in Vibrational Spectroscopy
- Vol. 141
Open AccessInfluence of the computational methods in the conformational assignment of experimental infrared spectra
- Open Access
- 23 May 2025
- Published by American Chemical Society (ACS) in ACS Earth and Space Chemistry
- Vol. 9 (6) , 1580-1592
Infrared Irradiation of H2O:CO2 Ice: A Combined Experimental and Computational Study of the Dissipation of CO2 Vibrational Excitations
In interstellar ices, the ice matrix can have a great influence on the chemical reactions. The hydrogen-bonding network in pure water ices facilitates fast energy dissipation that, for example, stabilizes the HOCO complex, a crucial step in the formation of CO2. To better understand the energy dynamics and its possible influence on the processes in the ice, we investigated a H2O:CO2 1:4 ice mixture exposed to infrared irradiation on-resonance with the CO2 vibrations. Experimentally, we find changes in the OH stretch of H2O after irradiating the asymmetric stretch of CO2 for several minutes with the intense monochromatic light of the FELIX free electron lasers. Using molecular dynamics simulations, we found that an excitation of the asymmetric stretch of CO2 readily dissipates to other asymmetric stretches in the environment, but only dissipates to the CO2 libration and H2O twist modes after roughly 2 ns because of its minimal anharmonicity and coupling with other modes. This is significantly longer than the off-time between laser pulses of 1 ns, suggesting ladder climbing or that the stacking of the excitation boosts the experimentally observed changes. For infrared excitation of the CO2 bending vibration, the simulations reveal a fast distribution of energy and coupling to the intermolecular interactions that lead to thermal heating of the H2O vibrational modes. This is not observed on the time scale of the experiments. Still, both simulations and experiments reveal nonthermal annealing of the H2O component of the mixed ice when exposed to infrared irradiation on-resonance with the CO2 vibrations.
- 23 May 2025
- Published by American Chemical Society (ACS) in ACS Earth and Space Chemistry
- Vol. 9 (6) , 1580-1592
Open AccessInfrared Irradiation of H2O:CO2 Ice: A Combined Experimental and Computational Study of the Dissipation of CO2 Vibrational Excitations
In interstellar ices, the ice matrix can have a great influence on the chemical reactions. The hydrogen-bonding network in pure water ices facilitates fast energy dissipation that, for example, stabilizes the HOCO complex, a crucial step in the formation of CO2. To better understand the energy dynamics and its possible influence on the processes in the ice, we investigated a H2O:CO2 1:4 ice mixture exposed to infrared irradiation on-resonance with the CO2 vibrations. Experimentally, we find changes in the OH stretch of H2O after irradiating the asymmetric stretch of CO2 for several minutes with the intense monochromatic light of the FELIX free electron lasers. Using molecular dynamics simulations, we found that an excitation of the asymmetric stretch of CO2 readily dissipates to other asymmetric stretches in the environment, but only dissipates to the CO2 libration and H2O twist modes after roughly 2 ns because of its minimal anharmonicity and coupling with other modes. This is significantly longer than the off-time between laser pulses of 1 ns, suggesting ladder climbing or that the stacking of the excitation boosts the experimentally observed changes. For infrared excitation of the CO2 bending vibration, the simulations reveal a fast distribution of energy and coupling to the intermolecular interactions that lead to thermal heating of the H2O vibrational modes. This is not observed on the time scale of the experiments. Still, both simulations and experiments reveal nonthermal annealing of the H2O component of the mixed ice when exposed to infrared irradiation on-resonance with the CO2 vibrations.
- Open Access
- 15 May 2025
- Published by Royal Society of Chemistry (RSC) in Chemical Communications
- Vol. 61 (49) , 8903-8906
Poly( N -acyliminophosphoranes): main chain functionalized poly(ylides) with pH-dependent degradation
Zwitterionic polymers are often non-degradable due to their polyolefin backbones. We present an Fe( ii )-catalyzed step-growth polymerization that incorporates N -acyliminophosphoranes into the main chain, enabling tailored pH-responsive degradability.
- 15 May 2025
- Published by Royal Society of Chemistry (RSC) in Chemical Communications
- Vol. 61 (49) , 8903-8906
Open AccessPoly( N -acyliminophosphoranes): main chain functionalized poly(ylides) with pH-dependent degradation
Zwitterionic polymers are often non-degradable due to their polyolefin backbones. We present an Fe( ii )-catalyzed step-growth polymerization that incorporates N -acyliminophosphoranes into the main chain, enabling tailored pH-responsive degradability.
- Open Access
- 24 March 2025
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 41 (13) , 8627-8636
Water-Soluble Sulfur-Ylide-Functionalized Polyacrylamides for Antibacterial Surface Applications
Surface fouling induced by biomolecules and microorganisms remains a persistent challenge in materials science, particularly in healthcare applications, where biofilm formation on medical devices may lead to infections and antimicrobial resistance. Antifouling strategies typically rely on the formation of either hydration layers or cytotoxic materials for direct antimicrobial effects. Recent advances in zwitterionic polymers derived from ylides offer a promising yet unexplored toolbox for the construction of antifouling and antimicrobial coatings. While N-oxide-based ylides have been extensively studied as building blocks for antifouling materials, sulfur-ylide-based materials, and the precise underlying mechanisms remain underexplored despite their broader chemical versatility. Here, we present a fully water-soluble acrylamide-based poly(sulfur ylide) and compare its properties to those of previously reported hydrophobic polystyrene-derived analogues. Notably, water-soluble poly(sulfur ylides) retain antimicrobial efficacy on surfaces but lose cytotoxicity in solution, unlike its hydrophobic counterpart. Computational studies reveal that the dipole moment of sulfur ylides is environmentally responsive, stabilizing in hydrophobic environments. Genetic analysis confirms outer membrane destabilization for both polymers but suggests that the hydrophobicity of the polystyrene backbone promotes stronger interactions. We suggest that future work should focus on elucidating additional interactions, including supramolecular behaviors of amphiphilic sulfur ylides, to refine their structure-property relationships and optimize their antifouling and antimicrobial properties.
- 24 March 2025
- Published by American Chemical Society (ACS) in Langmuir
- Vol. 41 (13) , 8627-8636
Open AccessWater-Soluble Sulfur-Ylide-Functionalized Polyacrylamides for Antibacterial Surface Applications
Surface fouling induced by biomolecules and microorganisms remains a persistent challenge in materials science, particularly in healthcare applications, where biofilm formation on medical devices may lead to infections and antimicrobial resistance. Antifouling strategies typically rely on the formation of either hydration layers or cytotoxic materials for direct antimicrobial effects. Recent advances in zwitterionic polymers derived from ylides offer a promising yet unexplored toolbox for the construction of antifouling and antimicrobial coatings. While N-oxide-based ylides have been extensively studied as building blocks for antifouling materials, sulfur-ylide-based materials, and the precise underlying mechanisms remain underexplored despite their broader chemical versatility. Here, we present a fully water-soluble acrylamide-based poly(sulfur ylide) and compare its properties to those of previously reported hydrophobic polystyrene-derived analogues. Notably, water-soluble poly(sulfur ylides) retain antimicrobial efficacy on surfaces but lose cytotoxicity in solution, unlike its hydrophobic counterpart. Computational studies reveal that the dipole moment of sulfur ylides is environmentally responsive, stabilizing in hydrophobic environments. Genetic analysis confirms outer membrane destabilization for both polymers but suggests that the hydrophobicity of the polystyrene backbone promotes stronger interactions. We suggest that future work should focus on elucidating additional interactions, including supramolecular behaviors of amphiphilic sulfur ylides, to refine their structure-property relationships and optimize their antifouling and antimicrobial properties.
- Open Access
- 10 December 2024
- Published by Royal Society of Chemistry (RSC) in Physical Chemistry Chemical Physics
- Vol. 27 (1) , 550-563
Assignment of IR spectra of ethanol at Brønsted sites of H-ZSM-5 to monomer adsorption using a Fermi resonance model
Understanding how alcohol molecules interact with the Brønsted acid sites (BAS) of zeolites is a prerequisite to the design of zeolite catalysts and catalytic processes. Here, we report IR spectra for the adsorption of ethanol on a highly crystalline sample of H-ZSM-5 zeolites exposed to ethanol gas at increasing pressure. We use density functional theory in combination with a FERMI resonance model to assign the measured spectra to a single adsorbed ethanol molecule per BAS. Specifically, we assign the bands at 2450 cm−1 and 1670 cm−1 to a FERMI resonance between the fundamental (Z)O–H stretching band of a single-ethanol-loaded BAS and the first overtone of the (Z)O–H out-of-plane bending. We conclude that adsorbed dimers do not contribute in a noticeable way up to a concentration of almost one ethanol molecule per BAS site. We further show that hybrid functionals (B3LYP) are required to get a close match between the predicted and experimental spectra, whereas commonly used generalized gradient type functionals such as PBE incorrectly describe the potential energy surface. They overestimate the redshift of the OH stretching band on hydrogen bond formation which results in an erroneous assignment of the IR bands.
- 10 December 2024
- Published by Royal Society of Chemistry (RSC) in Physical Chemistry Chemical Physics
- Vol. 27 (1) , 550-563
Open AccessAssignment of IR spectra of ethanol at Brønsted sites of H-ZSM-5 to monomer adsorption using a Fermi resonance model
Understanding how alcohol molecules interact with the Brønsted acid sites (BAS) of zeolites is a prerequisite to the design of zeolite catalysts and catalytic processes. Here, we report IR spectra for the adsorption of ethanol on a highly crystalline sample of H-ZSM-5 zeolites exposed to ethanol gas at increasing pressure. We use density functional theory in combination with a FERMI resonance model to assign the measured spectra to a single adsorbed ethanol molecule per BAS. Specifically, we assign the bands at 2450 cm−1 and 1670 cm−1 to a FERMI resonance between the fundamental (Z)O–H stretching band of a single-ethanol-loaded BAS and the first overtone of the (Z)O–H out-of-plane bending. We conclude that adsorbed dimers do not contribute in a noticeable way up to a concentration of almost one ethanol molecule per BAS site. We further show that hybrid functionals (B3LYP) are required to get a close match between the predicted and experimental spectra, whereas commonly used generalized gradient type functionals such as PBE incorrectly describe the potential energy surface. They overestimate the redshift of the OH stretching band on hydrogen bond formation which results in an erroneous assignment of the IR bands.