The ultimate goal of materials chemistry is the understanding of the macroscopic properties of materials in terms of the microscopic molecular interactions. This is a common theme in all of the natural sciences. The differences between the various scientific disciplines (chemistry, biology, physics, materials sciences, etc.) often disappear as the traditional macroscopic phenomenology is replaced by a more molecular approach. At present, most of what is known about the chemical and physical properties of materials is still largely empirical, especially in the case of amorphous, macromolecular or interfacial systems. The development of a microscopic theoretical description for a variety of such complex systems is the primary focus of our research. Several representative projects are listed below. They are all projects in which great advantage will be gained by exploiting the natural symbiosis between analytical and computer-simulation techniques.

Snapshot from a molecular dynamics
simulation of a premelting at a Al-Pb solid-liquid interface
[Phys. Rev. Lett. 110, 096102 (2013)]


I. Structure, thermodynamics and kinetics of condensed phase interfaces:

Understanding the structure, thermodynamics and kinetics of the interfaces between condensed phase systems is of paramount importance in a wide array of technologically important fields, including crystal growth and nucleation, wetting, material strength and stability, Experimental study is often difficult as such interfaces lie sandwiched between two dense phases. Because experimental data is lacking or often incomplete, molecular modeling studies of such systems can offer valuable insight into their detailed molecular phenomenology. Much of our work in this area in the past has focused on the solid-liquid interfaces of simple, one-component model systems. In recent years, our focused has shifted to interfaces in complex systems such as multicomponent systems (e.g., metal alloys), systems with molecular interactions and grain boundaries. For such studies, we combine state of the art molecular simulation techniques and applied statistical mechanics to a variety of realistic and model systems. Recent publications include

  1. "Thermodynamics and Intrinsic Structure of the Al-Pb Liquid-Liquid Interface: A Molecular-Dynamics Simulation Study", Y. Yang and B.B. Laird, J. Phys. Chem. B (Jim Skinner Festschrift), 118 8373-8380 (2014).
  2. "Calculation of the interfacial free energy of a binary hard-sphere liquid at a planar hard wall", J.L. Kern and B.B. Laird, J. Chem. Phys., 140, 024703:1-7 (2014).
  3. "Solid-liquid interfacial premelting", Y. Yang, M. Asta and B.B. Laird, Phys. Rev. Lett. 110, 096102 (2013).
  4. "Interfacial free energy of a hard-sphere fluid in contact with curved hard surfaces", B.B. Laird, A. Hunter and R.L. Davidchack, Phys. Rev. E, 86 060602:1-5 (2012).

II. Phase-Equilibrium and Transport properties of Novel Solvent Media for Environmentally Beneficial Catalysis:

For the last several years we have been working closely with KUs Center for Environmentally Beneficial Catalysis, an engineering research center whose mission is the development and application of environmentally friendly catalytic processes for industry. Much of our effort has gone into molecular modeling studies to support the development of novel  "green" solvent media for catalysis. One class of solvents that is of interest are so called gas-expanded liquids (GXLs), in which a large amount of gas (usually carbon dioxide) is dissolved into a traditional organic solvent, expanding its volume appreciably. The resulting solvents are environmentally beneficial in that they reduce the amount of potentially harmful organic solvents necessary to carry out a particular catalytic reaction. These solvent systems also have enhanced mass transport and reactant solubility properties over the neat organic solvents and over super-critical carbon dioxide. We have been using a number of molecular simulation techniques (Gibbs Ensemble and Grand Canonical Monte Carlo and Molecular Dynamics Simulation) to study the phase-equilibrium, structure and transport properties of these novel solvents.

Image of a Monte Carlo pore-filling calculation for a silica mesopore
with hydrophilic OH-groups, illustrating the concentration and
distribution of ethylene and solvent within the pore

Recently, we have entered into a new collaboration with the CEBC funded by the National Science Foundation program Networks for Materials Synthesis and Design. This is a highly multidisciplinary and multiscale effort involving several KU research groups: two in Chemical Engineering (Subramaniam and Choudhary), two in Computational Chemistry (Laird and Thompson) and one Organic Chemistry group (Tunge).This highly collaborative effort uses molecular modeling, synthesis, reactor engineering, kinetic studies and life-cycle assessments to develop and optimize the design of industrial catalytic processes for minimal environmental impact. This framework is being applied to two examples of high-volume chemical processes. One is a safer, phosgene-free route to dimethyl carbonate, which is used to manufacture polycarbonates and polyurethanes. The other is a cleaner, more efficient route from butadiene to adipic acid, used to make nylon. The discoveries and methodologies developed for these two chemicals will also find broader impacts for other major chemical processes. Our theoretical/computational modeling will provide a detailed understanding of the molecular processes involved in the complex reactions, which are often inaccessible through experiment alone.  For example, in a new synthetic route to DMC, ethylene is oxidized to ethylene oxide inside a catalytic metal-exchanged amorphous silica mesopore.  Molecular simulation allows us to measure the nano-confined reactant system more easily than typical experiments allow.  Generally, we are able to study the phase behavior, structure, and transport properties of these novel reaction environments.

  1. "Prediction of phase equilibria and transport properties in carbon-dioxide expanded solvents by molecular simulation", Y. Houndonougbo, B.B.Laird and K. Kuczera, Mol. Simul. 33, 861-869 (2007).
  2. "Transport properties in carbon dioxide-expanded acetonitrile from molecular-dynamics simulations", Y. Houndonougbo, B.B. Laird and K. Kuczera", J. Chem. Phys., 126 074507 (8 pages) (2007).
  3. "Phase Equilibria in Carbon Dioxide Expanded Solvents: Experiments and Molecular Simulations", Y, Houndonougbo, J. Hong, B. Rajagopalan, K. Wong, K. Kuczera, B. Subramaniam, and B.B. Laird, J. Phys. Chem. B 110, 13195-13202 (2006).
  4. "Monte Carlo Simulations of carbon-dioxide-expanded Acetonitrile", Y. Houndonougbo, J-X Guo, G.H. Lushington and B. Laird, Mol. Phys. 104, 2955--2961 (2006).

III. Development of new algorithms for molecular modeling

Molecular dynamics computer simulation has become an invaluable tool in chemistry, chemical engineering, physics,materials science and biology; however, its uses are still limited by the relatively small system sizes and short time scales that can be simulated atpresent. Progress in this area therefore comes from advances in computer technology and in the development of efficient and stable algorithms. The latter is the goal of an  multidisciplinary project involving the Laird group in collaboration with Prof. Ben Leimkuhler, an applied mathematician at the University of Edinburgh.

One central property of Hamiltonian systems used in molecular dynamics simulations is that it is sympletic, that is, the sum of oriented areas in phase space is a conserved quantity. A well- known consequence of this property is the conservation of phase space volume in Hamiltonian dynamics - otherwise known as Liouville's theorem. It has been shown that stable numerical integrators for molecular dynamics can be constructed by requiring that these approximate integrators are themselves sympletic. We have applied symplectic numerical integration techniques to a number of problems in molecular simulation ranging from rigid body motion to extended Hamiltonian methods for constant-temperature and/or constant pressure molecular dynamics. Representative publications are

  1. "Generating Generalized Distributions from Dynamical Simulation", E.J. Barth, B.B. Laird and B.J. Leimkuhler, J.Chem. Phys., J. Chem. Phys. 118 5759-5767 (2003)
  2. "Generalized dynamical thermostatting", B.B. Laird and B.J. Leimkuhler, Phys. Rev. E 68, 016704 (2003).
  3. "Constant-temperature molecular-dynamics algorithms for mixed hard-core/continuous potentials" , Y.A. Houndonougbo and B.B. Laird, J. Chem. Phys. 117, 1001 ( 2002).
  4. "Symplectic algorithm for constant-pressure molecular-dynamics using a Nose-Poincare thermostat" J.B. Sturgeon and B.B. Laird, J. Chem. Phys. 112 , 3474--3482 (2000)
  5. "Molecular dynamics algorithms for mixed hard-core/continuous potentials'' Y.A. Houndonougbo, B.B. Laird and B.J. Leimkuhler, Mol. Phys., 98 309--316 (1999).
  6. "The Nose-Poincare method for constant temperature molecular dynamics'', S.D. Bond, B.J. Leimkuhler and B.B. Laird, J. Comp. Phys. 151, 114-134 (1999).

Laird Group News

August 5-11, 2017 Zeke Piskulich attended the Gordon Research Seminar (GRS) and Conference (GRC) on the Physics and Chemistry of Liquids and was elected co-chair of the 2019 GRS on this topic. 

July 27, 2017 Summer 2017 Laird Group REU student Trevor Neal (Westminster College) presented a poster entitled "Increasing the efficiency of electrostatic calculations for vapor-liquid equilibnrium at the 2017 KU Undergraduate Research Poster Session. [Photo]

July, 2017 Seth Martin attended the CCP5 Molecular Dynamics Workshop in Lancaster England. 

June 2017 Former graduate student, Isaak Daniels, started a postoc at McGill University in Montreal, Canada in the group of Prof. Alejandro Rey in Chemical Engineering.

June 2017  Brian Laird visited former postdoc Yang Yang at East China Normal University (ECNU) in Shanghai, China and gave talks in both Physics and Chemsitry Departments at ECNU and Jiao Tong University. [Photo]

Dec. 15, 2016  Seth Martin has joined the group as a new graduate student. Welcome!

Dec. 15, 2016  Dr. Isaak Daniel's paper with Zhenxing Wang and Brian Laird "Dielectic Properties of Organic Solvents in an Electric Field" has been accepted to the Journal of Physical Chemistry C. 

August 2016 Dr. Pablo Palafox-Hernandez's paper with Brian Laird "Orientation dependence of heterogeneous nucleation at the Cu–Pb solid-liquid interface" was accept to the Journal of Chemical Physics. 

July 4, 2016  Dr. Pansy Patel's paper with Brian Laird and Ward Thompson "A Density Functional Theory Study of Ethylene Epoxidation Catalyzed by Niobium-doped Silica" has been accepted in the Journal of Molecular Catalysis. 

June 2016  Brian Laird has been appointed as an External Senior Fellow of the Freiburg Institute for Advanced Studies in Freiburg, Germany. He will be there during his sabbatical in Summer and Fall 2016. He will also be working with the group of Roland Roth at the University of Tübingen Physics Department

May 23, 2016 Our paper "Electric potential calculation in molecular simulation of electric double layer capacitors" by Zhenxing Wang, David Olmsted, Mark Asta and Brian B. Laird has been accepted to a Special Issue on the Physics of Supercapacitors and Electroactuators in J. Phys.: Condensed Matter. 

May 4, 2016 Isaak Daniels successfully defended his Ph.D. dissertation "Energy Applications of Molecular Simulation"

March 2016  Dr. Peter Barry has joined the group to work on an NSF funded project on chemically heterogeneous solid-liquid interfaces.

March 2016 B. Laird's paper with Ruslan Davidchack and Roland Roth "Hard spheres at a planar hard wall: Simulations and density functional theory" [Cond. Matt. Phys. 19 23001:1-10 (2016)]  was published as part of  the Tony Haymet Festschrift.  

March 2016 Our paper on "Molecular simulation of ethylene-expanded methanol: Phase behavior, structure, and transport properties" has been published in Fluid Phase Equilibria (Published online Dec. 12, 2015). 

February 2016 Our paper "Tunability of Gas-Expanded Liquids under Confinement: Phase Equilibrium and Transport Properties of Ethylene-Expanded Methanol in Mesoporous Silica" [J. Phys. Chem. C, 120, 5010-5019 (2016)] was published online. 

KU Today
One of 34 U.S. public institutions in the prestigious Association of American Universities
44 nationally ranked graduate programs.
—U.S. News & World Report
Top 50 nationwide for size of library collection.
23rd nationwide for service to veterans —"Best for Vets," Military Times