Halogen bonding as a new driving force for layer-by-layer assembly

A new approach for fabricating a layer-by-layer polymer film was explored, which was based on the halogen bonding between poly(4-(4-iodo-2,3,5,6-tetrafluorophenoxy)-butyl acrylate) and poly(4-vinylpyridine). Layer-by-layer assembly of two polymers was confirmed by UV-vis spectroscopy and quartz crystal microbalance measurements. The interaction between the two polymers was identified as halogen bonding by X-ray photoelectron spectroscopy. The surface of the multilayer film is flat, and the thickness of one bilayer is about 1.3 nm. We also compared the stability of a halogen-bonded multilayer film in methanol with that of a hydrogen-bonded multilayer film.     

Halogen bonding in solution

Halogen bonding is the electron density donation based weak interaction of halogens with Lewis bases. Its applicability for molecular recognition processes long remained unappreciated and has so far mostly been studied in silico and in solid state. As most physiological processes and chemical reactions take place in solution, investigations in solutions are of highest relevance for its use in the pharmaceutical and material scientific toolboxes. Following a short discussion of the phenomenon of halogen bonding, this tutorial review presents an overview of the methods hitherto applied for gaining an improved understanding of its behaviour in solutions and summarizes the gained knowledge in order to indicate the scope of the techniques and to facilitate exciting future developments.     

Halogen bonding in supramolecular chemistry

Halogen bonding is the noncovalent interaction where halogen atoms function as electrophilic species. The energetic and geometrical features of the interaction are described along with the atomic characteristics that confer molecules with the specific ability to interact through this interaction. Halogen bonding has an impact on all research fields where the control of intermolecular recognition and self-assembly processes plays a key role. Some principles are presented for crystal engineering based on halogen-bonding interactions. The potential of the interaction is also shown by applications in liquid crystals, magnetic and conducting materials, and biological systems.     

Halogen bonding in DNA base pairs

Halogen bonding (R-X···Y) is a qualitative analogue of hydrogen bonding that may prove useful in the rational design of artificial proteins and nucleotides. We explore halogen-bonded DNA base pairs containing modified guanine, cytosine, adenine and thymine nucleosides. The structures and stabilities of the halogenated systems are compared to the normal hydrogen bonded base pairs. In most cases, energetically stable, coplanar structures are identified. In the most favorable cases, halogenated base pair stabilities are within 2 kcal mol(-1) of the hydrogen bonded analogues. Among the halogens X = Cl, Br, and I, bromine is best suited for inclusion in these biological systems because it possesses the best combination of polarizability and steric suitability. We find that the most stable structures result from a single substitution of a hydrogen bond for a halogen bond in dA:dT and dG:dC base pairs, which allows 1 or 2 hydrogen bonds, respectively, to complement the halogen bond.     

Halogen bonding in mono- and dihydrated halobenzene

Density functional theory calculations were performed on halogen-bonded and hydrogen-bonded systems consisting of a halobenzene (XPh; X = F, Cl, Br, I, and At) and one or two water molecules, using the M06-2X density functional with the 6-31+G(d) (for C, H, F, Cl, and Br) and aug-cc-pVDZ-PP (for I, At) basis sets. The counterpoise procedure was performed to counteract the effect of basis set superposition error. The results show halogen bonds form in the XPh-H₂O system when X > Cl. There is a trend toward stronger halogen bonding as the halogen group is descended, as assessed by interaction energy and X•••O_w internuclear separation (where O_w is the water oxygen). For all XPh-H₂O systems hydrogen-bonded systems exist, containing a combination of CH•••O_w and O_wH_w•••X hydrogen bonds. For all systems except X = At the X•••H_w hydrogen-bonding interaction is stronger than the X•••O_w halogen bond. In the XPh-(H₂O)₂ system halogen bonds form only for X > Br. The two water molecules prefer to form a water dimer, either located around the C H bond (for X = Br, At, and I) or located above the benzene ring (for all halogens). Thus, even in the absence of competing strong interactions, halogen bonds may not form for the lighter halogens due to (1) competition from cooperative weak interactions such as C H•••O and OH•••X hydrogen bonds, or (2) if the formation of the halogen bond would preclude the formation of a water dimer. © 2018 Wiley Periodicals, Inc.     

Key role of the polarization anisotropy of water in modeling classical polarizable force fields

We have evaluated the extent to which classical polarizable force fields, based either on the chemical potential equalization principle or on distributed polarizabilities in the framework of the Sum of Interactions Between Fragments Ab initio computed (SIBFA), can reproduce the ab initio polarization energy and the dipole moment of three distinct water oligomers: bifurcated chains, transverse hydrogen-bonded chains, and longitudinal hydrogen-bonded chains of helical shape. To analyze the many-body polarization effect, chains of different size, i.e., from 2 to 12 water monomers, have been considered. Although the dipole moment is a well-defined quantity in both classical polarizable models and quantum mechanical methods, polarization energy can be defined unequivocally only in the former type of approaches. In this study we have used the Kitaura-Morokuma (KM) procedure. Although the KM approach is on the one hand known to overestimate the polarization energy for strongly interacting molecules, on the other hand it can account for the many-body polarization effectively, whereas some other procedures do not. Our data show that, if off-centered lone pair polarizabilities are explicitly represented, classical polarizable force fields can afford a close agreement with the ab initio results, both in terms of polarization energy and in terms of dipole moment.     

Interaction between benzenedithiolate and gold: classical force field for chemical bonding

We have constructed a group of classical potentials based on ab initio density-functional theory (DFT) calculations to describe the chemical bonding between benzenedithiolate (BDT) molecule and gold atoms, including bond stretching, bond angle bending, and dihedral angle torsion involved at the interface between the molecule and gold clusters. Three DFT functionals, local-density approximation (LDA), PBE0, and X3LYP, have been implemented to calculate single point energies (SPE) for a large number of molecular configurations of BDT-1, 2 Au complexes. The three DFT methods yield similar bonding curves. The variations of atomic charges from Mulliken population analysis within the molecule/metal complex versus different molecular configurations have been investigated in detail. We found that, except for bonded atoms in BDT-1, 2 Au complexes, the Mulliken partial charges of other atoms in BDT are quite stable, which significantly reduces the uncertainty in partial charge selections in classical molecular simulations. Molecular-dynamics (MD) simulations are performed to investigate the structure of BDT self-assembled monolayer (SAM) and the adsorption geometry of S adatoms on Au (111) surface. We found that the bond-stretching potential is the most dominant part in chemical bonding. Whereas the local bonding geometry of BDT molecular configuration may depend on the DFT functional used, the global packing structure of BDT SAM is quite independent of DFT functional, even though the uncertainty of some force-field parameters for chemical bonding can be as large as approximately 100%. This indicates that the intermolecular interactions play a dominant role in determining the BDT SAMs global packing structure.     

Halogen bonding: a paradigm in supramolecular chemistry

The term halogen bonding describes the tendency of halogen atoms to interact with lone pair possessing atoms. The binding features and structural properties of halogen bonding are discussed and applied to drive the intermolecular self-assembly of hydrocarbons and perfluorocarbons in chemo-, site-, and enantioselective supramolecular synthesis. The halogen bonding is thus an effective and reliable tool in crystal engineering at the disposal of the supramolecular chemist.     

On the representation of force fields for chemically reacting systems

“Relaxed” force constants are uniquely defined for systems involving redundant coordinates, in contradistinction to the usual “rigid” force constants, and thus their use allows meaningful correlations to be made between force fields calculated for reactants, transition state, and product of a chemical reaction, for example formaldehyde hydration.     

Halogen bonding driven self-assembly of fluorocarbons and hydrocarbons

Halogen bonding, the interaction between halogen atoms and lone pair possessing atoms, overcomes the low affinity existing between perfluorocarbon and hydrocarbon derivatives. This interaction effectively drives the self-assembly of materials wherein perfluorocarbon and hydrocarbon layers alternate thanks to a remarkable module segregation. The approach has been exploited in supramolecular polymer coating.     

Halogen bonding in complexes of proteins and non-natural amino acids

In this work, we have analyzed the influence of halogen bonding to the stability of 44 complexes of proteins and non-natural amino acids. Fluorine- and chlorine-containing non-natural amino acids are more prevalent in the dataset, and an even larger number of contacts made by iodine-containing ligands are found. Only few halogen bonds with the hydroxyl oxygens and carboxylate side chains are found in the dataset. Halogen bonds with the nitrogen-containing side chains have higher occurrence than other acceptors. Backbone carbonyl oxygens and nitrogens are to a substantial extent involved in our dataset. We have observed a small percentage of interactions involving water as hydrogen bond donors. Additionally, most of the interacting residues comprising the interfaces also show a great degree of conservation. There is a clear interaction hot spot at distances of 3.5–3.7 Å and Θ₁ angles of 100–120°. There is also a cluster of contacts featuring short distances (2.6–2.9 Å) but only nearly optimal Θ₁ angles (140–160°). 51.3% of stabilizing residues are involved in building halogen bonds with the non-natural amino acids. We discovered three types of structural motifs significantly over-represented: beta-turn-ir, beta-turn-il and niche-4r. The halogen-bonding statistics of the dataset do not show any preference for α-helices (36%), β-sheets (36%), or turns/coils (28%) structures. Most of the amino acid residues that were involved in halogen bonds prefer to be in the solvent excluded environment (buried). Furthermore, we have shown that in amino acid–protein complexes halogen atoms can sometimes be involved in hydrogen bonding interactions with hydrogen bonding-donors. The results from this study might be used for the rational design of halogenated ligands as inhibitors and drugs, and in biomolecular engineering.     

Development of a classical force field for the oxidized Si surface: application to hydrophilic wafer bonding

We have developed a classical two- and three-body interaction potential to simulate the hydroxylated, natively oxidized Si surface in contact with water solutions, based on the combination and extension of the Stillinger-Weber potential and of a potential originally developed to simulate SiO(2) polymorphs. The potential parameters are chosen to reproduce the structure, charge distribution, tensile surface stress, and interactions with single water molecules of a natively oxidized Si surface model previously obtained by means of accurate density functional theory simulations. We have applied the potential to the case of hydrophilic silicon wafer bonding at room temperature, revealing maximum room temperature work of adhesion values for natively oxidized and amorphous silica surfaces of 97 and 90 mJm(2), respectively, at a water adsorption coverage of approximately 1 ML. The difference arises from the stronger interaction of the natively oxidized surface with liquid water, resulting in a higher heat of immersion (203 vs 166 mJm(2)), and may be explained in terms of the more pronounced water structuring close to the surface in alternating layers of larger and smaller densities with respect to the liquid bulk. The computed force-displacement bonding curves may be a useful input for cohesive zone models where both the topographic details of the surfaces and the dependence of the attractive force on the initial surface separation and wetting can be taken into account.     

Halogen bonding in polymer science: towards new smart materials

The halogen bond is a special non-covalent interaction, which can represent a powerful tool in supramolecular chemistry. Although the halogen bond offers several advantages compared to the related hydrogen bond, it is currently still underrepresented in polymer science. The structural related hydrogen bonding assumes a leading position in polymer materials containing supramolecular interactions, clearly indicating the high potential of using halogen bonding for the design of polymeric materials. The current developments regarding halogen bonding containing polymers include self-assembly, photo-responsive materials, self-healing materials and others. These aspects are highlighted in the present perspective. Furthermore, a perspective on the future of this rising young research field is provided.

The incorporation of halogen bonding into polymer architectures is a new approach for the design of functional materials. This perspective emphasizes the current development in the field of halogen bonding featuring polymer materials.     

Simulations of vibrational spectra from classical trajectories: calibration with ab initio force fields

An algorithm allowing simulating vibrational spectra from classical time-dependent trajectories was applied for infrared absorption, vibrational circular dichroism, Raman, and Raman optical activity of model harmonic systems. The implementation of the theory within the TINKER molecular dynamics (MD) program package was tested with ab initio harmonic force fields in order to determine the feasibility for more extended MD simulations. The results suggest that sufficiently accurate frequencies can be simulated with integration time steps shorter than about 0.5 fs. For a given integration time step, lower vibrational frequencies ( approximately 0-2000 cm(-1)) could be reproduced with a higher accuracy than higher-frequency vibrational modes (e.g., O-H and C-H stretching). In principle, the algorithm also provides correct intensities for ideal systems. In applied simulations, however, the intensity profiles are affected by an unrealistic energy distribution between normal modes and a slow energy relaxation. Additionally, the energy fluctuations may cause weakening of the intensities on average. For ab initio force fields, these obstacles could be overcome by an arbitrary normal mode energy correction. For general MD simulations, averaging of many shorter MD trajectories started with randomly distributed atomic velocities provided the best spectral shapes. alpha-pinene, D-gluconic acid, formaldehyde dimer, and the acetylprolineamide molecule were used in the tests.     

Atom-atom partitioning of total (super)molecular energy: the hidden terms of classical force fields

Classical force fields describe the interaction between atoms that are bonded or nonbonded via simple potential energy expressions. Their parameters are often determined by fitting to ab initio energies and electrostatic potentials. A direct quantum chemical guide to constructing a force field would be the atom-atom partitioning of the energy of molecules and van der Waals complexes relevant to the force field. The authors used the theory of quantum chemical topology to partition the energy of five systems [H2, CO, H2O, (H2O)2, and (HF)2] in terms of kinetic, Coulomb, and exchange intra-atomic and interatomic contributions. The authors monitored the variation of these contributions with changing bond length or angle. Current force fields focus only on interatomic interaction energies and assume that these purely potential energy terms are the only ones that govern structure and dynamics in atomistic simulations. Here the authors highlight the importance of self-energy terms (kinetic and intra-atomic Coulomb and exchange).     

Halogen bonding and pi...pi stacking control reactivity in the solid state

The halogen bonding and the pi...pi stacking interactions induce the noncovalent self-assembly of modules into photoreactive supramolecular architecture. The pi...pi interaction pre-organizes the template, and the halogen bonding aligns the olefins to conform to the topochemical principle for photoreaction. The UV irradiation of the crystal resulted in a cyclization product with quantitative yield and stereospecificity.     

Halogen bonding with the halogenabenzene bird structure,halobenzene, and halocyclopentadiene

The ability of the “bird-like” halogenabenzene molecule, referred to as X-bird (XCl to At), to form halogen-bonded complexes with the nucleophiles H₂O and NH₃ was investigated using double-hybrid density functional theory and the aug-cc-pVTZ/aug-cc-pVTZ-PP basis set. The structures and interaction energies were compared with 5-halocyclopenta-1,3-diene (halocyclopentadiene; an isomer of halogenabenzene) and halobenzene, also complexed with H₂O and NH₃. The unusual structure of the X-bird, with the halogen bonded to two carbon atoms, results in two distinct σ-holes, roughly at the extension of the C-X bonds. Based on the behavior of the interaction energy (which increases for heavier halogens) and van der Waals (vdW) ratio (which decreases for heavier halogens), it is concluded that the X-bird forms proper halogen bonds with H₂O and NH₃. The interaction energies are larger than those of the halogen-bonded complexes involving halobenzene and halocyclopentadiene, presumably due to the presence of a secondary interaction. © 2019 Wiley Periodicals, Inc.