Source author record

Gang Su

Gang Su appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

64works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

64 published item(s)

preprint2026arXiv

Dominant Kitaev Interaction and Field-induced Quantum Disordered Phase in the Cobaltate Na$_2$Co$_2$TeO$_6$

The identification of quantum spin liquid phases in Kitaev candidate materials remains a major experimental challenge. Since most Kitaev candidates develop antiferromagnetic (AFM) order at low temperatures, currently there are great interest on the field-induced magnetic disordered phase in these compounds, that are distinct from (partially) polarized states. Recently, a cobaltate Na$_2$Co$_2$TeO$_6$ has emerged as a promising Kitaev candidate with high-spin $t^{5}_{2g}e^2_g$ configuration and spin-orbit entangled $J_{\rm eff} = 1/2$ honeycomb lattice system. There are intensive studies on field-induced magnetic states and phase transitions under in-plane magnetic fields. In this study, we propose an intermediate disordered phase induced by an out-of-plane field along the $c$-axis, through high-field magnetization and magnetocaloric effect measurements. To explain the high-field behavior of Na$_2$Co$_2$TeO$_6$, we develop an effective $K$-$J$-$Γ$-$Γ^{\prime}$ spin model featuring a dominant AFM Kitaev interaction. This framework uncovers an intermediate quantum spin liquid phase, establishing the material as a unique platform for exploring Kitaev physics and field-induced quantum-disordered states.

preprint2023arXiv

Dipolar Spin Liquid Ending with Quantum Critical Point in a Gd-based Triangular Magnet

By performing experiment and model studies on a triangular-lattice dipolar magnet KBaGd(BO$_3$)$_2$ (KBGB), we find the highly frustrated magnet with a planar anisotropy hosts a strongly fluctuating dipolar spin liquid (DSL), which originates from the intriguing interplay between dipolar and Heisenberg interactions. The DSL constitutes an extended regime in the field-temperature phase diagram, which gets lowered in temperature as field increases and eventually ends with an unconventional quantum critical point (QCP) at $B_c\simeq 0.75$~T. Based on dipolar Heisenberg model calculations, we identify the DSL as a Berezinskii-Kosterlitz-Thouless (BKT) phase with emergent U(1) symmetry. Due to the tremendous entropy accumulation that can be related to the strong BKT and quantum fluctuations, unprecedented magnetic cooling effects are observed in the DSL regime and particularly near the QCP, making KBGB a superior dipolar coolant to commercial Gd-based refrigerants. We establish the phase diagram for triangular-lattice dipolar quantum magnets where emergent symmetry plays an essential role, and provide a basis and opens an avenue for their applications in sub-Kelvin refrigeration.

preprint2022arXiv

Flat electronic band structure and anisotropic optical, mechanical, and thermoelectric properties of two-dimensional fullerene networks

Nanoclusters like fullerenes as the unit to build intriguing two-dimensional topological structures is of great challenge. Here we propose three bridged fullerene monolayers and comprehensively investigate the novel fullerene monolayer as synthesized experimentally Zheng et al.,[Nature 606, 507-510 (2022)] by state of the art first principles calculations. Our results show that alpha-C60-2D has a direct bandgap of 1.49 eV owing to a flat conduction band bottom close to the experimental value, the optical linear dichroism with strong absorption in long-wave ultraviolet region, a small anisotropic Youngs modulus, the large hole mobility, and the ultrahigh Seebeck coefficient at middle low temperatures. Moreover, Li ions are found to migrate easily along the X path in alpha-C60-2D. It is unveiled that the anisotropic optical, mechanical, electrical, and thermoelectric properties of alpha-C60-2D originate from the asymmetric bridging arrangements between C60 clusters. Our study promises potential applications of monolayer fullerene networks in diverse fields.

preprint2022arXiv

Titanium-based kagome superconductor CsTi_3Bi_5 and topological states

Since the discovery of a new family of vanadium-based kagome superconductor AV3Sb5 (A=K, Rb, and Cs) with topological band structures, extensive effort has been devoted to exploring the origin of superconducting states and the intertwined orders. Meanwhile, searching for new types of superconductors with novel physical properties and higher superconducting transition temperatures has always been a major thread in the history of superconductor research. Here we report a successful fabrication and the topological states of a Titanium-based kagome metal CsTi3Bi5 (CT3B5) crystal. The as-grown CT3B5 crystal is of high quality and possesses a perfect two-dimensional kagome net of Titanium. The superconductivity of the CT3B5 crystal shows that the critical temperature Tc is of ~4.8 K. First-principle calculations predict that the CT3B5 has robust topological surface states, implying that CT3B5 is a Z2 topological kagome superconductor. This finding provides a new type of superconductors and the base for exploring the origin of superconductivity and topological states in kagome superconductors.

preprint2022arXiv

Topological superconductivity and large spin Hall effect in the kagome family Ti6X4 (X = Bi, Sb, Pb, Tl, In)

Topological superconductors (TSC) become a focus of research due to the accompanying Majorana fermions. However, the experimentally reported TSC are extremely rare. The recent experiments reported the kagome TSC AV$_{3}$Sb$_{5}$ (A=K, Rb, Cs), which exhibit unique superconductivity, topological surface states (TSS), and Majorana bound states. More recently, the first titanium-based kagome superconductor CsTi$ _{3} $Bi$ _{5} $ with nontrivial topology was successfully synthesized as a perspective TSC. Given that Cs contributes little to the electronic structures of CsTi$ _{3} $Bi$ _{5} $ and binary compounds may be easier to be synthesized, here, by density functional theory calculations, we predict five stable non-magnetic kagome compounds Ti$ _{6} $X$ _{4} $ (X = Bi, Sb, Pb, Tl, In) which exhibit superconductivity with critical temperature Tc = 3.8$ - $5.1 K, nontrivial $\mathbb{Z}$$_2$ band topology, and TSS close to the Fermi level. In addition, the large intrinsic spin Hall effect is obtained in Ti$ _{6} $X$ _{4} $, which is caused by gapped Dirac nodes and nodal lines due to a strong spin-orbit coupling. This work offers new platforms for TSC and spintronic devices.

preprint2021arXiv

Electric field induced topological phase transition and large enhancements of spin-orbit coupling and Curie temperature in two-dimensional ferromagnetic semiconductors

Tuning topological and magnetic properties of materials by applying an electric field is widely used in spintronics. In this work, we find a topological phase transition from topologically trivial to nontrivial states at an external electric field of about 0.1 V/A in MnBi$_2$Te$_4$ monolayer that is a topologically trivial ferromagnetic semiconductor. It is shown that when electric field increases from 0 to 0.15 V/A, the magnetic anisotropy energy (MAE) increases from about 0.1 to 6.3 meV, and the Curie temperature Tc increases from 13 to about 61 K. The increased MAE mainly comes from the enhanced spin-orbit coupling due to the applied electric field. The enhanced Tc can be understood from the enhanced $p$-$d$ hybridization and decreased energy difference between $p$ orbitals of Te atoms and $d$ orbitals of Mn atoms. Moreover, we propose two novel Janus materials MnBi$_2$Se$_2$Te$_2$ and MnBi$_2$S$_2$Te$_2$ monolayers with different internal electric polarizations, which can realize quantum anomalous Hall effect (QAHE) with Chern numbers $C$=1 and $C$=2, respectively. Our study not only exposes the electric field induced exotic properties of MnBi2Te4 monolayer, but also proposes novel materials to realize QAHE in ferromagnetic Janus semiconductors with electric polarization.

preprint2021arXiv

Enhanced superconductivity in C-S-H compounds at high pressure

Recently, the superconducting transition temperature Tc = 287 K has been experimentally obtained in the material composed of carbon, sulfur, and hydrogen under the high pressure of 267 GPa. The material structure is unknown yet, where the carbon and sulfur were added at a molar ratio of 1:1. Here, fixing the molar ratio of C : S = 1:1, we studied several possible C-S-H structures, and found a new stable structure C2S2H4 using the first-principles calculations. The C2S2H4 shows an insulator-to-metal transition and the superconducting ground state at the pressure of 64 GPa, and its Tc can reach 16.5 K at 300 GPa. In addition, we found another stable structure of C2S3H4, whose Tc is 47.4 K at 300 GPa. The calculations show that the added S atom in C2S3H4 breaks part of C-H bonds in C2S2H4, makes the vibration of H atom at a lower frequency, and thus enhances the electro-phonon coupling and Tc. Our results suggest that the molar ratio of C:S lower than 1:1 in the C-S-H systems may be favorable to enhance Tc. This can be useful to figure out the structure of the C-S-H material with room temperature Tc in the recent experiment.

preprint2021arXiv

Fundamental Distinction between Intrinsic and Extrinsic Nonlinear Thermal Hall Effects

We theoretically investigated the fundamental distinction between intrinsic and extrinsic nonlinear thermal Hall effect in the presence of disorder at the second-order response to the temperature gradient in terms of the semi-classical Boltzmann equation. We found that, at low temperatures, the intrinsic contribution of the nonlinear thermal Hall conductivity is proportional to the square of temperature, whereas the extrinsic contributions (side-jump and skew-scattering) are independent of temperature. This distinct dependency on temperature provide a new approach to readily distinguish the intrinsic and extrinsic contributions. Specifically, we analysed the nonlinear thermal Hall effect for a tilted two-dimensional massive Dirac material. In particular, we showed that when the Fermi energy is located at the Dirac point, the signal is solely from the intrinsic mechanism; when the Fermi energy is higher, the extrinsic contributions are dominant, which are uncovered to be two to three orders of magnitude larger than the intrinsic contribution.

preprint2021arXiv

T-carbon: experiments, properties, derivatives and potential applications

Carbon is an extremely versatile element and carbon allotropes are very useful in all aspects of life and scientific research. T-carbon is a novel carbon allotrope with many appealing properties. Since the proposal of T-carbon, there has been a lot of intensive studies devoted to its physical, chemical, optical, magnetic, thermoelectrical and topological properties and possible applications in diverse areas in recent years. In this review, we provide a comprehensive review on the advances of the experiments, intriguing properties and various potential applications of T-carbon in energy storage, optoelectronics, thermoelectrics, topological states, etc. As intercalation and doping are effective methods to modify the electronic structure of materials or even to convert into sparkly different new structures or phases, we also discuss different atom doped T-carbon, which exhibit more intriguing properties and lead to promising potential applications in solar cells, photocatalysis, magnetism, superconductivity and so on. In addition, it is interesting to mention that the hydrogenated T-carbon molecules werefound to be possibly related to the physical origin of the UV extinction feature in interstellar medium that has been the half-century long unsolved puzzle. Many novel derivative structures either derived or inspired from T-carbon are included as well. Finally, we give prospects and outlook for future directions of study on T-carbon and related structures.

preprint2021arXiv

Two-dimensional topological superconductivity candidate in van der Waals layered material

Two-dimensional (2D) topological superconductors are highly desired because they not only offer opportunities for exploring novel exotic quantum physics, but also possesses potential applications in quantum computation. However, there are few reports on 2D superconductors, let alone topological superconductors. Here, we find a 2D monolayer W$_2$N$_3$, which can be exfoliated from its real van der Waals bulk material with much lower exfoliation energy than MoS$_2$, to be a topological metal with exotic topological states at different energy level. Due to the Van Hove singularities, the density of states near Fermi level are high, making the monolayer a compensate metal. Moreover, the monolayer W$_2$N$_3$ is unveiled to be a superconductor with the superconducting transition temperature Tc $\sim$ 22 K and a superconducting gap of about 5 meV based on the anisotropic Migdal-Eliashberg formalism, arising from the strong electron-phonon coupling around the $Γ$ point. Because of the strong electron and lattice coupling, the monolayer displays a non-Fermi liquid behavior in its normal states at temperatures lower than 80 K, where the specific heat exhibit T$^3$ behavior and the Wiedemann-Franz law dramatically violates. Our findings not only provide the platform to study the emergent phenomena in 2D topological superconductors, but also open a door to discover more 2D high-temperature topological superconductors in van der Waals materials.

preprint2020arXiv

Antiferromagnetic and Electric Polarized States in Two-Dimensional Janus Semiconductor Fe2Cl3I3

Two-dimensional (2D) Janus semiconductors with mirror asymmetry can introduce novel properties, such as large spin-orbit coupling (SOC) and normal piezoelectric polarization, which have attracted a great interest for their potential applications. Inspired by the recently fabricated 2D ferromagnetic (FM) semiconductor CrI3, a stable 2D (in x-y plane) antiferromagnetic (AFM) Janus semiconductor Fe2Cl3I3 with normal sublattice magnetization (m//z) is obtained by density functional theory calculations. By applying a tensile strain, the four magnetic states sequentially occur: AFM with m//z of sublattice, AFM with m//xy of sublattice, FM with m//xy, and FM with m//z. Such novel magnetic phase diagram driven by strain can be well understood by the spin-spin interactions including the third nearest-neighbor hoppings with the single-ion anisotropy, in which the SOC of I atoms is found to play an essential role. In addition, the electric polarization of Fe2Cl3I3 preserves with strain due to the broken inversion symmetry. Our results predict the rare Janus material Fe2Cl3I3 as an example of 2D semiconductors with both spin and charge polarizations, and reveal the highly sensitive strain-controlled magnetic states and magnetization direction, which highlights the 2D magnetic Janus semiconductor as a new platform to design spintronic materials.

preprint2020arXiv

Flat Band and Hole-induced Ferromagnetism in a Novel Carbon Monolayer

In recent experiments, superconductivity and correlated insulating states were observed in twisted bilayer graphene (TBG) with small magic angles, which highlights the importance of the flat bands near Fermi energy. However, the moiré pattern of TBG consists of more than ten thousand carbon atoms that is not easy to handle with conventional methods. By density functional theory calculations, we obtain a flat band at E$_F$ in a novel carbon monolayer coined as cyclicgraphdiyne with the unit cell of eighteen atoms. By doping holes into cyclicgraphdiyne to make the flat band partially occupied, we find that cyclicgraphdiyne with 1/8, 1/4, 3/8 and 1/2 hole doping concentration shows ferromagnetism (half-metal) while the case without doping is nonmagnetic, indicating a hole-induced nonmagnetic-ferromagnetic transition. The calculated conductivity of cyclicgraphdiyne with 1/8, 1/4 and 3/8 hole doping concentration is much higher than that without doping or with 1/2 hole doping. These results make cyclicgraphdiyne really attractive. By studying several carbon monolayers, we find that a perfect flat band may occur in the lattices with both separated or corner-connected triangular motifs with only including nearest-neighboring hopping of electrons, and the dispersion of flat band can be tuned by next-nearest-neighboring hopping. Our results shed insightful light on the formation of flat band in TBG. The present study also poses an alternative way to manipulate magnetism through doping flat band in carbon materials.

preprint2020arXiv

Hexagonal Warping Induced Nonlinear Planar Nernst Effect in Nonmagnetic Topological Insulators

We propose theoretically a new effect, i.e. nonlinear planar Nernst effect (NPNE), in nonmagnetic topological insulator (TI) Bi2Te3 in the presence of an in-plane magnetic field. We find that the Nernst current scales quadratically with temperature gradient but linearly with magnetic field and exhibits a cosine dependence of the orientation of the magnetic field with respect to the direction of the temperature gradient. The NPNE has a quantum origin arising from the conversion of a nonlinear transverse spin current to a charge current due to a joint result of hexagonal warping effect, spin-momentum locking, and the time-reversal symmetry breaking induced by the magnetic field.

preprint2020arXiv

Kagome quantum anomalous Hall effect with high Chern number and large band gap

Due to the potential applications in the low-power-consumption spintronic devices, the quantum anomalous Hall effect (QAHE) has attracted tremendous attention in past decades. However, up to now, QAHE was only observed experimentally in topological insulators with Chern numbers C= 1 and 2 at very low temperatures. Here, we propose three novel two-dimensional stable kagome ferromagnets Co3Pb3S2, Co3Pb3Se2and Co3Sn3Se2that can realize QAHE with high Chern number of |C|=3. Monolayers Co3Pb3S2, Co3Pb3Se2 and Co3Sn3Se2 possess the large band gap of 70, 77 and 63 meV with Curie temperature TC of 51, 42 and 46 K, respectively. By constructing a heterostructure Co3Sn3Se2/MoS2, its TC is enhanced to 60 K and the band gap keeps about 60 meV due to the tensile strain of 2% at the interface. For the bilayer compound Co6Sn5Se4, it becomes a half-metal, with a relatively flat plateau in its anomalous Hall conductivity corresponding to |C| = 3 near the Fermi level. Our results provide new topological nontrivial systems of kagome ferromagnetic monolayers and heterostructrues possessing QAHE with high Chern number |C| = 3 and large band gaps.

preprint2020arXiv

Large Family of Two-Dimensional Ferroelectric Metals Discovered via Machine Learning

Ferroelectricity and metallicity are usually believed not to coexist because conducting electrons would screen out static internal electric fields. In 1965, Anderson and Blount proposed the concept of 'ferroelectric metal', however, it is only until recently that very rare ferroelectric metals were reported. Here, by combining high-throughput ab initio calculations and data-driven machine learning method with new electronic orbital based descriptors, we systematically investigated a large family (2,964) of two-dimensional (2D) bimetal phosphates, and discovered 60 stable ferroelectrics with out-of-plane polarization, including 16 ferroelectric metals and 44 ferroelectric semiconductors that contain seven multiferroics. The ferroelectricity origins from spontaneous symmetry breaking induced by the opposite displacements of bimetal atoms, and the full-d-orbital coinage metal elements cause larger displacements and polarization than other elements. For 2D ferroelectric metals, the odd electrons per unit cell without spin polarization may lead to a half-filled energy band around Fermi level and is responsible for the metallicity. It is revealed that the conducting electrons mainly move on a single-side surface of the 2D layer, while both the ionic and electric contributions to polarization come from the other side and are vertical to the above layer, thereby causing the coexistence of metallicity and ferroelectricity. Van der Waals heterostructures based on ferroelectric metals may enable the change of Schottky barrier height or the Schottky-Ohmic contact type and induce a dramatic change of their vertical transport properties. Our work greatly expands the family of 2D ferroelectric metals and will spur further exploration of 2D ferroelectric metals.

preprint2020arXiv

Superconductivity in sodium-doped T-carbon

T-carbon has been proposed as a new carbon allotrope in 2011, which was successfully synthesized in recent experiments. Because of its fluffy structure, several kinds of atoms can be intercalated into T-carbon, making it a versatile candidate in various applications such as hydrogen storage, perovskite solar cells, lithium ion batteries, thermoelectrics, photocatalyst, etc. Here we show that superconductivity can appear in Na-doped T-carbon with superconducting transition temperature Tc of 11 K at ambient pressure, and Tc can be enhanced to 19 K under pressure of 14 GPa, which results from an enhancement of the electron-phonon coupling due to the shift of the phonon spectral weight to lower frequencies with the increase of pressure. The calculations on specific heat and electrical and thermal conductivities show that the normal state of the Na-doped T-carbon superconductor reveals a non-Fermi liquid behavior. The prediction of superconductivity in Na-doped T-carbon would spur great interest both experimentally and theoretically to explore novel carbon-based superconductors.

preprint2020arXiv

Tangent-Space Gradient Optimization of Tensor Network for Machine Learning

The gradient-based optimization method for deep machine learning models suffers from gradient vanishing and exploding problems, particularly when the computational graph becomes deep. In this work, we propose the tangent-space gradient optimization (TSGO) for the probabilistic models to keep the gradients from vanishing or exploding. The central idea is to guarantee the orthogonality between the variational parameters and the gradients. The optimization is then implemented by rotating parameter vector towards the direction of gradient. We explain and testify TSGO in tensor network (TN) machine learning, where the TN describes the joint probability distribution as a normalized state $\left| ψ\right\rangle $ in Hilbert space. We show that the gradient can be restricted in the tangent space of $\left\langle ψ\right.\left| ψ\right\rangle = 1$ hyper-sphere. Instead of additional adaptive methods to control the learning rate in deep learning, the learning rate of TSGO is naturally determined by the angle $θ$ as $η= \tan θ$. Our numerical results reveal better convergence of TSGO in comparison to the off-the-shelf Adam.

preprint2020arXiv

Tinselenidene: a Two-dimensional Auxetic Material with Ultralow Lattice Thermal Conductivity and Ultrahigh Hole Mobility

By means of extensive ab initio calculations, a new two-dimensional (2D) atomic material tin selenide monolayer (coined as tinselenidene) is predicted to be a semiconductor with an indirect gap (1.45 eV) and a high hole mobility (of order 10000 cm2V-1S-1), and will bear an indirect-direct gap transition under a rather low strain (<0.5 GPa). Tinselenidene has a very small Young's modulus (20-40 GPa) and an ultralow lattice thermal conductivity (<3 Wm-1K-1 at 300 K), making it probably the most flexible and most heat-insulating material in known 2D atomic materials. In addition, tinseleniden has a large negative Poisson's ratio of -0.17, thus could act as a 2D auxetic material. With these intriguing properties, tinselenidene could have wide potential applications in thermoelectrics, nanomechanics and optoelectronics.

preprint2020arXiv

Toward understanding physical origin of 2175Å extinction bump in interstellar medium

The 2175 Å ultraviolet (UV) extinction bump in interstellar medium (ISM) of the Milky Way was discovered in 1965. After intensive exploration of more than a half century, however, its exact origin still remains a big conundrum that is being debated. Here we propose a mixture model by which the extinction bump in ISM is argued possibly relevant to the clusters of hydrogenated T-carbon (HTC) molecules (C40H16) that have intrinsically a sharp absorption peak at the wavelength 2175 Å. By linearly combining the calculated absorption spectra of HTC mixtures, graphite, MgSiO3 and Fe2SiO4, we show that the UV extinction curves of optional six stars can be nicely fitted. This present work poses an alternative explanation toward understanding the physical origin of the 2175 Å extinction bump in ISM of the Milky Way.

preprint2020arXiv

Two-dimensional magnetic semiconductors with room Curie temperatures

We propose two-dimensional (2D) Ising-type ferromagnetic semiconductors TcSiTe3, TcGeSe3, and TcGeTe3 with high Curie temperatures around 200-0500 K. Owing to large spin-orbit couplings, the large magnetocrystalline anisotropy energy (MAE), large anomalous Hall conductivity, and large magneto-optical Kerr effect were discovered in these intriguing 2D materials. By comparing all possible 2D MGeTe3 materials (M = 3d, 4d, 5d transition metals), we found a large orbital moment around 0.5 $μ$B per atom and a large MAE for TcGeTe3. The large orbital moments are revealed to be from the comparable crystal fields and electron correlations in these Tc-based 2D materials. The microscopic mechanism of the high Curie temperature is also addressed. Our findings reveal the unique magnetic behaviors of 2D Tc-based materials and present a family of 2D ferromagnetic semiconductors with large MAE and Kerr rotation angles that would have wide applications in designing spintronic devices.

preprint2019arXiv

Accelerated Discovery of Two-Dimensional Optoelectronic Octahedral Oxyhalides via High-Throughput Ab Initio Calculations and Machine Learning

Traditional trial-and-error methods are obstacles for large-scale searching of new optoelectronic materials. Here, we introduce a method combining high-throughput ab initio calculations and machine-learning approaches to predict two-dimensional octahedral oxyhalides with improved optoelectronic properties. We develop an effective machine-learning model based on an expansive dataset generated from density functional calculations including the geometric and electronic properties of 300 two-dimensional octahedral oxyhalides. Our model accelerates the screening of potential optoelectronic materials of 5,000 two-dimensional octahedral oxyhalides. The distorted stacked octahedral factors proposed in our model play essential roles in the machine-learning prediction. Several potential two-dimensional optoelectronic octahedral oxyhalides with moderate band gaps, high electron mobilities, and ultrahigh absorbance coefficients are successfully hypothesized.

preprint2019arXiv

Generative Tensor Network Classification Model for Supervised Machine Learning

Tensor network (TN) has recently triggered extensive interests in developing machine-learning models in quantum many-body Hilbert space. Here we purpose a generative TN classification (GTNC) approach for supervised learning. The strategy is to train the generative TN for each class of the samples to construct the classifiers. The classification is implemented by comparing the distance in the many-body Hilbert space. The numerical experiments by GTNC show impressive performance on the MNIST and Fashion-MNIST dataset. The testing accuracy is competitive to the state-of-the-art convolutional neural network while higher than the naive Bayes classifier (a generative classifier) and support vector machine. Moreover, GTNC is more efficient than the existing TN models that are in general discriminative. By investigating the distances in the many-body Hilbert space, we find that (a) the samples are naturally clustering in such a space; and (b) bounding the bond dimensions of the TN's to finite values corresponds to removing redundant information in the image recognition. These two characters make GTNC an adaptive and universal model of excellent performance.

preprint2019arXiv

Lecture Notes of Tensor Network Contractions

Tensor network (TN), a young mathematical tool of high vitality and great potential, has been undergoing extremely rapid developments in the last two decades, gaining tremendous success in condensed matter physics, atomic physics, quantum information science, statistical physics, and so on. In this lecture notes, we focus on the contraction algorithms of TN as well as some of the applications to the simulations of quantum many-body systems. Starting from basic concepts and definitions, we first explain the relations between TN and physical problems, including the TN representations of classical partition functions, quantum many-body states (by matrix product state, tree TN, and projected entangled pair state), time evolution simulations, etc. These problems, which are challenging to solve, can be transformed to TN contraction problems. We present then several paradigm algorithms based on the ideas of the numerical renormalization group and/or boundary states, including density matrix renormalization group, time-evolving block decimation, coarse-graining/corner tensor renormalization group, and several distinguished variational algorithms. Finally, we revisit the TN approaches from the perspective of multi-linear algebra (also known as tensor algebra or tensor decompositions) and quantum simulation. Despite the apparent differences in the ideas and strategies of different TN algorithms, we aim at revealing the underlying relations and resemblances in order to present a systematic picture to understand the TN contraction approaches.

preprint2019arXiv

Quantum Compressed Sensing with Unsupervised Tensor-Network Machine Learning

We propose tensor-network compressed sensing (TNCS) by combining the ideas of compressed sensing, tensor network (TN), and machine learning, which permits novel and efficient quantum communications of realistic data. The strategy is to use the unsupervised TN machine learning algorithm to obtain the entangled state $|Ψ\rangle$ that describes the probability distribution of a huge amount of classical information considered to be communicated. To transfer a specific piece of information with $|Ψ\rangle$, our proposal is to encode such information in the separable state with the minimal distance to the measured state $|Φ\rangle$ that is obtained by partially measuring on $|Ψ\rangle$ in a designed way. To this end, a measuring protocol analogous to the compressed sensing with neural-network machine learning is suggested, where the measurements are designed to minimize uncertainty of information from the probability distribution given by $|Φ\rangle$. In this way, those who have $|Φ\rangle$ can reliably access the information by simply measuring on $|Φ\rangle$. We propose q-sparsity to characterize the sparsity of quantum states and the efficiency of the quantum communications by TNCS. The high q-sparsity is essentially due to the fact that the TN states describing nicely the probability distribution obey the area law of entanglement entropy. Testing on realistic datasets (hand-written digits and fashion images), TNCS is shown to possess high efficiency and accuracy, where the security of communications is guaranteed by the fundamental quantum principles.

preprint2019arXiv

Reentrance of Topological Phase in Spin-1 Frustrated Heisenberg Chain

For the Haldane phase, the magnetic field usually tends to break the symmetry and drives the system into a topologically trivial phase. Here, we report a novel reentrance of the Haldane phase at zero temperature in the spin-1 antiferromagnetic Heisenberg model on sawtooth chain. A partial Haldane phase is induced by the magnetic field, which is the combination of the Haldane state in one sublattice and a ferromagnetically ordered state in the other sublattice. Such a partial topological order is a result of the zero-temperature entropy due to quantum fluctuations caused by geometrical frustration.

preprint2016arXiv

Diverse anisotropy of phonon transport in two-dimensional IV-VI compounds: A comparative study

New classes two-dimensional (2D) materials beyond graphene, including layered and non-layered, and their heterostructures, are currently attracting increasing interest due to their promising applications in nanoelectronics, optoelectronics and clean energy, where thermal transport property is one of the fundamental physical parameters. In this paper, we systematically investigated the phonon transport properties of 2D orthorhombic group IV-VI compounds of $GeS$, $GeSe$, $SnS$ and $SnSe$ by solving the Boltzmann transport equation (BTE) based on first-principles calculations. Despite the similar puckered (hinge-like) structure along the armchair direction as phosphorene, the four monolayer compounds possess diverse anisotropic properties in many aspects, such as phonon group velocity, Young's modulus and lattice thermal conductivity ($κ$), etc. Especially, the $κ$ along the zigzag and armchair directions of monolayer $GeS$ shows the strongest anisotropy while monolayer $SnS$ and $SnSe$ shows an almost isotropy in phonon transport. The origin of the diverse anisotropy is fully studied and the underlying mechanism is discussed in detail. With limited size, the $κ$ could be effectively lowered, and the anisotropy could be effectively modulated by nanostructuring, which would extend the applications in nanoscale thermoelectrics and thermal management. Our study offers fundamental understanding of the anisotropic phonon transport properties of 2D materials, and would be of significance for further study, modulation and aplications in emerging technologies.

preprint2016arXiv

Layer dependence of geometric, electronic and piezoelectric properties of SnSe

By means of first-principles calculations, we explore systematically the geometric, electronic and piezoelectric properties of multilayer SnSe. We find that these properties are layer-dependent, indicating that the interlayer interaction plays an important role. With increasing the number of SnSe layers from 1 to 6, we observe that the lattice constant decreases from 4.27 $\mathring{A}$ to 4.22 $\mathring{A}$ along zigzag direction, and increases from 4.41 $\mathring{A}$ to 4.51 $\mathring{A}$ along armchair direction close to the bulk limit (4.21 $\mathring{A}$ and 4.52 $\mathring{A}$, respectively); the band gap decreases from 1.45 eV to 1.08 eV, approaching the bulk gap 0.95 eV. Although the monolayer SnSe exhibits almost symmetric geometric and electronic structures along zigzag and armchair directions, bulk SnSe is obviously anisotropic, showing that the stacking of layers enhances the anisotropic character of SnSe. As bulk and even-layer SnSe have inversion centers, they cannot exhibit piezoelectric responses. However, we show that the odd-layer SnSe have piezoelectric coefficients much higher than those of the known piezoelectric materials, suggesting that the odd-layer SnSe is a good piezoelectric material.

preprint2016arXiv

Quantum phase transition, universality and scaling behaviors in the spin-1/2 Heisenberg model with ferromagnetic and antiferromagnetic competing interactions on honeycomb lattice

The quantum phase transition, scaling behaviors, and thermodynamics in the spin-1/2 quantum Heisenberg model with antiferromagnetic coupling $J>0$ in armchair direction and ferromagnetic interaction $J'<0$ in zigzag direction on a honeycomb lattice are systematically studied using the continuous-time quantum Monte Carlo method. By calculating the Binder ratio $Q_{2}$ and spin stiffness $ρ$ in two directions for various coupling ratio $α=J'/J$ under different lattice sizes, we found that a quantum phase transition from the dimerized phase to the stripe phase occurs at the quantum critical point $α_c=-0.93$. Through the finite-size scaling analysis on $Q_{2}$, $ρ_{x}$ and $ρ_{y}$, we determined the critical exponent related to the correlation length $ν$ to be 0.7212(8), implying that this transition falls into a classical Heisenberg O(3) universality. A zero magnetization plateau is observed in the dimerized phase, whose width decreases with increasing $α$. A phase diagram in the coupling ratio $α$-magnetic field $h$ plane is obtained, where four phases, including dimerized, stripe, canted stripe and polarized phases are identified. It is also unveiled that the temperature dependence of the specific heat $C(T)$ for different $α$'s intersects precisely at one point, similar to that of liquid $^{3}$He under different pressures and several magnetic compounds under various magnetic fields. The scaling behaviors of $Q_{2}$, $ρ$ and $C(T)$ are carefully analyzed. The susceptibility is well compared with the experimental data to give the magnetic parameters of both compounds.

preprint2016arXiv

Spin ordered ground state and thermodynamic behaviors of the spin-3/2 kagome Heisenberg antiferromagnet

Three different tensor network optimization algorithms are employed to accurately determine the ground state and thermodynamic properties of the spin-3/2 kagome Heisenberg antiferromagnet. We found that the $\sqrt{3} \times \sqrt{3}$ state, rather than the $q = 0$ state, is the ground state of this system, and such an ordered state is melted at any finite temperature, thereby clarifying the existing experimental controversies. A 1/3-magnetization plateau in the magnetic curve is observed, being consistent with the experimental observation. The absence of a zero-magnetization plateau indicates a gapless spin excitation that is further supported by the thermodynamic asymptotic behaviors of the susceptibility and specific heat. At low temperatures, the specific heat is shown to exhibit a $\sqrt{T}$ behavior, and the susceptibility approaches a finite constant as $T\rightarrow 0$. Our TN results of thermodynamic properties are compared with those from high temperature series expansion. In addition, we observe a quantum phase transition between $q = 0$ and $\sqrt{3}\times\sqrt{3}$ states in a spin-3/2 kagome XXZ model at the critical point $Δ_c = 0.408$. This study provides reliable and useful information for further explorations on high spin kagome physics.

preprint2015arXiv

Anisotropic intrinsic lattice thermal conductivity of phosphorene from first principles

Phosphorene, the single layer counterpart of black phosphorus, is a novel two-dimensional semiconductor with high carrier mobility and a large fundamental direct band gap, which has attracted tremendous interest recently. Its potential applications in nano-electronics and thermoelectrics call for a fundamental study of the phonon transport. Here, we calculate the intrinsic lattice thermal conductivity of phosphorene by solving the phonon Boltzmann transport equation (BTE) based on first-principles calculations. The thermal conductivity of phosphorene at $300\,\mathrm{K}$ is $30.15\,\mathrm{Wm^{-1}K^{-1}}$ (zigzag) and $13.65\,\mathrm{Wm^{-1}K^{-1}}$ (armchair), showing an obvious anisotropy along different directions. The calculated thermal conductivity fits perfectly to the inverse relation with temperature when the temperature is higher than Debye temperature ($Θ_D = 278.66\,\mathrm{K}$). In comparison to graphene, the minor contribution around $5\%$ of the ZA mode is responsible for the low thermal conductivity of phosphorene. In addition, the representative mean free path (MFP), a critical size for phonon transport, is also obtained.

preprint2015arXiv

Highly efficient light management for perovskite solar cells

Organic-inorganic halide perovskite solar cells have enormous potential to impact the existing photovoltaic industry. As realizing a higher conversion efficiency of the solar cell is still the most crucial task, a great number of schemes were proposed to minimize the carrier loss by optimizing the electrical properties of the perovskite solar cells. Here, we focus on another significant aspect that is to minimize the light loss by optimizing the light management to gain a high efficiency for perovskite solar cells. In our scheme, the slotted and inverted prism structured SiO2 layers are adopted to trap more light into the solar cells, and a better transparent conducting oxide layer is employed to reduce the parasitic absorption. For such an implementation, the efficiency and the serviceable angle of the perovskite solar cell can be promoted impressively. This proposal would shed new light on developing the high-performance perovskite solar cells.

preprint2015arXiv

Identification of Degeneracy and Criticality of Two-Dimensional Statistical and Quantum Systems by the Boundary States of Tensor Networks

We propose a systematic scheme to reach the properties of two-dimensional (2D) statistical and quantum systems by studying the effective (1+1)-dimensional theory that is constructed from the tensor network representation. On on hand, we discover that the degeneracy of the 2D system can be determined by the purity of the boundary thermal state, which is the density operator of the effective theory at zero (effective) temperature. On the other hand, we find that the gapped (or critical) 2D system leads to a gapped (or critical) effective (1+1)-dimensional theory whose criticality can be accessed by the entanglement entropy $S$ of its ground state dubbed as boundary pure state. We also uncover that for the critical systems, $S$ obeys the same logarithmic law as that found in the critical 1D quantum chains, which reads $S = (κc/6)\log_2 D + const.$, with $c$ the central charge and $κ$ a constant related to the scaling property of the correlation length $ξ$ as $ξ\sim D^κ$. Such a scaling law presents an efficient way to characterize the critical universality class of the original 2D systems. An important implication of our work is that many well-established theories for 1D quantum chains become available for studying 2D systems with the help of the proposed lower dimensional correspondence.

preprint2015arXiv

Low cost and high performance light trapping structure for thin-film solar cells

Nano-scaled dielectric and metallic structures are popular light tapping structures in thin-film solar cells. However, a large parasitic absorption in those structures is unavoidable. Most schemes based on such structures also involve the textured active layers that may bring undesirable degradation of the material quality. Here we propose a novel and cheap light trapping structure based on the prism structured SiO2 for thin-film solar cells, and a flat active layer is introduced purposefully. Such a light trapping structure is imposed by the geometrical shape optimization to gain the best optical benefit. By examining our scheme, it is disclosed that the conversion efficiency of the flat a-Si:H thin-film solar cell can be promoted to exceed the currently certified highest value. As the cost of SiO2-based light trapping structure is much cheaper and easier to fabricate than other materials, this proposal would have essential impact and wide applications in thin-film solar cells.

preprint2015arXiv

Multivariable Scaling for the Anomalous Hall Effect

We derive a general scaling relation for the anomalous Hall effect in ferromagnetic metals involving multiple competing scattering mechanisms, described by a quadratic hypersurface in the space spanned by the partial resistivities. We also present experimental findings, which show strong deviation from previously found scaling forms when different scattering mechanism compete in strength but can be nicely explained by our theory.

preprint2015arXiv

Phase diagram and exotic spin-spin correlations of anisotropic Ising model on the Sierpiński gasket

The anisotropic antiferromagnetic Ising model on the fractal Sierpiński gasket is intensively studied, and a number of exotic properties are disclosed. The ground state phase diagram in the plane of magnetic field-interaction of the system is obtained. The thermodynamic properties of the three plateau phases are probed by exploring the temperature-dependence of magnetization, specific heat, susceptibility and spin-spin correlations. No phase transitions are observed in this model. In the absence of a magnetic field, the unusual temperature dependence of the spin correlation length is obtained with $0 \leq$J$_b/$J$_a<1$, and an interesting crossover behavior between different phases at J$_b/$J$_a=1$ is unveiled, whose dynamics can be described by the J$_b/$J$_a$-dependence of the specific heat, susceptibility and spin correlation functions. The exotic spin-spin correlation patterns that share the same special rotational symmetry as that of the Sierpiński gasket are obtained in both the $1/3$ plateau disordered phase and the $5/9$ plateau partially ordered ferrimagnetic phase. Moreover, a quantum scheme is formulated to study the thermodynamics of the fractal Sierpiński gasket with Heisenberg interactions. We find that the unusual temperature dependence of the correlation length remains intact in a small quantum fluctuation.

preprint2015arXiv

Quantum phase transitions and string orders in the spin-1/2 Heisenberg-Ising alternating chain with Dzyaloshinskii-Moriya interaction

Quantum phase transitions (QPTs) and the ground-state phase diagram of the spin-1/2 Heisenberg-Ising alternating chain (HIAC) with uniform Dzyaloshinskii-Moriya (DM) interaction are investigated by a matrix-product-state (MPS) method. By calculating the odd- and even-string order parameters, we recognize two kinds of Haldane phases, i.e., the odd- and even-Haldane phases. Furthermore, doubly degenerate entanglement spectra on odd and even bonds are observed in odd- and even-Haldane phases, respectively. A rich phase diagram including four different phases, i.e., an antiferromagnetic (AF), AF stripe, odd- and even-Haldane phases, is obtained. These phases are found to be separated by continuous QPTs: the topological QPT between the odd- and even-Haldane phases is verified to be continuous and corresponds to conformal field theory with central charge $c$=1; while the rest phase transitions in the phase diagram are found to be $c$=1/2. We also revisit, with our MPS method, the exactly solvable case of HIAC model with DM interactions only on odd bonds, and find that the even-Haldane phase disappears, but the other three phases, i.e., the AF, AF stripe, and odd-Haldane phases, still remain in the phase diagram. We exhibit the evolution of the even-Haldane phase by tuning the DM interactions on the even bonds gradually.

preprint2015arXiv

Simplex valence-bond crystal in the spin-1 kagome Heisenberg antiferromagnet

We investigate the ground state properties of a spin-1 kagome antiferromagnetic Heisenberg model using tensor-network (TN) methods. We obtain the energy per site {$e_0=-1.41090(2)$ with $D^*=8$ multiplets retained (i.e., a bond dimension of $D=24$), and $e_0=-1.4116(4)$ from large-$D$ extrapolation,} by accurate TN calculations directly in the thermodynamic limit. The symmetry between the two kinds of triangles is spontaneously broken, with a relative energy difference of $δ\approx$ 19\%, i.e, there is a trimerization (simplex) valence-bond order in the ground state. The spin-spin, dimer-dimer, and chirality-chirality correlation functions are found to decay exponentially with a rather short correlation length, showing that the ground state is gapped. We thus identify the ground state be a simplex valence-bond crystal (SVBC). We also discuss the spin-1 bilinear-biquadratic Heisenberg model on a kagome lattice, and determine its ground state phase diagram. Moreover, we implement non-abelian symmetries, here spin SU(2), in the TN algorithm, which improves the efficiency greatly and provides insight into the tensor structures.

preprint2014arXiv

Evidence of the side jump mechanism in the anomalous Hall effect in paramagnets

Persistent confusion has existed between the intrinsic (Berry curvature) and the side jump mechanisms of anomalous Hall effect (AHE) in ferromagnets. We provide unambiguous identification of the side jump mechanism, in addition to the skew scattering contribution in epitaxial paramagnetic Ni$_{34}$Cu$_{66}$ thin films, in which the intrinsic contribution is by definition excluded. Furthermore, the temperature dependence of the AHE further reveals that the side jump mechanism is dominated by the elastic scattering.

preprint2014arXiv

Featureless quantum spin liquid, 1/3-magnetization plateau state and exotic thermodynamic properties of spin-1/2 frustrated Heisenberg antiferromagnet on an infinite Husimi lattice

By utilizing tensor-network-based methods, we investigate the zero- and finite-temperature properties of the spin-1/2 Heisenberg antiferromagnetic (HAF) model on an infinite Husimi lattice that contains 3/2 sites per triangle. The ground state of this model is found to possess vanishing local magnetization and is featureless; the spin-spin and dimer-dimer correlation functions are verified to decay exponentially; and its ground-state energy per site is determined to be $e_0=-0.4343(1)$, which is very close to that [$e_0=-0.4386(5)$] of the intriguing kagome HAF model. The magnetization curve shows the absence of a zero-magnetization plateau, implying a gapless excitation. A $1/3$-magnetization plateau with spin up-up-down state is observed, which is selected and stabilized by quantum fluctuations. A ground state phase diagram under magnetic fields is presented. Moreover, both magnetic susceptibility and the specific heat are studied, whose low-temperature behaviors reinforce the conclusion that the HAF model on the infinite Husimi lattice owns a gapless and featureless spin liquid ground state.

preprint2014arXiv

Hinge-like structure induced unusual properties of black phosphorus and new strategies to improve the thermoelectric performance

We systematically investigated the geometric, electronic and thermoelectric (TE) properties of bulk black phosphorus (BP) under strain. The hinge-like structure of BP brings unusual mechanical responses such as anisotropic Young's modulus and negative Poisson's ratio. A sensitive electronic structure of BP makes it transform among metal, direct and indirect semiconductors under strain. The maximal figure of merit $ZT$ of BP is found to be 0.72 at $800\,\mathrm{K}$ that could be enhanced to 0.87 by exerting an appropriate strain, revealing BP could be a potential medium-high temperature TE material. Such strain-induced enhancements of TE performance are often observed to occur at the boundary of the direct-indirect band gap transition, which can be attributed to the increase of degeneracy of energy valleys at the transition point. By comparing the structure of BP with SnSe, a family of potential TE materials with hinge-like structure are suggested. This study not only exposes various novel properties of BP under strain, but also proposes effective strategies to seek for better TE materials.

preprint2014arXiv

New strategy to promote conversion efficiency using high-index nanostructures in thin-film solar cells

Nano-scaled metallic or dielectric structures may provide various ways to trap light into thin-film solar cells for improving the conversion efficiency. In most schemes, the textured active layers are involved into light trapping structures that can provide perfect optical benefits but also bring undesirable degradation of electrical performance. Here we propose a novel approach to design high-performance thin-film solar cells. In our strategy, a flat active layer is adopted for avoiding electrical degradation, and an optimization algorithm is applied to seek for an optimized light trapping structure for the best optical benefit. As an example, we show that the efficiency of a flat a-Si:H thin-film solar cell can be promoted close to the certified highest value. It is also pointed out that, by choosing appropriate dielectric materials with high refractive index (>3) and high transmissivity in wavelength region of 350nm-800nm, the conversion efficiency of solar cells can be further enhanced.

preprint2013arXiv

Electronic and magnetic properties of twisted graphene nanoribbon and Möbius strips: first-principles calculations

The geometrical, electronic, and magnetic properties of twisted zigzag-edged graphene nanoribbons (ZGNRs) and novel graphene Möbius strips (GMS) are systematically investigated using first-principles density functional calculations. The structures of ZGNRs and GMS are optimized, and their stabilities are examined. The molecular energy levels and the spin polarized density of states are calculated. It is found that for twisted ZGNRs, the atomic bonding energy decreases quadratically with the increase of the twisted angle, and the HOMO-LUMO gap are varying in a sine-like behavior with the twisted angle. The calculated spin densities reveal that the ZGNRs and GMS have antiferromagnetic ground states, which persist during the twisting. The spin flips on the zigzag edges of GMS are observed at some positions.

preprint2013arXiv

Honeycomb Heisenberg Spin Ladder: Unusual Ground State and Thermodynamic Properties

The unusual ground state and thermodynamic properties of spin-1/2 two-leg honeycomb (HC) spin ladder are systematically studied by jointly utilizing various analytical and numerical methods. The HC spin ladder is found to exhibit behaviors dramatically different from those of the conventional square spin ladder. A strong relevant term nJ and a half saturation magnetization plateau that can be attributed to the formation of diluted dimer states are observed in the HC ladder, both of which are absent in the square ladder. The ground state phase diagram of the HC spin ladder is identified, and the thermodynamic properties of the specific heat and susceptibility for different couplings are thoroughly explored, where two kinds of excitations are unveiled. The distinct Wilson ratios for both spin ladders at the lower critical fields are also obtained. Our calculated result is well fitted to the experimental data of the two-leg HC spin ladder compound [Cu2L1(N3)4]n.

preprint2013arXiv

Kosterlitz-Thouless phase transition and reentrance in an anisotropic 3-state Potts model on the generalized Kagome lattice

The unusual reentrant phenomenon is observed in the anisotropic 3-state Potts model on a gen- eralized Kagome lattice. By employing the linearized tensor renormalization group method, we find that the reentrance can appear in the region not only under a partial ordered phase as commonly known but also a phase without a local order parameter, which is uncovered to fall into the uni- versality of the Kosterlitz-Thouless (KT) type. The region of the reentrance depends strongly on the ratios of the next nearest couplings α = J2 /|J1 | and β = J3 /|J1 |. The phase diagrams in the plane of temperature versus β for different α are obtained. Through massive calculations, it is also revealed that the quasi-entanglement entropy can be used to accurately detect the KT transition temperature.

preprint2013arXiv

Superstructures of Cubic and Hexagonal Diamonds Comprising a Family of Novel $sp^3$ Superhard Carbon Allotropes

Superstructures of cubic and hexagonal diamonds (h- and c-diamond) comprising a family of stable diamond-like $sp^3$ hybridized novel carbon allotropes are proposed, which are referred to as U$_n$-carbon where $n \geq 2$ denotes the number of structural layers in a unit cell. The conventional h- and c-diamond are included in this family as members with $n=2$ and 3, respectively. U$_n$-carbon ($n=4-6$), which are unveiled energetically and thermodynamically more stable than h-diamond and possess remarkable kinetic stabilities, are shown to be insulators with indirect gaps of $5.6 \sim 5.8$ eV, densities of $ 3.5 \sim 3.6$ g/cm$^3$, bulk modulus of $4.3 \sim 4.4 \times 10^{2}$ GPa, and Vickers hardness of $92.9 \sim 97.5$ GPa even harder than h- and c-diamond. The simulated x-ray diffraction and Raman spectra are presented for experimental characterization. These new structures of carbon would have a compelling impact in physics, chemistry, materials science and geophysics.

preprint2013arXiv

Theory of Network Contractor Dynamics for Exploring Thermodynamic Properties of Two-dimensional Quantum Lattice Models

Based on the tensor network state representation, we develop a nonlinear dynamic theory coined as network contractor dynamics (NCD) to explore the thermodynamic properties of two-dimensional quantum lattice models. By invoking the rank-$1$ decomposition in the multi-linear algebra, the NCD scheme makes the contraction of the tensor network of the partition function be realized through a contraction of a local tensor cluster with vectors on its boundary. An imaginary-time-sweep algorithm for implementation of the NCD method is proposed for practical numerical simulations. We benchmark the NCD scheme on the square Ising model, which shows a great accuracy. Besides, the results on the spin-1/2 Heisenberg antiferromagnet on honeycomb lattice are disclosed in good agreement with the quantum Monte Carlo calculations. The quasi-entanglement entropy $S$, Lyapunov exponent $I^{lya}$ and loop character $I^{loop}$ are introduced within the dynamic scheme, which are found to display the ``nonlocality" near the critical point, and can be applied to determine the thermodynamic phase transitions of both classical and quantum systems.

preprint2012arXiv

Charged States and Band-Gap Narrowing in Codoped ZnO Nanowires for Enhanced Photoelectrochemical Responses

By means of first-principles calculations within the density functional theory, we study the structural and optical properties of codoped ZnO nanowires and compare them with those of the bulk and film. It is disclosed that the low negatively charged ground states of nitrogen related defects play a key role in the optical absorption spectrum tail that narrows the band-gap and enhances the photoelectrochemical response significantly. A strategy of uncompensated N, P and Ga codoping in ZnO nanowires is proposed to produce a red-shift of the optical absorption spectra further than the exclusive N doping and to get a proper formation energy with a high defect concentration and a suppressed recombination rate. In this way, the absorption of the visible light can be improved and the photocurrent can be raised. These observations are consistent with the existing experiments, which will be helpful to improve the photoelectrochemical responses for the wide-band-gap semiconductors especially in water splitting applications.

preprint2012arXiv

Low-energy effective theory and two distinct critical phases in a spin-1/2 frustrated three-leg spin tube

Motivated by the crystal structures of [(CuCl2tachH)3Cl]Cl2 and Ca3Co2O6, we develop a low-energy effective theory using the bosonization technique for a spin-1/2 frustrated three-leg spin tube with trigonal prism units in two limit cases. The features obtained with the effective theory are numerically elucidated by the density matrix renormalization group method. Three different quantum phases in the ground state of the system, say, one gapped dimerized phase and two distinct gapless phases, are identified, where the two gapless phases are found to have the conformal central charge c=1 and 3/2, respectively. Spin gaps, spin and dimer correlation functions, and the entanglement entropy are obtained. In particular, it is disclosed that the critical phase with c=3/2 is the consequence of spin frustrations, which might belong to the SU(2)k=2 Wess-Zumino-Witten-Novikov universality class, and is induced by the twist term in the bosonized Hamiltonian density.

preprint2012arXiv

Matrix Product State and Quantum Phase Transitions in the One-Dimensional Extended Quantum Compass Model

The matrix product state (MPS) is utilized to study the ground state properties and quantum phase transitions (QPTs) of the one-dimensional quantum compass model (QCM). The MPS wavefunctions are argued to be very efficient descriptions of QCM ground states, and are numerically determined by imaginary time projections. The ground state energy, correlations, quantum entanglement and its spectrum, local and nonlocal order parameters, etc., are calculated and studied in details. It is revealed that the bipartite and block entanglement entropies, as well as the nearest neighbor correlation functions can be used to detect the second-order QPTs, but not the first-order ones, while fidelity detections can recognize both. The entanglement spectrum is extracted from the MPS wavefunction, and found to be doubly degenerate in disordered phases of QCM, where non-local string order parameters exist. Moreover, with linearized tensor renormalization group method, the specific heat curves are evaluated and their low temperature behaviors are investigated.

preprint2012arXiv

Octagraphene as a Versatile Carbon Atomic Sheet for Novel Nanotubes, Unconventional Fullerenes and Hydrogen Storage

We study a versatile structurally favorable periodic $sp^2$-bonded carbon atomic planar sheet with $C_{4v}$ symmetry by means of the first-principles calculations. This carbon allotrope is composed of carbon octagons and squares with two bond lengths and is thus dubbed as octagraphene. It is a semimetal with the Fermi surface consisting of one hole and one electron pocket, whose low-energy physics can be well described by a tight-binding model of $π$-electrons. Its Young's modulus, breaking strength and Poisson's ratio are obtained to be 306 $N/m$, 34.4 $N/m$ and 0.13, respectively, which are close to those of graphene. The novel sawtooth and armchair carbon nanotubes as well as unconventional fullerenes can also be constructed from octagraphene. It is found that the Ti-absorbed octagraphene can be allowed for hydrogen storage with capacity around 7.76 wt%.

preprint2012arXiv

Optimized Decimation of Tensor Networks with Super-orthogonalization for Two-Dimensional Quantum Lattice Models

A novel algorithm based on the optimized decimation of tensor networks with super-orthogonalization (ODTNS) that can be applied to simulate efficiently and accurately not only the thermodynamic but also the ground state properties of two-dimensional (2D) quantum lattice models is proposed. By transforming the 2D quantum model into a three-dimensional (3D) closed tensor network (TN) comprised of the tensor product density operator and a 3D brick-wall TN, the free energy of the system can be calculated with the imaginary time evolution, in which the network Tucker decomposition is suggested for the first time to obtain the optimal lower-dimensional approximation on the bond space by transforming the TN into a super-orthogonal form. The efficiency and accuracy of this algorithm are testified, which are fairly comparable with the quantum Monte Carlo calculations. Besides, the present ODTNS scheme can also be applicable to the 2D frustrated quantum spin models with nice efficiency.

preprint2012arXiv

Phase Diagrams, Distinct Conformal Anomalies and Thermodynamics of Spin-1 Bond-Alternating Heisenberg Antiferromagnetic Chain in Magnetic Fields

The ground state phase diagrams and thermodynamic properties of a spin-1 bond-alternating Heisenberg antiferromagnetic chain with a single-ion anisotropy in longitudinal and transverse magnetic fields are investigated jointly by means of the infinite time evolving block decimation, the linearized tensor renormalization group and the density matrix renormalization group methods. It is found that in the magnetic field-bond alternating ratio plane six phases such as singlet dimer, Haldane, two Tomonaga-Luttinger liquid, 1/2 magnetization plateau and spin polarized phases are identified in a longitudinal field, while in a transverse field there are fives phases including singlet dimer, Haldane, Z$_2$ symmetry breaking, quasi-1/2 magnetization plateau and quasi spin polarized phases. A reentrant behavior of the staggered magnetization in a transverse field is observed. The quantum critical behaviors in longitudinal and transverse magnetic fields are disclosed to fall into different universality classes with corresponding conformal field central charge $c=1$ and 1/2, respectively. The experimental data of the compound NTENP under both longitudinal and transverse magnetic fields are nicely fitted, and the Luttinger liquid behavior of low temperature specific heat experimentally observed is also confirmed.

preprint2012arXiv

Thermoelectric Properties of Silicon Carbide Nanowires with Nitrogen Dopants and Vacancies

The thermoelectric properties of cubic zincblend silicon carbide nanowires (SiCNWs) with nitrogen impurities and vacancies along [111] direction are theoretically studied by means of atomistic simulations. It is found that the thermoelectric figure of merit ZT of SiCNWs can be significantly enhanced by doping N impurities together with making Si vacancies. Aiming at obtaining a large ZT, we study possible energetically stable configurations, and disclose that, when N dopants locate at the center, a small number of Si vacancies at corners are most favored for n-type nanowires, while a large number of Si vacancies spreading into the flat edge sites are most favored for p-type nanowires. For the SiCNW with a diameter of 1.1 nm and a length of 4.6 nm, the ZT value for the n-type is shown capable of reaching 1.78 at 900K. The conditions to get higher ZT values for longer SiCNWs are also addressed.

preprint2011arXiv

Global phase diagram of three-dimensional extended Boson Hubbard model - a continuous time Quantum Monte Carlo study

We present the global phase diagram of the extended boson Hubbard model on a simple cubic lattice by quantum Monte Carlo simulation with worm update algorithm. Four kinds of phases are supported by this model, including superfluid, supersolid, Mott, and charge density wave (CDW) states, which are identified in the phase diagram of chemical potential $μ$ versus nearest neighbor interaction V . By changing the chemical potential, a continuous transition is found from the Mott phase to a superfluid phase without breaking the translational symmetry. For an insulating CDW state, adding particles to it gives rise to a continuous transition to a supersolid phase, while removing particles usually leads to a first-order one to either supersolid or superfluid phase. By tuning the nearest neighbor interaction, one can realize the transition between two insulating phases, Mott and CDW with the same particle density, which turns out to be of the first-order. We also demonstrate that a supersolid phase with average particle density less than 1/2 can exist in a small region of $μ$ - V phase diagram.

preprint2011arXiv

Linearized Tensor Renormalization Group Algorithm for Thermodynamics of Quantum Lattice Models

A linearized tensor renormalization group (LTRG) algorithm is proposed to calculate the thermodynamic properties of one-dimensional quantum lattice models, that is incorporated with the infinite time-evolving block decimation technique, and allows for treating directly the two-dimensional transfer-matrix tensor network. To illustrate its feasibility, the thermodynamic quantities of the quantum XY spin chain are calculated accurately by the LTRG, and the precision is shown to be comparable with (even better than) the transfer matrix renormalization group (TMRG) method. Unlike the TMRG scheme that can only deal with the infinite chains, the present LTRG algorithm could treat both finite and infinite systems, and may be readily extended to boson and fermion quantum lattice models.

preprint2011arXiv

Low-energy properties of anisotropic two-dimensional spin-1/2 Heisenberg models in staggered magnetic fields

We present a systematic study of the anisotropic spin-1/2 Heisenberg model in staggered magnetic fields in two dimensions (2D). To mimic real materials, we consider a system of coupled, antiferromagnetic chains, whose interchain interaction can be either ferro- or antiferromagnetic. When the staggered field is commensurate with the magnetic interactions, an energy gap opens immediately and follows a power law as a function of the applied field, similar to the situation in 1D. When the field competes with the interactions, a quantum phase transition (QPT) occurs from a gapless, magnetically ordered phase at low fields to a gapped, disordered regime. We use a continuous-time Monte Carlo method to compute the staggered moment of the ordered phases and the spin gap of the disordered phases. We deduce the phase diagrams as functions of the anisotropy ratio and the applied field, and calculate the scaling behavior of the models in both quantities. We show that in the competitive case, the staggered field acts to maintain a regime of quasi-1D behavior around the QPT, and we discuss as a consequence the nature of the crossover from 1D to 2D physics.

preprint2011arXiv

Modified Tucker Decomposition for Tensor Network and Fast Linearized Tensor Renormalization Group Algorithm for Two-Dimensional Quantum Spin Lattice Systems

We propose a novel algorithm with a modified Tucker decomposition for tensor network that allows for efficiently and precisely calculating the ground state and thermodynamic properties of two-dimensional (2D) quantum spin lattice systems, and is coined as the fast linearized tensor renormalization group (fLTRG). Its amazing efficiency and precision are examined by studying the spin-1/2 anisotropic Heisenberg antiferromagnet on a honeycomb lattice, and the results are found to be fairly in agreement with the quantum Monte Carlo calculations. It is also successfully applied to tackle a quasi-2D spin-1/2 frustrated bilayer honeycomb Heisenberg model, where a quantum phase transition from an ordered antiferromagnetic state to a gapless quantum spin liquid phase is found. The thermodynamic behaviors of this frustrated spin system are also explored. The present fLTRG algorithm could be readily extended to other quantum lattice systems.

preprint2011arXiv

T-Carbon: A Novel Carbon Allotrope

A structurally stable crystalline carbon allotrope is predicted by means of the first-principles calculations. This allotrope can be derived by substituting each atom in diamond with a carbon tetrahedron, and possesses the same space group Fd^1 3m as diamond, which is thus coined as T- carbon. The calculations on geometrical, vibrational and electronic properties reveal that T-carbon, with a considerable structural stability and a much lower density 1.50 g/cm3, is a semiconductor with a direct band gap about 3.0 eV, and has a Vickers hardness 61.1 GPa lower than diamond but comparable with cubic boron nitride. Such a form of carbon, once obtained, would have wide applications in photocatalysis, adsoption, hydrogen storage and aerospace materials.

preprint2010arXiv

Magnetism and thermodynamics of spin-(1/2,1) decorated Heisenberg chain with spin-1 pendants

The magnetic and thermodynamic properties of a new ferrimagnetic decorated spin-(1/2,1) Heisenberg chain with spin-1 pendant spins are investigated for three cases: (A) J1,J2>0; (B) J1>0, J2<0; and (C) J1<0, J2>0, where J1 and J2 are the exchange couplings between spins in the chain and along the rung, respectively. The low-lying and magnetic properties are explored jointly by the real-space renormalization group, spin wave, and density-matrix renormalization group methods, while the transfer-matrix renormalization group method is invoked to study the thermodynamics. It is found that the magnon spectra consist of a gapless and two gapped branches. Two branches in case (C) have intersections. The coupling dependence of low-energy gaps are analyzed. In a magnetic field, an m=3/2 (m is the magnetization per unit cell) plateau is observed for case (A), while two plateaux at m=1/2 and 3/2 are observed for cases (B) and (C). Between the two plateaux in cases (B) and (C), the sublattice magnetizations for the spins coupled by ferromagnetic interactions have novel decreasing regions with increasing the magnetic field. At finite temperature, the zero-field susceptibility temperature product chi*T and specific heat exhibit distinct exotic features with varying the couplings and temperature for different cases. chi*T is found to converge as T approaches zero, which is different from the divergent behavior in the spin-(1/2,1) mixed-spin chain without pendants. The observed thermodynamic behaviors are also discussed with the help of their low-lying excitations.

preprint2010arXiv

Magnetization plateau and incommensurate spin modulation in Ca3Co2O6

The magnetic properties of a trigonal prism unit of the spin-2 frustrated compound Ca3Co2O6 are studied by means of the density-matrix renormalization group method. A magnetization plateau at $ms/3$ ($ms$ is the saturation magnetization) with ferrimagnetic structure is observed. By fitting the experimental data of magnetic curve, an estimation of the couplings gives J1=-26.84K, J_{2}=0.39K, and J_{3}=0.52K. The local magnetic moments are unveiled to exhibit an incommensurate sinusoidally modulation along the three chains of the trigonal prism, which gives a strong theoretical support to the experimentally observed incommensurate partially disordered antiferromagnetic state for Ca3Co2O6. The present result suggests that the modulation indeed originates from the competition of antiferromagnetic and ferromagnetic couplings.

preprint2010arXiv

Phase transitions and thermodynamics of the two-dimensional Ising model on a distorted Kagomé lattice

The two-dimensional Ising model on a distorted Kagomé lattice is studied by means of exact solutions and the tensor renormalisation group (TRG) method. The zero-field phase diagrams are obtained, where three phases such as ferromagnetic, ferrimagnetic and paramagnetic phases, along with the second-order phase transitions, have been identified. The TRG results are quite accurate and reliable in comparison to the exact solutions. In a magnetic field, the magnetization ($m$), susceptibility and specific heat are studied by the TRG algorithm, where the $m=1/3$ plateaux are observed in the magnetization curves for some couplings. The experimental data of susceptibility for the complex Co(N$_3$)$_2$(bpg)$\cdot$ DMF$_{4/3}$ are fitted with the TRG results, giving the couplings of the complex $J=22K$ and $J'=33K$.

preprint2010arXiv

Spin transfer in a ferromagnet-quantum dot and tunnel barrier coupled Aharonov-Bohm ring system with Rashba spin-orbit interactions

The spin transfer effect in ferromagnet-quantum dot (insulator)-ferromagnet Aharonov-Bohm (AB) ring system with Rashba spin-orbit (SO) interactions is investigated by means of Keldysh nonequilibrium Green function method. It is found that both the magnitude and direction of the spin transfer torque (STT) acting on the right ferromagnet electrode can be effectively controlled by changing the magnetic flux threading the AB ring or the gate voltage on the quantum dot. The STT can be greatly augmented by matching a proper magnetic flux and an SO interaction at a cost of low electrical current. The STT, electrical current, and spin current are uncovered to oscillate with the magnetic flux. The present results are expected to be useful for information storage in nanospintronics.

preprint2002arXiv

Effect of Spin-Flip Scattering on Electrical Transport in Magnetic Tunnel Junctions

By means of the nonequilibrium Green function technique, the effect of spin-flip scatterings on the spin-dependent electrical transport in ferromagnet-insulator-ferromagnet (FM-I-FM) tunnel junctions is investigated. It is shown that Julliere's formula for the tunnel conductance must be modified when including the contribution from the spin-flip scatterings. It is found that the spin-flip scatterings could lead to an angular shift of the tunnel conductance, giving rise to the junction resistance not being the largest when the orientations of magnetizations in the two FM electrodes are antiparallel, which may offer an alternative explanation for such a phenomenon observed previously in experiments in some FM-I-FM junctions. The spin-flip assisted tunneling is also observed.

preprint1995arXiv

Classical and Quantum Chaotic Behaviors of Two Colliding Harmonic Oscillators

We have systematically studied both classical and quantum chaotic behaviors of two colliding harmonic oscillators. The classical case falls in Kolmogorov-Arnold-Moser class. It is shown that there exists an energy threshold, above which the system becomes nonintegrable. For some values of the initial energy near the threshold, we have found that the ratio of frequencies of the two oscillators affects the Poincaré sections significantly. The largest Lyapunov character exponent depends linearly on the ratio of frequencies of the two oscillators away from the energy threshold in some chaotic regions, which shows that the chaotic behaviors of the system are mainly determined by the ratio. In the quantum case, for certain parameters, the distribution of the energy level spacings also varies with the ratio of frequencies of the two oscillators. The relation between the energy spectra and the ratio of frequencies of the two oscillators, the interaction constant, and the semi-classical quantization constant, is also investigated respectively.