Gyorgy Hantal | Universtät Wien (original) (raw)
Papers by Gyorgy Hantal
The Journal of Physical Chemistry C, 2010
Substantial improvements in the molecular level understanding of fluid interfaces have recently b... more Substantial improvements in the molecular level understanding of fluid interfaces have recently been achieved by recognizing the importance of detecting the intrinsic surface of the coexisting condensed phases in computer simulations (i.e., after the removal of corrugations caused by capillary waves) and by developing several methods for identifying the molecules that are indeed located at the boundary of the two phases. In our previous paper [J. Phys. Chem. C 2010, 114, 11169], we critically compared those methods in terms of reliability, robustness, and computation speed. Once the intrinsic surface of a given phase is detected, various profiles, such as the density profiles of the components, can be calculated relative to this intrinsic surface rather than to the macroscopically planar Gibbs dividing surface. As a continuation of our previous study, here we present a detailed and critical comparison of various methods that can be used to calculate intrinsic density profiles once the full set of truly interfacial molecules has been identified. Two of the methods, the Fourier function and the Voronoi tessellation, are already described in the literature; two other methods, the covering surface and the triangular interpolation, are newly proposed algorithms; one method, the modified grid-based intrinsic profile (GIP) method, is an improvement over an existing procedure. The different methods are again compared in terms of accuracy and computational cost. On the basis of this comparison, we propose a fast and accurate protocol to be routinely used for intrinsic surface analyses in computer simulations.
The Journal of Physical Chemistry C, 2009
ABSTRACT
Langmuir, 2010
Adsorption study of benzaldehyde on ice surfaces is performed by combining experimental and theor... more Adsorption study of benzaldehyde on ice surfaces is performed by combining experimental and theoretical approaches. The experiments are conducted over the temperature range 233-253 K using a coated wall flow tube coupled to a mass spectrometric detector. Besides the experimental way, the adsorption isotherm is also determined by performing a set of grand canonical Monte Carlo simulations at 233 K. The experimental and calculated adsorption isotherms show a very good agreement within the corresponding errors. Besides, both experimental and theoretical studies permit us to derive the enthalpy of adsorption of benzaldehyde on ice surfaces ΔH ads , which are in excellent agreement: ΔH ads =-61.4 ( 9.7 kJ/mol (experimental) and ΔH ads =-59.4 ( 5.1 kJ/mol (simulation). The obtained results indicate a much stronger ability of benzaldehyde of being adsorbed at the surface of ice than that of small aliphatic aldehydes, such as formaldehyde or acetaldehyde. At low surface coverages the adsorbed molecules exclusively lie parallel with the ice surface. With increasing surface coverage, however, the increasing competition of the adsorbed molecules for the surface area to be occupied leads to the appearance of two different perpendicular orientations relative to the surface. In the first orientation, the benzaldehyde molecule turns its aldehyde group toward the ice phase, and, similarly to the molecules in the lying orientation, forms a hydrogen bond with a surface water molecule. In the other perpendicular orientation the aldehyde group turns to the vapor phase, and its O atom interacts with the delocalized π system of the benzene ring of a nearby lying benzaldehyde molecule of the second molecular layer. In accordance with this observed scenario, the saturated adsorption layer, being stable in a roughly 1 kJ/mol broad range of chemical potentials, contains, besides the first molecular layer, also traces of the second molecular layer of adsorbed benzaldehyde. *Corresponding authors. B€ urger, V.; Droste-Fanke, B.; Grimm, F.; Krieger, A.; Schneider, J.; Stilp, T. Geophys. Res. Lett. 1997, 24, 3017. (5) Mari, C.; Sa€ ut, C.; Jacob, D.; Ravetta, F.; Anderson, B.; Avery, M. A.; Blake, D. R.; Brune, W. H.; Faloona, I.; Gregory, G. L.; Heikes, B. G.; Sachse, G. W.; Sandholm, S. T.; Singh, H. B.; Talbot, R. W.; Tan, D.; Vay, S. J. Geophys. Res. 2003, 108-D2, 8229. (6) Singh, H.; Chen, Y.; Staudt, A.; Jacob, D.; Blake, D.; Heikes, B.; Snow, J. Nature 2001, 410, 1078. (7) Singh, H.; Chen, Y.; Tabazadeh, A.; Fukui, Y.; Bey, I.; Yantosca, R.; Jacob, D.; Arnold, F.; Wohlfrom, K.; Atlas, E.; Flocke, F.; Blake, D.; Heikes, B.; Snow, J.; Talbot, R.; Gregory, G.; Sachse, G.; Vay, S.; Kondo, Y.
The Journal of Chemical Physics, 2010
The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water mole... more The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water molecules on different types of model soot particles. These soot models are constructed by first removing atoms from onion-fullerene structures in order to create randomly distributed pores inside the soot, and then performing molecular dynamics simulations, based on the reactive adaptive intermolecular reactive empirical bond order ͑AIREBO͒ description of the interaction between carbon atoms, to optimize the resulting structures. The obtained results clearly show that the main driving force of water adsorption on soot is the possibility of the formation of new water-water hydrogen bonds with the already adsorbed water molecules. The shape of the calculated water adsorption isotherms at 298 K strongly depends on the possible confinement of the water molecules in pores of the carbonaceous structure. We found that there are two important factors influencing the adsorption ability of soot. The first of these factors, dominating at low pressures, is the ability of the soot of accommodating the first adsorbed water molecules at strongly hydrophilic sites. The second factor concerns the size and shape of the pores, which should be such that the hydrogen bonding network of the water molecules filling them should be optimal. This second factor determines the adsorption properties at higher pressures.
The Journal of Chemical Physics, 2009
A general formalism for introducing nuclear quantum effects in the expression of the quantum time... more A general formalism for introducing nuclear quantum effects in the expression of the quantum time correlation function of an operator in a multilevel electronic system is presented in the adiabatic limit. The final formula includes the nuclear quantum time correlation functions of the operator matrix elements, of the energy gap, and their cross terms. These quantities can be inferred and evaluated from their classical analogs obtained by mixed quantum-classical molecular dynamics simulations. The formalism is applied to the absorption spectrum of a hydrated electron, expressed in terms of the time correlation function of the dipole operator in the ground electronic state. We find that both static and dynamic nuclear quantum effects distinctly influence the shape of the absorption spectrum, especially its high energy tail related to transitions to delocalized electron states. Their inclusion does improve significantly the agreement between theory and experiment for both the low and high frequency edges of the spectrum. It does not appear sufficient, however, to resolve persistent deviations in the slow Lorentzian-like decay part of the spectrum in the intermediate 2-3 eV region.
The Journal of Physical Chemistry C, 2010
Substantial improvements in the molecular level understanding of fluid interfaces have recently b... more Substantial improvements in the molecular level understanding of fluid interfaces have recently been achieved by recognizing the importance of detecting the intrinsic surface of the coexisting condensed phases in computer simulations (i.e., after the removal of corrugations caused by capillary waves) and by developing several methods for identifying the molecules that are indeed located at the boundary of the two phases. In our previous paper [J. Phys. Chem. C 2010, 114, 11169], we critically compared those methods in terms of reliability, robustness, and computation speed. Once the intrinsic surface of a given phase is detected, various profiles, such as the density profiles of the components, can be calculated relative to this intrinsic surface rather than to the macroscopically planar Gibbs dividing surface. As a continuation of our previous study, here we present a detailed and critical comparison of various methods that can be used to calculate intrinsic density profiles once the full set of truly interfacial molecules has been identified. Two of the methods, the Fourier function and the Voronoi tessellation, are already described in the literature; two other methods, the covering surface and the triangular interpolation, are newly proposed algorithms; one method, the modified grid-based intrinsic profile (GIP) method, is an improvement over an existing procedure. The different methods are again compared in terms of accuracy and computational cost. On the basis of this comparison, we propose a fast and accurate protocol to be routinely used for intrinsic surface analyses in computer simulations.
The Journal of Physical Chemistry C, 2009
ABSTRACT
Langmuir, 2010
Adsorption study of benzaldehyde on ice surfaces is performed by combining experimental and theor... more Adsorption study of benzaldehyde on ice surfaces is performed by combining experimental and theoretical approaches. The experiments are conducted over the temperature range 233-253 K using a coated wall flow tube coupled to a mass spectrometric detector. Besides the experimental way, the adsorption isotherm is also determined by performing a set of grand canonical Monte Carlo simulations at 233 K. The experimental and calculated adsorption isotherms show a very good agreement within the corresponding errors. Besides, both experimental and theoretical studies permit us to derive the enthalpy of adsorption of benzaldehyde on ice surfaces ΔH ads , which are in excellent agreement: ΔH ads =-61.4 ( 9.7 kJ/mol (experimental) and ΔH ads =-59.4 ( 5.1 kJ/mol (simulation). The obtained results indicate a much stronger ability of benzaldehyde of being adsorbed at the surface of ice than that of small aliphatic aldehydes, such as formaldehyde or acetaldehyde. At low surface coverages the adsorbed molecules exclusively lie parallel with the ice surface. With increasing surface coverage, however, the increasing competition of the adsorbed molecules for the surface area to be occupied leads to the appearance of two different perpendicular orientations relative to the surface. In the first orientation, the benzaldehyde molecule turns its aldehyde group toward the ice phase, and, similarly to the molecules in the lying orientation, forms a hydrogen bond with a surface water molecule. In the other perpendicular orientation the aldehyde group turns to the vapor phase, and its O atom interacts with the delocalized π system of the benzene ring of a nearby lying benzaldehyde molecule of the second molecular layer. In accordance with this observed scenario, the saturated adsorption layer, being stable in a roughly 1 kJ/mol broad range of chemical potentials, contains, besides the first molecular layer, also traces of the second molecular layer of adsorbed benzaldehyde. *Corresponding authors. B€ urger, V.; Droste-Fanke, B.; Grimm, F.; Krieger, A.; Schneider, J.; Stilp, T. Geophys. Res. Lett. 1997, 24, 3017. (5) Mari, C.; Sa€ ut, C.; Jacob, D.; Ravetta, F.; Anderson, B.; Avery, M. A.; Blake, D. R.; Brune, W. H.; Faloona, I.; Gregory, G. L.; Heikes, B. G.; Sachse, G. W.; Sandholm, S. T.; Singh, H. B.; Talbot, R. W.; Tan, D.; Vay, S. J. Geophys. Res. 2003, 108-D2, 8229. (6) Singh, H.; Chen, Y.; Staudt, A.; Jacob, D.; Blake, D.; Heikes, B.; Snow, J. Nature 2001, 410, 1078. (7) Singh, H.; Chen, Y.; Tabazadeh, A.; Fukui, Y.; Bey, I.; Yantosca, R.; Jacob, D.; Arnold, F.; Wohlfrom, K.; Atlas, E.; Flocke, F.; Blake, D.; Heikes, B.; Snow, J.; Talbot, R.; Gregory, G.; Sachse, G.; Vay, S.; Kondo, Y.
The Journal of Chemical Physics, 2010
The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water mole... more The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water molecules on different types of model soot particles. These soot models are constructed by first removing atoms from onion-fullerene structures in order to create randomly distributed pores inside the soot, and then performing molecular dynamics simulations, based on the reactive adaptive intermolecular reactive empirical bond order ͑AIREBO͒ description of the interaction between carbon atoms, to optimize the resulting structures. The obtained results clearly show that the main driving force of water adsorption on soot is the possibility of the formation of new water-water hydrogen bonds with the already adsorbed water molecules. The shape of the calculated water adsorption isotherms at 298 K strongly depends on the possible confinement of the water molecules in pores of the carbonaceous structure. We found that there are two important factors influencing the adsorption ability of soot. The first of these factors, dominating at low pressures, is the ability of the soot of accommodating the first adsorbed water molecules at strongly hydrophilic sites. The second factor concerns the size and shape of the pores, which should be such that the hydrogen bonding network of the water molecules filling them should be optimal. This second factor determines the adsorption properties at higher pressures.
The Journal of Chemical Physics, 2009
A general formalism for introducing nuclear quantum effects in the expression of the quantum time... more A general formalism for introducing nuclear quantum effects in the expression of the quantum time correlation function of an operator in a multilevel electronic system is presented in the adiabatic limit. The final formula includes the nuclear quantum time correlation functions of the operator matrix elements, of the energy gap, and their cross terms. These quantities can be inferred and evaluated from their classical analogs obtained by mixed quantum-classical molecular dynamics simulations. The formalism is applied to the absorption spectrum of a hydrated electron, expressed in terms of the time correlation function of the dipole operator in the ground electronic state. We find that both static and dynamic nuclear quantum effects distinctly influence the shape of the absorption spectrum, especially its high energy tail related to transitions to delocalized electron states. Their inclusion does improve significantly the agreement between theory and experiment for both the low and high frequency edges of the spectrum. It does not appear sufficient, however, to resolve persistent deviations in the slow Lorentzian-like decay part of the spectrum in the intermediate 2-3 eV region.