Carbon Monosulfide Nanostructures: Chains Arrays, Monolayers, and Thin Films

Theoretical predictions have lead to the experimental synthesis of new low dimensional layered structures. Herein we show for the very first time that compounds of carbon monosulfide exhibit a great variety of layered nanostructures, such as chains arrays, monolayers, and thin films. We find that the chains arrays are the most stable because they are mainly dimensionality-driven by the $sp^2$ hybridization of sulfur and carbon orbitals. Furthermore, the chains arrays are direct gap semiconductors. In contrast to thin films, the monolayers are stable at room temperature with a semiconductor phase followed in energy by a metallic phase. Then, we achieve a semiconductor-to-metal phase transition in carbon sulfur monolayers, which can be driven by strain engineering controlling conductivity and carrier mobility.


I. INTRODUCTION
New phenomena in condensed matter physics are based on two-dimensional materials.An example is a single layer of carbon atoms known as graphene 1 .Graphene characterization shows outstanding structural and electronic properties such as high stiffness or flexibility, and massless Dirac fermions with low resistivity [2][3][4][5][6][7] , which allow for many applications.The main drawback is the absence of band gap in the graphene electronic structure and the small on/off ratio 4,8 , a band-gap may be opened via different techniques, such as chemical tailoring, inducing strains, and geometric patterning [9][10][11][12][13] .Although graphene looks promising for the future, research is today focusing on other carbon nanostructured materials and other two-dimensional materials 14 .
The possibility to exfoliate materials and the existence of accurate experimental techniques to study single-layer materials [15][16][17] supposed a big boost to study other layered materials.Layers with other elements in the carbon group have been studied to produce counterparts to graphene, such as silicene, germanene or stanene [18][19][20][21] .Graphene is flat because there is sp 2 hybridization, and these other compounds adopt a buckled structure because they prefer sp 3 hybridization.Other layered nanostructures isoelectronic to graphene are obtained either combining III-V group elements, such as in hexagonal boron nitride 22 .2D materials research is recently widen to other new materials based on phosphorus, whose isoelectronic counterparts are a central issue nowadays [23][24][25] .Black phosphorus is formed by single layers mediated by van der Waals interaction, so that analogously to graphene, few layers are isolated to form phosphorene 26,27 .Phosphorene has a two dimensional honeycomb puckered structure where each atom is bonded to three neighbors.It displays a direct band gap, tunable attending to the number of layers and the stacking 28,29 .
More and more systems are being proposed and deeply analyzed to obtain 2D materials with improved stability and new properties 30 ; one of the hottest topic nowadays focuses on phosphorene isoelectronics compounds.
They are currently undergoing characterization, for instance, by studying doping 42 , point defects 42,43 , and oxygen contamination 44 .The group IV monochalcogenides includes a combination of light elements, such as in silicon monosulfide monolayers and thin films [23][24][25] , and a combination of heavy (Ge, Sn) and light (O, S) atoms 45 , such as in silicon telluride 46 and germanium monosulfide 47 monolayers, which exhibit strain-tunable indirect band gap.Nevertheless, none of the today proposed phosphorene isoelectronic compounds is based on carbon.In fact, the layered carbon compounds, such as carbon nitride 48 and carbon phosphide 49,50 , are not isoelectronic to phosphorene.Our aim is to design new layered nanostructures isoelectronic to black phosphorus, but still composed of carbon.
We present here for the first time stable carbon monosulfide nanostructures covering chains, monolayers, and thin films.We have performed phonon and molecular dynamics calculations to ensure the stability of those structures at room temperature.The carbon monosulfide nanostructures vary from metallics in thin films to semiconductors in chains and most monolayers.The chains and thin films are more stable than the monolayers because sulfur atoms have two bonds instead of three (or four) for the monolayers 25 .The large variety of nanostructures that we explored with either metallic or semiconductor character makes this material an interesting compound for several applications based on electronic transport within semiconductors devices.

II. COMPUTATIONAL DETAILS
We use the SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms) package to carry out density functional theory (DFT) calculations for carbon monosulfide nanostructures.For the exchange and correlation potentials we used the generalized gradient approximation (GGA) in the form of Perdew-Burke-Ernzenhof 51 .We fully relax atomic positions and unit cells until forces are below 0.006 eV/ Å.The used unit cells have large vectors about 24.5 Å in the perpendicular direction to the layers in order to avoid interaction with images.An electronic temperature of 25 meV and a meshcutoff of 250 Ry are used.The sampling of the Brillouin zone has a fine grid of 20×20×1 k points.The atomic cores are described by nonlocal norm-conserving Troullier-Martins 52 pseudopotentials factorized in the Kleynman-Bylander form.Orbitals are developed in basis sets on each atom, and the basis size is double zeta plus polarization orbitals.Details on the used pseudopotentials and basis were described previously 9,25,53 .Anyhow our main results are checked by repeating calculations with the VASP code using the projected augmented wave method (PAW) and within the PBE formalism for the exchange and correlation 54,55 , and the Quantum ESPRESSO method with equivalent parameters 56 .We estimate the van der Waals contribution to the total energy for a two dimensional array with well separated chains using the implementation within the VASP code 57,58 of the vdW-DF 59 functional.Last but not least, we even confirm the stability of such highly stable obtained structures by carrying out molecular dynamics calculations within the SIESTA code using the Nose thermostat at room temperature.We employed a Nose mass [60][61][62] of 10 Ry fs 2 , and a time step of 1 fs.We have chosen 3000 as final time step and the relaxation time to reach target temperature was 2500 fs.Phonon dispersion curves are obtained using the Quantum ESPRESSO code, with details included in the Supplemental Information.

III. RESULTS AND DISCUSSION
We study three different types of carbon monosulfide nanostructures, monolayers, thin films, and chains arrays, following the order from lower to higher stability.Molecular dynamics simulations show that thin films are not stable at room temperature.Comments and analysis of thin film structures are therefore presented in the Supplemental Information, together with further information on monolayers and chains.

A. Monolayers
We explore a wide variety of structures for carbon monosulfide monolayers.We next present the geometry and electronic structure of the four most stable structures, and finally consider the ground state in more detail.Figure 1 shows four monolayers of carbon monosulfide after structural relaxation including their relative energies.The most stable α monolayer presents a similar structure to that of a single layer of black phosphorus.This structure has also been reported recently for group IV monochalcogenides 23,31 , and group V semiconductors such as phosphorous-nitrogen or arsenic-nitrogen 40,41,63 .The results obtained with VASP method yield energy differences slightly larger, although the stability order is the same.Second in energy we find the δ monolayer that is a regular structure composed of squares looked from above.A side view reveals two different heights for sulfur atoms.Next there is a similar structure, the κ structure, with all the carbon atoms shifted towards one side producing long and short distances with the sulfur atoms.When the distortion is maximized the α structure is reached.The κ layer seems an intermediate structure between the α and δ monolayers, so that the three monolayer monosulfide-carbon structures seems related, to be discussed at length in next section III A 2. Phonon dispersions, shown in Supplemental Information, indicate that the α and δ structures are really stable, and the κ structure is unstable because there are negative frequencies.Furthermore, molecular dynamics simulations show that the α and δ structures are stable even at room temperature.Last we find the β structure which from a top view it displays a hexagonal pattern, and in a side view it has a buckled structure, reported for other phosphorene isoelectronic compounds 23 .The β structure is similar to that of blue phosphorene 64,65 .The phonon dispersion curves show that the carbon monosulfide β phase is unstable.
To gain insight on the electronic structure of the carbon monosulfide monolayers we consider the band structures shown in Fig. 2. On the one hand, three of the monolayers, labeled as α, κ and β, are semiconductors with indirect band gaps of 0.97, 0.77 and 1.64 eV, respectively.For the α and κ layers, the minimum of the conduction band is located at the gamma point.For the α case, the highest occupied valence band is well separated from lower valence bands with a pseudogap of 2 eV.The κ structure has still a pseudogap opening in the valence band, although the value is reduced to 0.5 eV.On the other hand, the δ structure is metallic.
The bands of the α, δ and κ structures appear similar because they share common structural features as commented above.For the β monolayer the bands are degenerated throughout the Brillouin zone due to the lattice symmetry.Most of the features of the bands are the same for the α and κ cases; for instance there are two occupied bands just below the Fermi energy, extended along all Brillouin zone, and well separated from the deeper states.However, the band structure of the δ case differs because one of the two occupied bands close to the Fermi energy becomes broad and metallized, going down to deep energies following a parabolic shape.The two occupied narrow bands for the α and κ cases are localized within a interval of 2 eV, while the bands for the δ case become broad expanding over more than 6 eV.We find two unoccupied conduction bands that mirror the behavior of two nearest valence bands with respect to the Fermi energy.There are even bands for the δ structure crossing the Fermi energy, resulting in a conducting monolayer.
We furthermore study the charge transfer between carbon and sulfur atoms using the Mulliken scheme.The charge values are 0.81, 0.92, 0.65 and 0.66 electrons for the α, δ, κ and β monolayers, respectively.We find that the charge transfer of almost one electron in some cases is large and goes from carbon to sulfur atoms.This result is surprising because carbon and sulfur have almost the same Pauling electronegativity.It seems that the larger size of sulfur atoms with respect to carbon atoms must be taken into account in order to explain the charge gain by sulfur atoms.

α monolayer
We focus on the ground state monolayer given by the α structure.Each carbon atom is bonded to three sulfur atoms and vice versa.There are two bonds of 1.82 Å and one of 1.90 Å.The angles between the bonds from a carbon atom are 113 • in two cases and 103 • in the other.It seems that carbon atoms present a sp 3 -like hybridization.
We discuss on the relative contribution between carbon and sulfur atoms to discuss stability.The projected density of states is shown in Fig. 3(a).The total states projected on sulfur contribute more than those on carbon because sulfur has a larger number of electrons than carbon.Nevertheless, the occupied peak close to the Fermi energy has a larger carbon contribution than sulfur.This particular peak is related to the isolated valence band, discussed above when commenting Fig. 2. The density of states (LDOS) in real space within the energy range of 2 eV is included in panel (b).Most of the states are localized around carbon atoms with orbital lobes suggesting a large p-type contribution.The peak LDOS globally looks like a dangling bond of the sp 3 hybridization.However, the wave function for the highest occupied orbital in Γ point, as shown in Fig. 3c, reveals a p z -like orbital over the carbon atoms, to be held responsible for the broadening of a band close to Fermi energy upon distortions.The projected sulfur LDOS shown in Fig. 3(a) has a non negligible d contribution, especially for the highest energy occupied peak.Therefore, spd hybridization on sulfur atoms seems crucial to describe the valence band close to the Fermi energy 25 .
We then compare carbon-sulfur and silicon-sulfur 23,25 compounds, both being isoelectronic to phosphorene, in the α structure.There are two well differentiated roles depending of the atom size, as shown in the α structure of Fig. 1.The smaller atom, either carbon for CS or sulfur for SiS, is the inner one, and the larger atom is the outer one.The electronic structure produces similar trends taking into account that the inner atom contributes with p z orbitals to the highest occupied orbital at the Γ point.It is noteworthy that sulfur atoms can adopt the two roles depending on the counterpart element involved in the isoelectronic compound to phosphorene.

Strain-driven transition between two stable monolayer phases
We try to shed light over the striking metallization of the δ monolayer.A distortion of the α atomic positions could end in the δ structure, so that the monolayer could undergo a semiconductor-metal transition driven by a strain field.The δ structure is isotropic and similar to a rippled square lattice.Sulfur atoms show spd hybridization as denoted by the projected density-of-states, similarly to the α structure.A detailed analysis of the localized density of states shows that the Fermi energy states are p-type on carbon, which could be seen as a square lattice of p z orbitals.Details on the geometry and the projected and localized density-of-states of the δ phase are found in the Supplemental Information.
We follow the transition between the α and δ monolayers by changing aspect ratio, as shown schematically in Fig. 4(a).To this aim we modify the unit cell shape step-by-step and transform the square cell into a rectangular one.The unit cell following the transformation is characterized by the aspect ratio between length and width, which is 1 for the δ phase and 1.4 for the α phase.For each calculation the unit cell remains fixed and the inner atomic positions are relaxed.We expected that the δ monolayer changing aspect ratio would transform into the α monolayer by relaxing the inner atomic positions.However, we find that the inner atomic structure is retained when deforming the unit cell, and the energy increases continuously with the aspect ratio, as shown in Fig. 4(b).We also find a similar behavior when deforming the α monolayer.The two curves of energy versus aspect ratio cross at a point with a barrier of about 0.2 eV between the monolayers, a fact sufficient to confirm the stability for the δ structure.The energy of the κ monolayer, included in Fig. 4(b), is very close to the crossing point, which supports the idea of κ being a metastable phase between α and δ monolayers.Figure 4(c) sketches a possible energy profile following this picture.

B. Chains arrays
We find other kind of structures for carbon monosulfide with chains arrays.The previously studied γ SiS monolayer 25 is fully relaxed by substituting silicon with carbon atoms.A layer with chains of carbon monosulfide is then stabilized by more than half electronvolt per atom when compared to the most stable α monolayer.We obtained two different chains arrays with hexagonal and pentagonal motifs.The hexa chains are more stable than the penta chains by more than 0.25 eV per atom, so we focus on the hexa chains below.The penta chains are commented in the Supplemental Information.The structure with hexa-motif chains, shown in Fig. 5(a), is 0.766 eV per atom more stable than the ground state α monolayer, which results in a large rise of stability.The distance between the closest atoms from two different chains is larger than 3.7 Å, which means that interaction between adjacent chains is mediated by dispersion forces.Thus, the bonds between chains are much weaker than those established within a single chain. 69e investigate in more detail the bonds between atoms in the chains.The top view in Fig. 5 (a) shows that due to kinks, the hexa chains combines long and short stretches.The kinked chains are more stable by an 70 meV per atom than the straight chains.We check the structural stability of this structure with chains using molecular dynamics simulation for a single chain at room temperature.Carbon atoms form covalent bonds with sulfur atoms.Each carbon atom is bonded to two sulfur atoms with two bond lengths of about 1.82 Å and to another carbon atom with a bond length of 1.38 Å.The short bond length between carbon atoms indicates the existence of a double bond with a σ bond and a π bond.The π bond is spatially shown by the highest occupied molecular orbital at the Γ point in Fig. 5(b), formed by two p-type orbitals from each pair of carbon atoms.Concerning the angles between the three bonds, there are two different carbon atoms depending on whether they are in a short stretch or in a long stretch; however, for all the cases the angles are close to 120 • .The fact that the three bonds are planar and form angles close to 120 • is assigned to a sp 2 hybridization.The remaining p z orbital is used to form the double π bond with the neighboring carbon atom.Focusing now on the sulfur atoms, each atom establishes two bonds with carbon.The chain adopts a zigzag pattern with two consecutive sulfur atoms in kinks followed by other in a straight stretch.The two types of sulfur atoms seem to have the same hybridization.They contribute with p-orbitals to the highest occupied orbital at the Γ point, as shown in Fig. 5(b).However, the kinks effect is to tilt the p z orbitals by an angle θ of approximately 25 • , as presented in Fig. 5(b).
The analysis of the projected density of states for the hexa chains array, attached in the Supplemental Information, shows that close to the Fermi energy there are mainly p-type contributions, in contrast to the α and δ monolayers.The charge reorganization under the Mulliken scheme results in a net transfer of approximately 0.6 electrons from each carbon atom towards the sulfur atoms.Furthermore, we analyze the band structure shown in Fig. 5(c) along two directions, parallel and perpendicular to the chains.We find that the hexa chains array is semiconductor with a direct band gap of 1.03 eV at the Γ point, a fact that could be interesting for electronic applications.

IV. FINAL REMARKS
Finally, we comment on the stability order between the different layered nanostructures of carbon monosulfide we have studied this far.We find chains, thin films and monolayers in order of increasing energy.The energy is explained by hybridization and the numbers of bonds for carbon and sulfur atoms in each structure, as summarized in Fig. 6.The α monolayer, taken as reference, presents a sp 3 hybridization for carbon atoms with three bonds and a dangling bond, and a spd hybridization for sulfur atoms with three bonds.The thin films have lower energies per atom because the carbon atoms establish a fourth bond and the sulfur atoms lost a bond.Note that sulfur atoms prefer to establish two bonds to fulfill the octet rule 24,25 .The thin films are unstable at room temperature in molecular simulations.Last but not least, the hexa chains array improves energy and stability with respect to thin films.The gain in stability for chains is ascribed to the carbon preference for the sp 2 hybridization with two single bonds and a double bond over the fourfold coordination of the sp 3 hybridization in thin films.
In summary, we present brand new stable nanostructures for carbon monosulfide, which include stable chains arrays and monolayers.On the one hand, the chains arrays, which are more stable than monolayers, present hexagonal and pentagonal patterns with carbon and sul-fur atoms displaying sp 2 hybridizations.These results pave the way so that this type of nanostructures with carbon monosulfide chains can be synthesized as freestanding chains, but probably better grown supported on some substrates as chains arrays.On the other hand, the monolayers show sp 3 and spd hybridizations for carbon and sulfur atoms, respectively.The ground state α monolayer has the structure of a single layer of black phosphorus with an indirect band gap of approximately 1 eV.We find that strain can induce a transition between the most stable α and δ monolayers.Because the α monolayer is semiconductor and the δ monolayer is metallic, the α to δ transition could just allow to control conductivity by applying strain, an effect that looks promising for the design of new electronic devices.We further study the striking metallization of the δ monolayer.A distortion of the α atomic positions could end in the δ structure, meaning that by engineering with a strain field the monolayer could undergo a semiconductor-metal transition.The δ structure is isotropic and similar to a rippled square lattice.All the carbon-sulfur bonds are equal with a distance of 1.93 Å.The consecutive carbon bonds have four angles of 96 • , the opposite bonds have two angles of 141 • .The angles between the sulfur bonds are four of 84 • .The right angles for sulfur atoms are related to d orbitals taken part in spd hybridization, as shown using the density of states in Fig. 7(a).However, the occupied states close to the Fermi energy present a major contribution from the states p of carbon and a small contribution with d-character from sulfur.The spatial localized density of states (LDOS) for two energy intervals is shown in Fig. 7 (b).Around to the Fermi energy, the states in interval 1 are localized on carbon atoms, like a square lattice of p z orbitals explaining metallic behavior.They become metallic along the ΓM direction in reciprocal space, equivalent to the (110) direction in real space among carbon atoms.Along the ΓX direction, there are also empty states very close to the Fermi energy, but with a small gap of 0.22 eV.These states plotted in interval 2 are now including lobules on sulfur towards the outer part of the δ monolayer.

B. Chains arrays
1. Density of states for the hexa chains Figure 8(a) shows the projected density of states resolved per atom site.We find that the s and p orbitals are hybridized at deep energies.Close to the Fermi energy all the states come from the p orbitals centered in carbon and sulfur atoms.The d orbitals for chains arrays give just residual contributions.

Penta chains arrays
Because the carbon-carbon bond seems the cause of stability we next investigate structures in which the C 2 unit has rotated, as in Stone-Wales rotations.Figure 8(b) shows a penta array with chains interacting through their edges.The presence of pentagons motivates the structure name.The distances between carbon and sulfur atoms between chains are larger than 3.5 Å.Therefore, there are few atoms interacting through dispersion forces, which means that van der Waals contribution would be lower than for the hexa motif.Focusing on a single penta chain, there are two remarkable facts that explain the loss of stability with respect to the hexa chains.First, the penta motif is similar to the hexa motif with a carboncarbon double bond rotated 90 • , so that instead of having two hexagons we have the two less stable pentagons.The double bonds are preserved with typical bond lengths between carbon atoms of 1.39 Å in most of the cases, and 1.45 Å close to a terrace.The carbon-sulfur bond length is around 1.8 Å.Second, the individual chains have steps each couple of two pentagons linked by a carbon double bond.The angles of carbon bonds are close to 120 • , a characteristic of sp 2 -like hybridization.The two bonds of sulfur atoms have angles about 95 • , except for the sulfur atoms close to terraces where the angle values ( e.g.99 • and 115 • ) are slightly larger.The charge transfer goes from carbon to sulfur with a value larger than 0.5 electrons.Thus, the penta array with chains presents hybridizations and a bonding scheme similar as those for the hexa motif.It seems that having arrays with chains formed by carbon molecules separated by sulfur atoms is crucial for having high stability with respect to monolayers.

C. Thin films
Carbon group elements form buckled monolayers, such as in silicene and germanene, because sp 3 hybridization is favored rather than sp 2 hybridization.Silicon nanostructures are expected to gain energy when higher silicon coordination is allowed, as reported for silicon monosulfide thin films 24,25 .We study carbon monosulfide following the silicon monosulfide thin films structures.We find that there is a large energy gain with respect to the α monolayer because a fourth bond is achieved for the carbon having sp 3 hybridization.On the one hand, the most stable monolayer do not exhibit a purely carbon sp 2 hybridization, but it has a strong sp 3 character with a dangling bond.On the other hand, sulfur prefers to have two bonds instead of three in order to fulfill the octet rule 24 .Consequently the energy per atom increases notably with respect to monolayers because the carbon dangling bond is passivated and a sulfur bond is lost.Figure 9 shows the detailed geometries of P ma2 and P mma structures.The coordination number, number of bonds, and distances are very similar for both structures.The difference between the P ma2 and P mma structures is the slight distortion around sulfur atoms, which requires considering a large unit cell.The unit cell of P ma2 is large than that of P mma, as shown by the unit cell vectors in Fig. 9.For the P ma2 structure, the sulfur angle between the two C-S bonds is 99 • .Carbon atoms bind to two carbon atoms located at 1.71 Å and two sulfur atoms at 1.86 Å.The angles between bonds established by carbon atoms are all close to 109.5 • , the typical angle for an sp 3 hybridization as for methane.The description for the P mma case reveals almost identical findings.
The band structures of the P ma2 and P mma structures are presented in Fig. 9.The electronic structure of these structures look similar as expected.However, the structural distortion of P ma2 is splitting the degener-ate bands of P mma.The conduction bands close to the Fermi energy expand within a broad range of more than 4 eV.We find that there is no band gap, which means the carbon monosulfide thin films are metallic.The states are metallic in the band crossing the Fermi energy around the Γ point and in the ΓX direction, regions to be held responsible for electronic conduction.
The projected density of states of the P ma2 and P mma structures shown in Fig. 9 for both cases are similar, as expected.For carbon we distinguish three zones: an interval [-15,-10] eV with sp hybridization, an interval [-10,-4] with almost all the p-type contribution, and an interval [-4,0] where the contribution from carbon is very small compared to that of sulfur.For sulfur there is large s-type contribution very deep in energy for the interval [-15,-11], and p-type contribution is largely predominant in energies around the Fermi level.In fact, the sulfur p-type states near the Fermi energy are mostly responsible for conduction.
Molecular dynamics simulations show that the two proposed thin films P ma2 and P mma are not stable at room temperature.The structures break into separated parts, showing preference for chain-like geometries in agreement with our results given above.The absence of a band gap for the P ma2 and P mma structures could be related to the low stability of carbon monosulfide thin films.

D. Phonon dispersion and molecular dynamics
We have carried out molecular dynamics simulations using the Nose-Hoover thermostat in order to study the stability of the carbon monosulfide chains, thin films, and monolayers at room temperature.We show the results for the two most stable monolayers α and δ and for the hexa chains arrays, which are found to be stable at room temperature, in Fig. 11 and Fig. 12.Molecular dynamics simulations of the thin film geometries P ma2 and P mma for carbon monosulfide show that they are not stable at room temperature.
Phonon dispersions are calculated using linear response in the Quantum ESPRESSO 56,66,67 method, a plane wave code that uses ultrasoft pseudopotentials 68 with a cutoff energy of 50 Ry and the PBE GGA.We have also computed the phonon dispersion curves for α and δ carbon monosulfide monolayers, as shown in Fig. 10.There are no negative frequencies.

FIG. 1 :
FIG.1: Geometries of the four carbon monosulfide monolayers in order of decreasing stability.Carbon (sulfur) atoms are represented by black (yellow) spheres through the article.The unit cell is shown in green for each case.Differences in the energy per atom using the SIESTA and VASP methods are shown on top and in parentheses for each structure, respectively.The energy differences are calculated relative to the energy per atom in the α structure.

FIG. 3 :
FIG. 3: Electronic properties of the ground α structure.(a) Density of states projected on carbon (up panel) and sulfur (down panel).(b) Spatial localized density of states for the energy peak marked up with a blue rectangle in (a).(c) Highest occupied orbital of the valence band at the Γ point.

Figure 3 (
FIG. 4: Stretching the α and δ monolayer as shown in panel (a) yields an increasing energy as shown in panel (b), suggesting a barrier between those structures.(c) Relative position in energy of the three most stable monolayers sketching a possible energy profile.

FIG. 5 :
FIG. 5: (a) Geometry of the hexa chains arrays, with the hexagon motif of the structure highlighted in blue color.The unit cell to compute the band structures are marked by vectors.The energies per atom obtained using SIESTA and VASP codes are given above referred to the ground state α monolayer.(b) Highest occupied molecular orbital (HOMO) at the Γ point for the hexa chains in the region marked with a black ellipse in panel (a).(c) Band structure of the hexa chains.

FIG. 6 :
FIG. 6: Comparison of energetic stability, number of bonds, and hybridizations for carbon monosulfide layered nanostructures.Note that the most stable structures are the chains arrays.

FIG. 7 :
FIG. 7: Electronic properties of the δ monolayer.(a) Density of states projected on carbon (top) and sulfur (bottom).(b) Spatial localized density of states for the two energy intervals marked in the band structure around the Fermi energy.

FIG. 8 :
FIG. 8: (a) Density of states projected on carbon (up) and sulfur (down) for the hexa chains array.(b) Geometry and energy of the penta chains array.Energy origin is relative to the α monolayer of Fig. 1.

FIG. 9 :
FIG. 9: Geometries for two carbon-sulfur thin films in the (a) Pma2 and (b) Pmma structures.The unit cells employed for the electronic analysis are shown in green.Energies per atom with the SIESTA and VASP methods relative to the α monolayer are given on top and within parenthesis, respectively.Band structures and projected density of states are shown in the lower part.

FIG. 11 :
FIG. 11: Energy versus time step for the α (a) and δ (b) monolayers in a Nose thermostat at 300K.Structural snapshots are shown below at the four steps marked up in the graph for each case.

TABLE I :
Lattice-vectors, first-neighbor distances, and thickness ∆ z (in Å) for the four carbon monosulfide monolayers.Band gap (in eV).