
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

71793
10.1038/s41598-024-71793-6
Article
Exploring hydrophobicity or hydrophilicity of borophene surface via reactive molecular dynamics simulation
Foroutan Masumeh foroutan@ut.ac.ir

Sababkar Mahnaz
Bavani Borhan Mostafavi
https://ror.org/05vf56z40 grid.46072.37 0000 0004 0612 7950 Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran, Tehran, Iran
13 9 2024
13 9 2024
2024
14 2143611 6 2024
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Borophene, a novel two-dimensional material unveiled in 1998, has garnered significant interest among researchers due to its distinct mechanical and electrical characteristics. Efforts to experimentally synthesize borophene continue to captivate researchers’ interest in recent years. Given the current lack of experimental studies on the interaction between water and the borophene surface, molecular dynamics simulation offers a valuable approach for predicting the substance’s reactivity with water. Additionally, such simulations can assess the hydrophilicity and hydrophobicity of borophene, providing valuable insights into its properties. In our current research, we utilized reactive molecular dynamics simulation to investigate the wetting behavior of borophene. Our findings reveal that the borophene surface exhibits hydrophobic characteristics, demonstrating anisotropic wettability. Specifically, the water contact angle was calculated to be 149.11° along the zigzag direction and 148.4° along the armchair direction. The contour map of the interaction energy between a water molecule and the borophene surface revealed a notable energy barrier in the zigzag direction. This barrier contributes to the asymmetric spreading of the water droplet on the surface. Density profiles and radial pair distribution function (RDF) diagrams of the water droplet on the borophene surface further corroborated the hydrophobic nature of borophene by indicating a significant distance between the water droplet and the surface. Moreover, analysis of the number of hydrogen bonds demonstrated that borophene efficiently utilizes nearly all its capacity to form hydrogen bonds. Additionally, we compared the wettability of borophene with that of other two-dimensional materials, such as various graphene allotropes and phosphorene, which have been subjects of recent investigation.

Keywords

Borophene
Wettability
Reactive molecular dynamic simulation
Subject terms

Physical chemistry
Surface chemistry
Theoretical chemistry
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Borophene is a two-dimensional surface of boron atoms1–3. Borophene has been the focus of a large number of researchers due to its various applications in various sciences. In comparison to other two-dimensional materials, borophene is highly significant, with many of its allotropes predicted by quantum studies prior to their experimental synthesis2,4–7. Finally in 2007, a few years after the prediction of the first structure of borophene surface, this structure was modified, and more stable structures of it were introduced7–9. boron is the fifth element in the periodic table, and its ground-state configuration is [He]2s22p1, which makes this atom to not follow the octet rule because the number of atomic orbitals are largen than valance electrons. This causes the boron atom to use unconventional three-center-two-electron bonds (3c-2e) besides the conventional two-center-two-electron bonds (2c–2e) in order to reach an equilibrium structureuse; as a result, this fact leads borophene surface to the formation of various allotrops10,11.

During recent years, some allotropes of borophene surface were experimentally synthesized on metal substrates such as Ag5,12–14, Al15, Au16, Cu14, Ir17. Recently, some of the borophene's mechanistic features18–20, optical transparency20,21, thermal conductivity19,22–25, biomedical applications26,27, and electrical properties28 are investigated. Additionally, sundry studies convey that borophene surface is a great candidate to design batteries29–32 and gas sensors33–35. Despite the extensive research on borophene’s properties, its wettability characteristics remain unexplored.

First-hand information on a surface hydrophilicity or hydrophobicity can be valuable across various industries such as self-cleaning surfaces36,37, anti‐icing properties38, development of super waterproof coatings39,40, non-stick material41, anti-fouling42, vapor chamber43, thermal diodes44, Separation of oil from water45, anti-corrosion coating46,47, and so forth.

Calculating the contact angle experimentally is challenging because the existence of impurities and defects on the surface is inevitable; as a result, this changes the behavior of the water droplet48. Molecular dynamics simulations provide a robust framework for examining the wettability of two-dimensional materials. Molecular dynamics simulation not only enables researchers to accurately investigate changes at the molecular level, but also macroscopic parameters such as density, surface tension, viscosity, contact angle, and layering mechanism can be studied in more detail49–51.

In the present work, we will study the wettability of borophene with molecular dynamics simulation approach. For this purpose, first, the contact angle of a water drop will be calculated with the help of drawn contour maps. In the following, the energy contour map of the borophene surface, the density profile diagrams of the water droplet, the RDF diagram, the number of hydrogen bonds and the way the water droplet moves on the borophene surface will be studied and analyzed.

Computational method

Reactive force field

Reactive molecular dynamics simulation is a bridge between the world of quantum calculations and the experimental force field. The reactive force field can model chemical reactions based on the dynamic distribution of charges and the concept of bond length. Equation (1) shows the total energy that the reactive force field calculates for the system.1 Esystem=Ebond+Elp+Eover+Eunder+Eval+Epen+Ecoa+EC2+Etriple+Etors+Econj+EHbonds+Evdwaals+Ecoulomb

In Eq. (1), the first term calculates the bond energy, and the second term calculates the Lone pair electron energy. In the following, the third and fourth terms calculate over-coordination and under-coordination energies for the system, respectively. Also, expressions Eval and Epen are used to calculate angular energy and penalty energy, respectively. Ecoa is the term related to three-body conjugation energy calculation. Ec2 and Etriple are used to calculate the energy correction for the C2 molecule and the triple bond energy, respectively. Moreover, Etors is the term related to the calculation of torsion angle energy and Econj corresponds to the conjugation effect term. Other terms are related to the calculation of hydrogen energy and include the contribution of hydrogen bond, van der Waals, and Coulomb energies respectively. To know more details about the Rexx force field, refer to references52,53. In the current work, Reax force field reference54 is used to study interatomic interactions.

Simulation details

In the present work, we investigate the wetting behavior of the two-dimensional borophene surface with the helping hand of molecular dynamics simulation. All the simulations have been performed with a laerge-scale atomic/molecular massively parallel simulator (LAMMPS) package55, and structures illustrated and photos taken by Visual Molecular Dynamics (VMD) software56. Periodic boundary conditions were considered periodic in all directions except the Z direction. Boron atoms are considered fixed during the simulation. The temperature of this simulation is 300 K in which the Nose–Hoover thermostat has been used in order to keep the temperature constant during the simulation. It goes without saying that the simulation is performed using an NVT ensemble with a time step of 0.2 fs. Besides, the total simulation time was 1.5 ns; this time was enough to the droplet to reach equilibrium. Also, a droplet with 500 water molecules was used in this simulation.

Details of borophene and water drop configuration

In this section, borophene surface and the initial configuration of the water drop on this surface will be introduced. The two lattice constants of this structure of borophene are 1.614 Å and 2.866 Å and its out-of-plane height (h) in the unit cell is 0.911 Å57. Figure 1 shows the surface and unit cell of borophene.Fig. 1 Structure of borophene surface along the Zig-zag and Armchair direction and the position of atoms in the unit cell.

According to Fig. 1, the borophene surface consists of two sub-layers, top and bottom. In this figure, the light and dark purple spheres represent atoms in the top and bottom sublayer, respectively. Also, the X and Y direction of the borophene surface has an armchair and zigzag structure, respectively.

Among the two-dimensional structures, phosphorene and penta-graphene are examples that, be close resemblance to borophene, are composed of sub-layers. The arrangement of atoms, the out-of-plane height, and the distance between atoms in each of these structures are different. For example, the out-of-plane height of phosphorene is 2.246 Å and penta-graphene is 1.20 Å.

To study the wettability of the borophene surface, a drop of water is placed in the center of this surface. Figure 2 shows the initial configuration of a water droplet on the borophene surface.Fig. 2 Initial configuration of water drop on borophene surface from (a) top view, (b) zig-zag direction, and (c) armchair direction of borophene. (Purple, red, and white spheres represent boron, oxygen, and hydrogen atoms, respectively).

Figure 2a shows the initial configuration of a water droplet on the borophene surface from the top view. In this figure, the radius of the water drop is 16 Å and the size of the borophene surface is 79 × 79 Å2. On the other hand, the water droplet is 23 Å away from the edges of the borophene surface. Figure 2b and c show the initial configuration from the zigzag and armchair directions of the surface. In these images, the distance of the water drop from the surface of borophene is 3 Å.

Contact angle measurement methods

Wettability is a substantial indicator to identify the properties of solid materials. In fact, wettability shows the degree of interaction and adhesion of a liquid such as water in contact with the surface of a solid. Calculating the contact angle of a liquid drop on a solid surface determines the degree of wetting of the solid surface. The Yang equation is used to measure the contact angle in order to describe the wettability properties58:2 cosθ=γSV-γSLγLV

In Eq. (2), γsv, γsl, and γlv represent the surface tensions of the solid–vapor, solid–liquid, and liquid–vapor surfaces, respectively. θ indicates the contact angle of the droplet on the solid surface.

There are assorted methods for calculating the contact angle of a water drop, one of them is fitting a circle. In this method, first of all, the water drop is meshed into small cubes, and then, the number of water molecules in each mesh is calculated. Accordingly, the density of water molecules at any point of the drop can be calculated. After calculating the density of water molecules in the droplet, parts of the water drop with a density of 0.4 g/cm3 are considered as the boundary between the liquid and vapor phases59. Then a circle is fitted on these points as the boundary between the two phases; and the angle of the circle with the surface is calculated. It should be noted that to follow this method, the coordinates of the atoms in the last 300 frames; where the water droplet reached equilibrium, were taken as output from the simulation, and then an average of their coordinates was used to calculate the contact angle.

Results and discussions

In this section, the wetting behavior of the borophene surface is investigated. For this purpose, first, in section “Contact angle”, the contact angle of the water droplet will be examined with the help of a density contour map. Then, in section “Energy contour map of borophene surface”, the surface energy contour map will be studied, and in the following sections, the density profile, the number of hydrogen bonds, and the RDF graph will be studied respectively.

Contact angle

In order to calculate the contact angle of the water drop on the borophene surface, we need to know how the atoms are placed in the equilibrium at the end of the simulation. Figure 3 shows the arrangement of water molecules on the borophene surface in an equilibrium state.Fig. 3 Equilibrium configuration of the water droplet on Borophene surface from (a) up view, (b) zigzag direction, and (c) armchair direction. (Purple, red, and white spheres represent boron, oxygen, and hydrogen atoms, respectively).

Figure 3 shows snapshots of the water droplet on the borophene in the equilibrium state. In order to be able to calculate the contact angle of the water drop on borophene, it is necessary to draw the density contour map of water molecules. In this regard, the water droplet has meshed with the size of 1.5 Å2, and then the number of water molecules in each meshing is calculated. Finally, each mesh of water droplet with the same density is marked with the same color. Figure 4 shows the density contour map of the water droplet.Fig. 4 Density contour map of water droplet from the (a) armchair direction and (b) zigzag direction.

According to the color bar in Fig. 4, the accumulation of particles in the center of the drop is more than in other areas, and by moving from the center of the drop to the outer shell of the it, the density of the particles decreases. In essence, according to the color bar in Fig. 4, the density of the water droplet in the center is 1, but as it moves away from the center of the drop, the density decreases and finally reaches zero.

The red circle, which is drawn from the density of 0.4 g/cm3 of water, in Fig. 4 defines the boundary between the liquid and the vapor phase. After that, the angle between the red circle and the borophene surface is calculated. This angle is 149.11° in the zigzag direction and 148.4° in the armchair direction. The contact angle of the water drop shows that the borophene surface is hydrophobic; to put it simply, water molecules do not have much tendency to spread on the borophene surface. The hydrophobicity of the borophene surface was also reported in previous studies60, but the exact contact angle of the water droplet was not investigated.

Additionally, density contour maps in Fig. 4 show that the distribution of water molecules near the surface is different along the zigzag and armchair directions. According to Fig. 4, in the armchair direction, the water molecules have a regular distribution near the surface, while in the zig-zag direction, have taken on a zig-zag structure to some extent. This behavior has not been observed in the central parts of the drop, which have a greater distance from the borophene surface. It seems that the reason for the erratic distribution of water molecules in the zigzag direction of the borophene surface is that the energy barrier in this direction is greater than The energy barrier in the armchair direction30. This causes the water droplet to be elliptical; so the contact angle of the water molecules in the zigzag and armchair directions is different.

Previous studies have shown that phosphorene, similar to borophene, has anisotropic wettability. unlike borophene, the contact angle of the water droplet on phosphorene in the zigzag direction is lower than in the armchair direction. Based on this study, the contact angle of a water drop on phosphorene61 in the zigzag and armchair direction is 130° and 132°, respectively. Phosphorus, vice versa borophene surface, has the barrier energy in the armchair direction. In layman's terms, the tendency of the water drop to spread in the armchair is less than the zigzag direction. The contact angle of water drop on penta-graphen62 for both directions is reported to be 134.63°, which shows that this surface does not have anisotropic behavior. Comparing the wetting behavior of the three surfaces of borophene, phosphorene, and penta-graphene shows that the most hydrophobic surface is borophene, phosphorene, and penta-graphene respectively.

Energy contour map of borophene surface

Since the unevenness of the borophene surface can affect the placement of water molecules, the energy profile of the borophene surface has been investigated in this section. Figure 5 shows the interaction energy contour map of a water molecule on the surface of borophene.Fig. 5 Interaction energy contour map of a water molecule on the borophene surface along the zigzag and armchair direction.

To draw the contour map shown in Fig. 5, a water molecule was scanned on the borophene surface and its interaction with each boron atom in the substrate is calculated. According to Fig. 5, different positions of the borophene surface have different interaction energies with water molecules. According to the color bar in Fig. 5, the interaction of a water molecule was strong in red regions and weak in blue regions.

As it is clear in this figure, the borophene surface has an energy barrier in the zigzag direction. Similar to previous studies, there are two different energy barriers in the zigzag direction, both of which are stronger than the energy barrier along the armchair direction30. This energy barrier has caused water molecules to have less tendency to spread in a zigzag direction. The difference in the energy barrier along the zigzag and armchair directions has led to the anisotropic wetting behavior of the borophene surface. This behavior is consistent with the contact angle studies that were done in the previous section. In the previous section, it was observed that the contact angle of the water drop on the borophene surface along the zigzag direction is lower than it is along the armchair direction.

Generally, spreading water droplets on two-dimensional surfaces is affected by their structure. For example, the diagram of the interaction energy of a water molecule on the phosphorene surface shows that the water molecule has a strong interaction with the phosphorus atoms in the armchair direction. To make it clearer, the energy barrier in the armchair phosphorene direction is higher than the energy barrier in the zigzag direction. As a result, the presence of this energy barrier in phosphorene gives less movement freedom to water molecules in this direction. Finally, this has caused the reduction of the contact angle along the armchair direction. Unlike borophene and phosphorene, in penta-graphene, the distance between the energy barriers is large, so the water drop can overcome the energy barrier and spread evenly in two directions.

Density profile

The equilibrium configuration of water drop on the borophene surface is shown in Fig. 3. Figure 6 shows the density profile diagram of this water droplet according to its height in the Z direction.Fig. 6 Density profile of water droplet on borophene surface along the Z direction.

In this figure, the highest atom of the substrate is considered as a reference atom for drawing the density profile diagram. According to Fig. 6, there is no separate layer structure of water drops on the borophene surface. The layering of water drops is tangled and overlapping. The height of the water drop on the borophene surface is about 27 Å. Moreover, the first peak of the density profile of water molecules on the borophene surface has been seen at 4 Å. In essence, the water drop is at a distance of 4 Å from the surface. Penta-graphene density profile diagrams62 show that water layers are completely separate from each other and the number of layering is much more than the layering of a water drop on the borophene surface.

Hydrogen bond

One of the most important parameters in determining the spatial structure of a water drop is calculating the number of hydrogen bonds. This parameter can be calculated both experimentally and using quantum mechanics. The number of hydrogen bonds was calculated according to the method reported in reference63. There is a hydrogen bond between two molecules i and j if first of all, the distance between the two oxygen atoms is less than 3.5 Å, second, the distance between the hydrogen atom i and the oxygen atom j is less than 2.45 Å, and finally, the angle between oxygen i and the OH bond Molecule j be less than 30°. Figure 7 shows the number of hydrogen bonds between water molecules on the borophene surface in the last 300 fs.Fig. 7 Number of hydrogen bond between water molecules relative to the time.

Figure 7 shows the number of hydrogen bonds in a water drop during 1.5 ns. The number of hydrogen bonds is affected by the surface structure. To make it clearer, water molecules form about 835 hydrogen bonds on the borophene surface. Since the number of hydrogen bonds that a water molecule can have is almost 2, it can be seen that the water drop on the borophene surface has almost used its full power to form hydrogen bonds. A water droplet with a higher number of hydrogen bonds between its molecules is more cohesive and spherical compared with other water droplets; as a result, its contact angle will increase. The high number of hydrogen bonds between water molecules indicates the hydrophobicity of the borophene surface.

The number of hydrogen bonds reported for a water molecule on the phosphorene surface is approximately 80061. On the other hand, the diagrams of hydrogen bonds of penta-graphene62 show that the water drop on penta-graphene has used a quarter of its container to form hydrogen bonds. Needless to say, the presence of more hydrogen bonds in the drop increaces the contact angle of the water drop with the surface63. This fact is evident in the contact angle of water drops on borophene, phosphorene, and penta-graphene.

Radial distribution function diagram

One of the other parameters that help to have a better understanding of droplet behavior is the radial distribution function diagram. Accordingly, we can understand how the atoms are arranged and placed in relation to each other. The diagram in Fig. 8 shows the radial distribution function of B–O, O–O, and O–H atoms.Fig. 8 Radial distribution function (RDF) analysis. (a) Full RDF curves for Boron-Oxygen (B-O) and Boron–Hydrogen (B–H) interactions. (b) Zoom-in of the initial part of the RDF curves from (a), and (c) RDF curves for Oxygen–Oxygen (O–O) and Oxygen-Hydrogen (O–H) interactions.

The radial distribution function is one of the most important features of determining the structure of matter. Based on this function, the probability of finding an atom around a hypothetical central atom can be calculated. According to Fig. 8a, which is related to boron-oxygen and boron hydrogen atoms, shows that the oxygens are far away from the boron atom. This issue is consistent with the water droplet distance seen in the density profile graphs. There are two minor peaks in the initial part of this diagram for B–O and B–H; the magnification of these peaks is shown in Fig. 8b. The place of appearance of peaks B–H and B–O is almost the same, but the peak intensity of B–O is almost twice the peak of B–H. This shows that the first layer of water molecules is formed at this distance from the borophene surface, and also the water molecules in this layer are facing the borophene surface from their oxygen side.

Three main peaks have been observed for the radial distribution function of oxygen-hydrogen atoms in Fig. 8c. The first peak of this graph occurred at 0.9 Å. This peak shows the bond length of oxygen-hydrogen atoms in the water molecule. The second peak of this graph occurred approximately at 1.9 Å, representing the length of the hydrogen bonds between water molecules. The third peak is related to the hydrogen bond of molecules that are located at a further distance from the hypothetical central atom.

Two relatively sharp peaks and one small peak have also been observed in the red diagram, which is related to oxygen–oxygen atoms. The first peak appeared approximately at 2.9 Å. This means that water molecules are facing each other by their hydrogen atoms. The second peak of this diagram is located at a distance of 3.4 Å from the hypothetical central atom. The presence of hydrogen bonds between water molecules makes oxygen atoms unable to move freely.

The large distance of oxygens from boron atoms causes low interaction between oxygen and boron atoms, and as a result, the surface structure has a delicate effect on water molecules layering. This indicates the hydrophobicity of the borophene surface which was studied in the previous sections. The small out-of-plane height of the borophene structure, which is introduced in sections “Introduction–Results and discussions”, prevents water molecules from being placed between the grooves in the structure. As a result, water molecules are inevitably placed at higher heights. The large distance of the water droplet from the surface of borophene is clearly evident in the graphs of Figs. 6 and 8.

Conclusion

Borophene is a new two-dimensional material that can be a more suitable candidate for industrial applications due to its unique characteristics compared to other two-dimensional materials. For better and more appropriate use of borophene in the industry, it is necessary to have a deep understanding of the surface properties. For this purpose, in the present work, the wettability of borophene surface has been studied with the helping hand of molecular dynamics simulation. Our results show that the borophene surface is hydrophobic and anisotropic. For this reason, two different contact angles have been calculated along the zigzag and armchair directions, which are 149.11° and 148.4°, respectively. The analysis of the energy contour map of the borophene surface showed that there is an energy barrier along the zigzag direction, and this has prevented water molecules from spreading along this direction. The density profile graph showed that the water layers on the borophene surface are not separate from each other and peaks overlap each other overlap. Besides them, the RDF diagram showed that the water drop is located at a great distance from the surface. Finally, the calculation of the number of hydrogen bonds between water molecules revealed that the water drop on the borophene surface maximizes its hydrogen bonding potential. To sum up, all the parameters that studies showed that the borophene surface is hydrophobic.

Author contributions

M.F. reviewed manuscript and analyzed data M.S. wrote the main manuscripts, prepared figures B.M. reviewed figures.

Data availability

Data is provided within the manuscript.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Piazza ZA Hu HS Li WL Zhao YF Li J Wang LS Planar hexagonal B36 as a potential basis for extended single-atom layer boron sheets Nat. Commun. 2014 5 1 6 10.1038/ncomms4113
Piazza, Z. A. et al. Planar hexagonal B36 as a potential basis for extended single-atom layer boron sheets. Nat. Commun. 5, 1–6 (2014).10.1038/ncomms4113
2. Tang H Ismail-Beigi S Novel precursors for boron nanotubes: The competition of two-center and three-center bonding in boron sheets Phys. Rev. Lett. 2007 99 115501 10.1103/PhysRevLett.99.115501 17930448
Tang, H. & Ismail-Beigi, S. Novel precursors for boron nanotubes: The competition of two-center and three-center bonding in boron sheets. Phys. Rev. Lett. 99, 115501 (2007).17930448 10.1103/PhysRevLett.99.115501
3. Tang H Ismail-Beigi S Self-doping in boron sheets from first principles: A route to structural design of metal boride nanostructures Phys. Rev. 2009 80 134113 10.1103/PhysRevB.80.134113
Tang, H. & Ismail-Beigi, S. Self-doping in boron sheets from first principles: A route to structural design of metal boride nanostructures. Phys. Rev. 80, 134113 (2009).10.1103/PhysRevB.80.134113
4. Gindulytė A Lipscomb WN Massa L Proposed boron nanotubes Inorg. Chem. 1998 37 6544 6545 10.1021/ic980559o 11670779
Gindulytė, A., Lipscomb, W. N. & Massa, L. Proposed boron nanotubes. Inorg. Chem. 37, 6544–6545 (1998).11670779 10.1021/ic980559o
5. Mannix AJ Zhou XF Kiraly B Wood JD Alducin D Synthesis of borophenes: Anisotropic, two-dimensional boron polymorphs Science 2015 2015 31513 31516
Mannix, A. J. et al. Synthesis of borophenes: Anisotropic, two-dimensional boron polymorphs. Science 2015, 31513–31516 (2015).
6. Kunstmann J Quandt A Broad boron sheets and boron nanotubes: An ab initio study of structural, electronic, and mechanical properties Phys. Rev. B 2006 74 035413 10.1103/PhysRevB.74.035413
Kunstmann, J. & Quandt, A. Broad boron sheets and boron nanotubes: An ab initio study of structural, electronic, and mechanical properties. Phys. Rev. B 74, 035413 (2006).10.1103/PhysRevB.74.035413
7. Lau KC Pandey R Stability and electronic properties of atomistically-engineered 2D boron sheets J. Phys. Chem. C 2007 111 2906 2912 10.1021/jp066719w
Lau, K. C. & Pandey, R. Stability and electronic properties of atomistically-engineered 2D boron sheets. J. Phys. Chem. C 111, 2906–2912 (2007).10.1021/jp066719w
8. Szwacki NG Matsuda I A historical review of theoretical boron allotropes in various dimensions 2D Boron Boraphene, Borophene, Boronene 2021 2021 1 25
Szwacki, N. G. & Matsuda, I. A historical review of theoretical boron allotropes in various dimensions. 2D Boron Boraphene, Borophene, Boronene 2021, 1–25 (2021).
9. Wu X Dai J Zhao Y Zhuo Z Yang J Two-dimensional boron monolayer sheets ACS Nano 2012 6 7443 7453 10.1021/nn302696v 22816319
Wu, X. et al. Two-dimensional boron monolayer sheets. ACS Nano 6, 7443–7453 (2012).22816319 10.1021/nn302696v
10. Zubarev DY Boldyrev AI Comprehensive analysis of chemical bonding in boron clusters J. Comput. Chem. 2007 28 251 268 10.1002/jcc.20518 17111395
Zubarev, D. Y. & Boldyrev, A. I. Comprehensive analysis of chemical bonding in boron clusters. J. Comput. Chem. 28, 251–268 (2007).17111395 10.1002/jcc.20518
11. Mannix AJ Zhang Z Guisinger NP Yakobson BI Hersam MC Borophene as a prototype for synthetic 2D materials development Nat. Nanotechnol. 2018 13 444 450 10.1038/s41565-018-0157-4 29875501
Mannix, A. J., Zhang, Z., Guisinger, N. P., Yakobson, B. I. & Hersam, M. C. Borophene as a prototype for synthetic 2D materials development. Nat. Nanotechnol. 13, 444–450 (2018).29875501 10.1038/s41565-018-0157-4
12. Feng B Zhang J Zhong Q Li W Li S Experimental realization of two-dimensional boron sheets Nat. Chem. 2016 8 563 568 10.1038/nchem.2491 27219700
Feng, B. et al. Experimental realization of two-dimensional boron sheets. Nat. Chem. 8, 563–568 (2016).27219700 10.1038/nchem.2491
13. Liu X Wang L Li S Rahn MS Yakobson BI Geometric imaging of borophene polymorphs with functionalized probes Nat. Commun. 2019 10 1 7 30602773
Liu, X. et al. Geometric imaging of borophene polymorphs with functionalized probes. Nat. Commun. 10, 1–7 (2019).30602773
14. Wu R Drozdov IK Eltinge S Zahl P Ismail-Beigi S Large-area single-crystal sheets of borophene on Cu (111) surfaces Nat. Nanotechnol. 2019 14 44 49 10.1038/s41565-018-0317-6 30510278
Wu, R. et al. Large-area single-crystal sheets of borophene on Cu (111) surfaces. Nat. Nanotechnol. 14, 44–49 (2019).30510278 10.1038/s41565-018-0317-6
15. Li W Kong L Chen C Gou J Sheng S Experimental realization of honeycomb borophene Sci. Bull. 2018 63 282 286 10.1016/j.scib.2018.02.006
Li, W. et al. Experimental realization of honeycomb borophene. Sci. Bull. 63, 282–286 (2018).10.1016/j.scib.2018.02.006
16. Kiraly B Liu X Wang L Zhang Z Mannix AJ Borophene synthesis on Au (111) ACS Nano 2019 13 3816 3822 10.1021/acsnano.8b09339 30844248
Kiraly, B. et al. Borophene synthesis on Au (111). ACS Nano 13, 3816–3822 (2019).30844248 10.1021/acsnano.8b09339
17. Vinogradov NA Lyalin A Taketsugu T Vinogradov AS Preobrajenski A Single-phase borophene on Ir (111): Formation, structure, and decoupling from the support ACS Nano 2019 13 14511 14518 10.1021/acsnano.9b08296 31790188
Vinogradov, N. A., Lyalin, A., Taketsugu, T., Vinogradov, A. S. & Preobrajenski, A. Single-phase borophene on Ir (111): Formation, structure, and decoupling from the support. ACS Nano 13, 14511–14518 (2019).31790188 10.1021/acsnano.9b08296
18. Zhang Z Yang Y Penev ES Yakobson BI Elasticity, Flexibility, and ideal strength of borophenes Adv. Funct. Mater. 2017 27 1605059 10.1002/adfm.201605059
Zhang, Z., Yang, Y., Penev, E. S. & Yakobson, B. I. Elasticity, Flexibility, and ideal strength of borophenes. Adv. Funct. Mater. 27, 1605059 (2017).10.1002/adfm.201605059
19. Liang T Zhang P Yuan P Zhai S Yang D A molecular dynamics study on the thermal conductivities of single-and multi-layer two-dimensional borophene Nano Futures 2019 3 015001 10.1088/2399-1984/aafc8c
Liang, T., Zhang, P., Yuan, P., Zhai, S. & Yang, D. A molecular dynamics study on the thermal conductivities of single-and multi-layer two-dimensional borophene. Nano Futures 3, 015001 (2019).10.1088/2399-1984/aafc8c
20. Le MQ Mortazavi B Rabczuk T Mechanical properties of borophene films: A reactive molecular dynamics investigation Nanotechnology 2016 27 445709 10.1088/0957-4484/27/44/445709 27678335
Le, M. Q., Mortazavi, B. & Rabczuk, T. Mechanical properties of borophene films: A reactive molecular dynamics investigation. Nanotechnology 27, 445709 (2016).27678335 10.1088/0957-4484/27/44/445709
21. Adamska L Sadasivam S Foley JJ IV Darancet P Sharifzadeh S First-principles investigation of borophene as a monolayer transparent conductor J. Phys. Chem. C 2018 122 4037 4045 10.1021/acs.jpcc.7b10197
Adamska, L., Sadasivam, S., Foley, J. J. IV., Darancet, P. & Sharifzadeh, S. First-principles investigation of borophene as a monolayer transparent conductor. J. Phys. Chem. C 122, 4037–4045 (2018).10.1021/acs.jpcc.7b10197
22. Sun H Li Q Wan XG First-principles study of thermal properties of borophene Phys. Chem. Chem. Phys. 2016 18 14927 14932 10.1039/C6CP02029A 27188523
Sun, H., Li, Q. & Wan, X. G. First-principles study of thermal properties of borophene. Phys. Chem. Chem. Phys. 18, 14927–14932 (2016).27188523 10.1039/C6CP02029A
23. Li D He J Ding G Tang QQ Ying Y Stretch-driven increase in ultrahigh thermal conductance of hydrogenated borophene and dimensionality crossover in phonon transmission Adv. Funct. Mater. 2018 28 1801685 10.1002/adfm.201801685
Li, D. et al. Stretch-driven increase in ultrahigh thermal conductance of hydrogenated borophene and dimensionality crossover in phonon transmission. Adv. Funct. Mater. 28, 1801685 (2018).10.1002/adfm.201801685
24. Carrete J Li W Lindsay L Broido DA Gallego LJ Physically founded phonon dispersions of few-layer materials and the case of borophene Mater. Res. Lett. 2016 4 204 211 10.1080/21663831.2016.1174163
Carrete, J. et al. Physically founded phonon dispersions of few-layer materials and the case of borophene. Mater. Res. Lett. 4, 204–211 (2016).10.1080/21663831.2016.1174163
25. Farahani A Foroutan M Jamshidi M Sababkar M Boudaghi A Introducing structures from hexagonal borophene to nitrophene and their thermal conductivity investigation using a reactive molecular dynamics simulation Langmuir 2023 39 14015 14024 10.1021/acs.langmuir.3c01789 37735112
Farahani, A., Foroutan, M., Jamshidi, M., Sababkar, M. & Boudaghi, A. Introducing structures from hexagonal borophene to nitrophene and their thermal conductivity investigation using a reactive molecular dynamics simulation. Langmuir 39, 14015–14024 (2023).37735112 10.1021/acs.langmuir.3c01789
26. Duo Y Xie Z Wang L Abbasi NM Yang T Borophene-based biomedical applications: Status and future challenges Coord. Chem. Rev. 2021 427 213549 10.1016/j.ccr.2020.213549
Duo, Y. et al. Borophene-based biomedical applications: Status and future challenges. Coord. Chem. Rev. 427, 213549 (2021).10.1016/j.ccr.2020.213549
27. Yadav S Sadique MA Kaushik A Ranjan P Khan R Borophene as an emerging 2D flatland for biomedical applications: Current challenges and future prospects J. Mater. Chem. B 2022 10 1146 1175 10.1039/D1TB02277F 35107476
Yadav, S. et al. Borophene as an emerging 2D flatland for biomedical applications: Current challenges and future prospects. J. Mater. Chem. B 10, 1146–1175 (2022).35107476 10.1039/D1TB02277F
28. Wang X Wu R Xu T Gao Y Mechanical and electrical properties of borophene and its band structure modulation via strain and electric fields: A first-principles study Mater. Res. Express 2021 8 065003 10.1088/2053-1591/ac04ec
Wang, X., Wu, R., Xu, T. & Gao, Y. Mechanical and electrical properties of borophene and its band structure modulation via strain and electric fields: A first-principles study. Mater. Res. Express 8, 065003 (2021).10.1088/2053-1591/ac04ec
29. Jiang HR Shyy W Liu M Ren YX Zhao TS Borophene and defective borophene as potential anchoring materials for lithium–sulfur batteries: A first-principles study J. Mater. Chem. A 2018 6 2107 2114 10.1039/C7TA09244J
Jiang, H. R., Shyy, W., Liu, M., Ren, Y. X. & Zhao, T. S. Borophene and defective borophene as potential anchoring materials for lithium–sulfur batteries: A first-principles study. J. Mater. Chem. A 6, 2107–2114 (2018).10.1039/C7TA09244J
30. Zhang Y Wu ZF Gao PF Zhang SL Wen YH Could borophene be used as a promising anode material for high-performance lithium ion battery? ACS Appl. Mater. Interfaces 2016 8 22175 22181 10.1021/acsami.6b05747 27487298
Zhang, Y., Wu, Z. F., Gao, P. F., Zhang, S. L. & Wen, Y. H. Could borophene be used as a promising anode material for high-performance lithium ion battery?. ACS Appl. Mater. Interfaces 8, 22175–22181 (2016).27487298 10.1021/acsami.6b05747
31. Grixti S Mukherjee S Singh CV Two-dimensional boron as an impressive lithium-sulphur battery cathode material Energy Storage Mater. 2018 13 80 87 10.1016/j.ensm.2017.12.024
Grixti, S., Mukherjee, S. & Singh, C. V. Two-dimensional boron as an impressive lithium-sulphur battery cathode material. Energy Storage Mater. 13, 80–87 (2018).10.1016/j.ensm.2017.12.024
32. Jiang HR Lu Z Wu MC Ciucci F Zhao TS Borophene: A promising anode material offering high specific capacity and high rate capability for lithium-ion batteries Nano Energy 2016 23 97 104 10.1016/j.nanoen.2016.03.013
Jiang, H. R., Lu, Z., Wu, M. C., Ciucci, F. & Zhao, T. S. Borophene: A promising anode material offering high specific capacity and high rate capability for lithium-ion batteries. Nano Energy 23, 97–104 (2016).10.1016/j.nanoen.2016.03.013
33. Shukla V Warna J Jena NK Grigoriev A Ahuja R Toward the realization of 2D borophene based gas sensor J. Phys. Chem. C 2017 121 26869 26876 10.1021/acs.jpcc.7b09552
Shukla, V., Warna, J., Jena, N. K., Grigoriev, A. & Ahuja, R. Toward the realization of 2D borophene based gas sensor. J. Phys. Chem. C 121, 26869–26876 (2017).10.1021/acs.jpcc.7b09552
34. Tian Y Yang H Li J Hu S Cao S A comprehensive first-principle study of borophene-based nano gas sensor with gold electrodes Front. Phys. 2022 17 1 7 10.1007/s11467-021-1094-5
Tian, Y. et al. A comprehensive first-principle study of borophene-based nano gas sensor with gold electrodes. Front. Phys. 17, 1–7 (2022).10.1007/s11467-021-1094-5
35. Hou C Tai G Liu Y Liu X Borophene gas sensor Nano Res. 2022 15 2537 2544 10.1007/s12274-021-3926-6
Hou, C., Tai, G., Liu, Y. & Liu, X. Borophene gas sensor. Nano Res. 15, 2537–2544 (2022).10.1007/s12274-021-3926-6
36. Blossey R Self-cleaning surfaces—virtual realities Nat. Mater. 2003 2 301 306 10.1038/nmat856 12728235
Blossey, R. Self-cleaning surfaces—virtual realities. Nat. Mater. 2, 301–306 (2003).12728235 10.1038/nmat856
37. Miljkovic N Enright R Nam Y Lopez K Dou N Jumping-droplet-enhanced condensation on scalable superhydrophobic nanostructured surfaces Nano Lett. 2013 13 179 187 10.1021/nl303835d 23190055
Miljkovic, N. et al. Jumping-droplet-enhanced condensation on scalable superhydrophobic nanostructured surfaces. Nano Lett. 13, 179–187 (2013).23190055 10.1021/nl303835d
38. Guo P Zheng Y Wen M Song C Lin Y Icephobic/anti-icing properties of micro/nanostructured surfaces Adv. Mater. 2012 24 2642 2648 10.1002/adma.201104412 22488894
Guo, P. et al. Icephobic/anti-icing properties of micro/nanostructured surfaces. Adv. Mater. 24, 2642–2648 (2012).22488894 10.1002/adma.201104412
39. Xu QF Wang JN A superhydrophobic coating on aluminium foil with an anti-corrosive property N. J. Chem. 2009 33 734 738 10.1039/B817130K
Xu, Q. F. & Wang, J. N. A superhydrophobic coating on aluminium foil with an anti-corrosive property. N. J. Chem. 33, 734–738 (2009).10.1039/B817130K
40. Cho KL Liaw II Wu AHF Lamb RN Influence of roughness on a transparent superhydrophobic coating J. Phys. Chem. C 2010 114 11228 11233 10.1021/jp103479k
Cho, K. L., Liaw, I. I., Wu, A. H. F. & Lamb, R. N. Influence of roughness on a transparent superhydrophobic coating. J. Phys. Chem. C 114, 11228–11233 (2010).10.1021/jp103479k
41. Boreyko JB Chen CH Self-propelled dropwise condensate on superhydrophobic surfaces Phys. Rev. Lett. 2009 103 184501 10.1103/PhysRevLett.103.184501 19905808
Boreyko, J. B. & Chen, C. H. Self-propelled dropwise condensate on superhydrophobic surfaces. Phys. Rev. Lett. 103, 184501 (2009).19905808 10.1103/PhysRevLett.103.184501
42. Miwa M Nakajima A Fujishima A Hashimoto K Watanabe T Effects of the surface roughness on sliding angles of water droplets on superhydrophobic surfaces Langmuir 2000 16 5754 5760 10.1021/la991660o
Miwa, M., Nakajima, A., Fujishima, A., Hashimoto, K. & Watanabe, T. Effects of the surface roughness on sliding angles of water droplets on superhydrophobic surfaces. Langmuir 16, 5754–5760 (2000).10.1021/la991660o
43. Boreyko JB Chen CH Vapor chambers with jumping-drop liquid return from superhydrophobic condensers Int. J. Heat Mass Transfer 2013 61 409 418 10.1016/j.ijheatmasstransfer.2013.01.077
Boreyko, J. B. & Chen, C. H. Vapor chambers with jumping-drop liquid return from superhydrophobic condensers. Int. J. Heat Mass Transfer 61, 409–418 (2013).10.1016/j.ijheatmasstransfer.2013.01.077
44. Boreyko JB Zhao Y Chen CH Planar jumping-drop thermal diodes Appl. Phys. Lett. 2011 99 234105 10.1063/1.3666818
Boreyko, J. B., Zhao, Y. & Chen, C. H. Planar jumping-drop thermal diodes. Appl. Phys. Lett. 99, 234105 (2011).10.1063/1.3666818
45. Zhang J Meng Z Liu J Schlaich C Yu Z Breath figure lithography for the construction of a hierarchical structure in sponges and their applications to oil/water separation J. Mater. Chem. A 2017 5 16369 16375 10.1039/C7TA02751F
Zhang, J. et al. Breath figure lithography for the construction of a hierarchical structure in sponges and their applications to oil/water separation. J. Mater. Chem. A 5, 16369–16375 (2017).10.1039/C7TA02751F
46. Wen N Peng S Yang X Long M Deng W A cycle-etching approach toward the fabrication of superamphiphobic stainless steel surfaces with excellent anticorrosion properties Adv. Eng. Mater. 2017 19 1600879 10.1002/adem.201600879
Wen, N. et al. A cycle-etching approach toward the fabrication of superamphiphobic stainless steel surfaces with excellent anticorrosion properties. Adv. Eng. Mater. 19, 1600879 (2017).10.1002/adem.201600879
47. Xia B Liu H Fan Y Zhu W Geng C Preparation of robust CuO/TiO2 superamphiphobic steel surface through chemical deposition and sol–gel methods Adv. Eng. Mater. 2017 19 1600572 10.1002/adem.201600572
Xia, B., Liu, H., Fan, Y., Zhu, W. & Geng, C. Preparation of robust CuO/TiO2 superamphiphobic steel surface through chemical deposition and sol–gel methods. Adv. Eng. Mater. 19, 1600572 (2017).10.1002/adem.201600572
48. Bernardin JD Mudawar I Walsh CB Franses EI Contact angle temperature dependence for water droplets on practical aluminum surfaces Int. J. Heat Mass Transfer 1997 40 1017 1033 10.1016/0017-9310(96)00184-6
Bernardin, J. D., Mudawar, I., Walsh, C. B. & Franses, E. I. Contact angle temperature dependence for water droplets on practical aluminum surfaces. Int. J. Heat Mass Transfer 40, 1017–1033 (1997).10.1016/0017-9310(96)00184-6
49. De Coninck J D'Ortona U Koplik J Banavar JR Terraced spreading of chain molecules via molecular dynamics Phys. Rev. Lett. 1995 74 928 10.1103/PhysRevLett.74.928 10058884
De Coninck, J., D’Ortona, U., Koplik, J. & Banavar, J. R. Terraced spreading of chain molecules via molecular dynamics. Phys. Rev. Lett. 74, 928 (1995).10058884 10.1103/PhysRevLett.74.928
50. D'Ortona U De Coninck J Koplik J Banavar JR Terraced spreading mechanisms for chain molecules Phys. Rev. E 1996 53 562 10.1103/PhysRevE.53.562
D’Ortona, U., De Coninck, J., Koplik, J. & Banavar, J. R. Terraced spreading mechanisms for chain molecules. Phys. Rev. E 53, 562 (1996).10.1103/PhysRevE.53.562
51. Cazabat AM Valignat MP Villette S De Coninck J Louche F The mechanism of spreading: A microscopic description Langmuir 1997 13 4754 4757 10.1021/la970242l
Cazabat, A. M., Valignat, M. P., Villette, S., De Coninck, J. & Louche, F. The mechanism of spreading: A microscopic description. Langmuir 13, 4754–4757 (1997).10.1021/la970242l
52. Rahnamoun A Kaymak MC Manathunga M Götz AW van Duin ACT ReaxFF/AMBER—a framework for hybrid reactive/nonreactive force field molecular dynamics simulations J. Chem. Theory Comput. 2020 16 7645 7654 10.1021/acs.jctc.0c00874 33141581
Rahnamoun, A. et al. ReaxFF/AMBER—a framework for hybrid reactive/nonreactive force field molecular dynamics simulations. J. Chem. Theory Comput. 16, 7645–7654 (2020).33141581 10.1021/acs.jctc.0c00874
53. Mueller JE Van Duin ACT Goddard WA III Development and validation of ReaxFF reactive force field for hydrocarbon chemistry catalyzed by nickel J. Phys. Chem. C 2010 114 4939 4949 10.1021/jp9035056
Mueller, J. E., Van Duin, A. C. T. & Goddard, W. A. III. Development and validation of ReaxFF reactive force field for hydrocarbon chemistry catalyzed by nickel. J. Phys. Chem. C 114, 4939–4949 (2010).10.1021/jp9035056
54. Verma A Zhang W Van Duin ACT ReaxFF reactive molecular dynamics simulations to study the interfacial dynamics between defective h-BN nanosheets and water nanodroplets Phys. Chem. Chem. Phys. 2021 23 10822 10834 10.1039/D1CP00546D 33908500
Verma, A., Zhang, W. & Van Duin, A. C. T. ReaxFF reactive molecular dynamics simulations to study the interfacial dynamics between defective h-BN nanosheets and water nanodroplets. Phys. Chem. Chem. Phys. 23, 10822–10834 (2021).33908500 10.1039/D1CP00546D
55. Plimpton S Fast parallel algorithms for short-range molecular dynamics J. Comput. Phys. 1995 117 1 19 10.1006/jcph.1995.1039
Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).10.1006/jcph.1995.1039
56. Humphrey W Dalke A Schulten K VMD: Visual molecular dynamics J. Mol. Graph. 1996 14 33 38 10.1016/0263-7855(96)00018-5 8744570
Humphrey, W., Dalke, A. & Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 14, 33–38 (1996).8744570 10.1016/0263-7855(96)00018-5
57. Zhou YP Jiang JW Molecular dynamics simulations for mechanical properties of borophene: Parameterization of valence force field model and Stillinger-Weber potential Sci. Rep. 2017 7 1 12 28127051
Zhou, Y. P. & Jiang, J. W. Molecular dynamics simulations for mechanical properties of borophene: Parameterization of valence force field model and Stillinger-Weber potential. Sci. Rep. 7, 1–12 (2017).28127051
58. Young T III An essay on the cohesion of fluids Philos. Trans. R. Soc. Lond. 1805 31 65 87
Young, T. III. An essay on the cohesion of fluids. Philos. Trans. R. Soc. Lond. 31, 65–87 (1805).
59. Foroutan M Torabi-Rad M Boudaghi A Ataeizadeh H The shape of two-dimensional and three-dimensional drops on flat and curved hydrophilic substrates: Variational, numerical and molecular dynamics simulation investigations J. Iran. Chem. Soc. 2021 19 1 11
Foroutan, M., Torabi-Rad, M., Boudaghi, A. & Ataeizadeh, H. The shape of two-dimensional and three-dimensional drops on flat and curved hydrophilic substrates: Variational, numerical and molecular dynamics simulation investigations. J. Iran. Chem. Soc. 19, 1–11 (2021).
60. Gu Z Duan M Tu Y Investigations on different two-dimensional materials as slit membranes for enhanced desalination J. Membr. Sci. 2020 523 115452
Gu, Z., Duan, M. & Tu, Y. Investigations on different two-dimensional materials as slit membranes for enhanced desalination. J. Membr. Sci. 523, 115452 (2020).
61. Foroutan M Bavani BM Boudaghi A Controlled hydrophilization of black phosphorene: A reactive molecular dynamics simulation approach Phys. Chem. Chem. Phys. 2022 24 27532 27547 10.1039/D2CP02382B 36345729
Foroutan, M., Bavani, B. M. & Boudaghi, A. Controlled hydrophilization of black phosphorene: A reactive molecular dynamics simulation approach. Phys. Chem. Chem. Phys. 24, 27532–27547 (2022).36345729 10.1039/D2CP02382B
62. Torabi-Rad M Foroutan M Wettability of penta-graphene: A molecular dynamics simulation approach J. Phys. Chem. C 2022 126 1590 1599 10.1021/acs.jpcc.1c06826
Torabi-Rad, M. & Foroutan, M. Wettability of penta-graphene: A molecular dynamics simulation approach. J. Phys. Chem. C 126, 1590–1599 (2022).10.1021/acs.jpcc.1c06826
63. Foroutan M Darvishi M Fatemi SM Structural and dynamical characterization of water on the Au (100) and graphene surfaces: A molecular dynamics simulation approach Phys. Rev. E 2017 96 033312 10.1103/PhysRevE.96.033312 29346917
Foroutan, M., Darvishi, M. & Fatemi, S. M. Structural and dynamical characterization of water on the Au (100) and graphene surfaces: A molecular dynamics simulation approach. Phys. Rev. E 96, 033312 (2017).29346917 10.1103/PhysRevE.96.033312
