Source author record

Manoranjan Kumar

Manoranjan Kumar 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

34works
11topics
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

34 published item(s)

preprint2022arXiv

Anomalous transport in itinerant van der Waals ferromagnets Fe$_n$GeTe$_2$ (\emph{n}=3, 4, 5)

Ferromagnetic (FM) semimetals Fe$_n$GeTe$_2$(n=3, 4, 5), exhibit several symmetry-protected band-crossing points or lines near the Fermi energy (E$_F$) and these topological properties of energy bands lead to interesting transport properties. We study these materials employing the first-principle calculations and the tight-binding Hamiltonian constructed by fitting the parameters of the first principles calculation. In the presence of spin-orbit coupling (SOC) for n=3,5 a large Berry curvature (BC) concentrated on the nodal lines is observed. The consequence of the correlation of the topological nodal line and magnetic moments on anomalous Hall conductivity (AHC) $σ_{xy}$ and anomalous Nernst conductivity (ANC) $α_{xy}$ have been investigated. We find $σ_{xy}=150$ S/cm for n=3, 295 S/cm for n=4, and 90 S/cm for n=5 at 0 K, while the ANC is observed as $α_{xy}=0.55$ A/Km for n=3, 0.10 A/Km for n=5, and 0.80 A/Km for n=4, at the E$_F$ at room temperature. Our calculated AHC values at 0 K, i.e., 150 S/cm for Fe$_3$GeTe$_2$ and 90 S/cm Fe$_5$GeTe$_2$, are consistent with the experimentally reported values. Also the experimentally reported value of ANC for Fe$_5$GeTe$_2$ is close to our calculated value at room temperature, i.e., 0.10 A/Km.

preprint2022arXiv

Band splitting induced Berry flux and intrinsic anomalous Hall conductivity in NiCoMnGa quaternary Heusler compound

The anomalous transport properties of Heusler compounds become a hotspot of research in recent years due to their unique band structure and possible application in spintronics. In this paper, we report the anomalous Hall effect in polycrystalline NiCoMnGa quaternary Heusler compound by experimental means and theoretical calculations. The experimental anomalous Hall conductivity (AHC) was found at about 256 S/cm at 10K with an intrinsic contribution of ~ 121 S/cm. The analysis of Hall data reveals the presence of both extrinsic and intrinsic contributions in AHE. Our theoretical calculations show that a pair of spin-orbit coupled band formed by the band splitting due to spin-orbit interaction (SOI) at the Fermi level produces a finite Berry flux in the system that provides the intrinsic AHC about 100 S/cm, which is in good agreement with the experiment.

preprint2022arXiv

Density matrix renormalization group approach to the low temperature thermodynamics of correlated 1D fermionic models

The low temperature thermodynamics of correlated 1D fermionic models with spin and charge degrees of freedom is obtained by exact diagonalization (ED) of small systems and followed by density matrix renormalization group (DMRG) calculations that target the lowest hundreds of states $\{E(N)\}$ at system size $N$ instead of the ground state. Progressively larger $N$ reaches $T < 0.05t$ in correlated models with electron transfer $t$ between first neighbors and bandwidth $4t$. The size dependence of the many-fermion basis is explicitly included for arbitrary interactions by scaling the partition function. The remaining size dependence is then entirely due to the energy spectrum $\{E(N)\}$ of the model. The ED/DMRG method is applied to Hubbard and extended Hubbard models, both gapped and gapless, with $N_e = N$ or $N/2$ electrons and is validated against exact results for the magnetic susceptibility $χ(T)$ and entropy $S(T)$ per site. Some limitations of the method are noted. Special attention is given to the bond-order-wave phase of the extended Hubbard model with competing interactions and low $T$ thermodynamics sensitive to small gaps.

preprint2022arXiv

Machine learning approach to study quantum phases of a frustrated one dimensional spin-1/2 system

Frustration driven quantum fluctuation leads to many exotic phases in the ground state and study of these quantum phase transitions is one of the most challenging areas of research in condensed matter physics. Here, a frustrated Heisenberg $J_1-J_2$ model of spin-1/2 chain with nearest exchange interaction $J_1$ and next nearest exchange interaction $J_2$ is studied using the principal component analysis (PCA) which is an unsupervised machine learning technique. In this method most probable spin configurations (MPSC) of ground-state (GS) and first excited state (FES) for different $J_2/J_1$ are used as the input in PCA to construct the co-variance matrix. The `quantified principal component' of the largest eigenvalue of co-variance matrix $p_1(J_2/J_1)$ is calculated and it is shown that the nature and variation of $p_1(J_2/J_1)$ can accurately predict the phase transitions and degeneracies in the GS. The $p_1(J_2/J_1)$ calculated from the MPSC of GS can only exhibit the signature of degeneracies in the GS, whereas, $p_1(J_2/J_1)$ calculated from MPSC of FES captures the gapless spin liquid (GSL)-dimer phase transition as well as all the degeneracies of the model system. We show that jump in $p_1(J_2/J_1)$ of FES at $J_2/J_1 \approx 0.241$, indicates the GSL-dimer phase transition, whereas its kinks give the signature of the GS degeneracies. The scatter plot of first two principal components of FES shows distinct band formation for different phases.

preprint2022arXiv

Spin Peierls transition of $J_{1}-J_{2}$ and extended models with ferromagnetic $J_{1}$.Sublattice dimerization and thermodynamics of zigzag chains in $β$-TeVO$_{4}$

The spin$-1/2$ chain with ferromagnetic exchange $J_1 < 0$ between first neighbors and antiferromagnetic $J_2 > 0$ between second neighbors supports two spin-Peierls (SP) instabilities depending on the frustration $α= J_2/\vert J_1\vert$. Instead of chain dimerization with two spins per unit cell, $J_1-J_2$ models with $α> 0.65$ and linear spin-phonon coupling are unconditionally unstable to sublattice dimerization with four spins per unit cell. Unequal $J_1$ to neighbors to the right and left extends the model to gapped ($γ> 0$) chains with conditional SP transitions at $T_{SP}$ to dimerized sublattices and a weaker specific heat $C(T)$ anomaly. The spin susceptibility $χ(T)$ and $C(T)$ are obtained in the thermodynamic limit by a combination of exact diagonalization of small systems with $α> 0.65$ and density matrix renormalization group (DMRG) calculations of systems up to $N \sim 100$ spins. Both $J_1-J_2$ and $γ> 0$ models account quantitatively for $χ(T)$ and $C(T)$ in the paramagnetic phase of $β$-TeVO$_{4}$ for $T > 8$ K, but lower $T$ indicates a gapped chain instead of a $J_1-J_2$ model as previously thought. The same parameters and $T_{SP} = 4.6$ K generate a $C(T)/T$ anomaly that reproduces the anomaly at the $4.6$ K transition of $β$-TeVO$_{4}$, but not the weak $χ(T)$ signature.

preprint2022arXiv

Study of the Berezinskii-Kosterlitz-Thouless transition: An unsupervised machine learning approach

The Berezinskii-Kosterlitz-Thouless (BKT) transition in magnetic system is an intriguing phenomena and an accurate estimation of the BKT transition temperature has been a long-standing problem. In this work we explore the anisotropic classical Heisenberg XY and XXZ models with ferromagnetic exchange on a square lattice and antiferromagnetic exchange on a triangular lattice using an unsupervised machine learning approach called principal component analysis (PCA). In earlier studies of the BKT transition, spin configurations and vorticities calculated from Monte Carlo method are used to determine the transition temperature $T_{BKT}$, but those methods fail to give any conclusive results by analyzing the principal components in the PCA approach. In this work vorticities are used as initial input to the PCA and curve of the first principal component with temperature is fitted with a function to determine an accurate value of $T_{BKT}$. This procedure works well for anisotropic classical Heisenberg with ferromagnetic exchange on square lattice as well as for frustrated antiferromagnetic exchange on a triangular lattice. The classical anisotropic Heisenberg antiferromagnetic model on the triangular lattice has two close transitions; the BKT at $T_{BKT}$ and Ising like phase transition for chirality at $T_c$ and it is difficult to separate these transition points. It is also noted that using the PCA method and manipulation of their first principal component, not only separation of transition points are possible but also transition temperature can be determined accurately.

preprint2021arXiv

Topological transitions to Weyl states in bulk Bi$_2$Se$_3$: Effect of hydrostatic pressure and doping

Bi$_2$Se$_3$, a layered three dimensional (3D) material, exhibits topological insulating properties due to presence of surface states and a band gap of 0.3 eV in the bulk. We study the effect hydrostatic pressure $P$ and doping with rare earth elements on the topological aspect of this material in bulk from a first principles perspective. Our study shows that under a moderate pressure of P$>$7.9 GPa, the bulk electronic properties show a transition from an insulating to a Weyl semi-metal state due to band inversion. This electronic topological transition may be correlated to a structural change from a layered van der Waals material to a 3D system observed at $P$=7.9 GPa. At large $P$ density of states have significant value at the Fermi-energy. Intercalating Gd with a small doping fraction between Bi$_2$Se$_3$ layers drives the system to a metallic anti-ferromagnetic state, with Weyl nodes below the Fermi-energy. At the Weyl nodes time reversal symmetry is broken due to finite local field induced by large magnetic moments on Gd atoms. However, substituting Bi with Gd induces anti-ferromagnetic order with an increased direct band gap. Our study provides novel approaches to tune topological transitions, particularly in capturing the elusive Weyl semimetal states, in 3D topological materials.

preprint2020arXiv

Fermion parity gap and exponential ground state degeneracy of the one-dimensional Fermi gas with intrinsic attractive interaction

We examine the properties of a one-dimensional (1D) Fermi gas with attractive intrinsic (Hubbard) interactions in the presence of spin-orbit coupling and Zeeman field by numerically computing the pair binding energy, excitation gap, and susceptibility to local perturbations using the density matrix renormalization group. Such a system can, in principle, be realized in a system of ultracold atoms confined in a 1D optical lattice. We note that, in the presence of spatial interfaces introduced by a smooth parabolic potential, the pair binding and excitation energy of the system decays exponentially with the system size, pointing to the existence of an exponential ground state degeneracy, and is consistent with recent works. However, the susceptibility of the ground state degeneracy of this number-conserving system to local impurities indicates that the energy gap vanishes as a power law with the system size in the presence of local perturbations. We compare this system with the more familiar system of an Ising antiferromagnet in the presence of a transverse field realized with Rydberg atoms and argue that the exponential splitting in the clean number-conserving 1D Fermi system is similar to a phase with only conventional order.

preprint2020arXiv

Haldane and Dimer phases in a frustrated spin chain: an exact groundstate and associated topological phase transition

A Heisenberg spin-$s$ chain with alternating ferromagnetic ($-J_1^F<0$) and antiferromagnetic ($J_1^A>0$) nearest-neighbor (NN) interactions, exhibits the Dimer and spin-$2s$ Haldane phases in the limits $J_1^F/J_1^A \rightarrow 0$ and $J_1^F/J_1^A \rightarrow \infty$ respectively. These two phases are understood to be topologically equivalent. Induction of the frustration through the next nearest-neighbor ferromagnetic interaction ($-J_2^F<0$) produces a very rich quantum phase diagram. With frustration, the whole phase diagram is divided into a ferromagnetic (FM) and a nonmagnetic (NM) phase. For $s=1/2$, the full NM phase is seen to be of Haldane-Dimer type, but for $s>1/2$, a spiral phase comes between the FM and the Haldane-Dimer phases. The study of a suitably defined string-order parameter and spin-gap at the phase boundary indicates that the Haldane-Dimer and spiral phases have different topological characters. We also find that, along the $J_2^F=\frac 12 J_1^F$ line in the NM phase, an NN dimer state is the {\it exact} groundstate, provided $J_1^A>J_C=κJ_1^F$ where $κ\le s + h$ for applied magnetic field $h$. Without magnetic field, the position of $J_C$ is on the FM-NM phase boundary when $s=1/2$, but for $s>1/2$, the location of $J_C$ is on the phase separation line between the Haldane-Dimer and spiral phases.

preprint2020arXiv

Magnetization plateaus of spin-$\mathbf{\frac{1}{2}}$ system on a 5/7 skewed ladder

Magnetization plateaus are some of the most striking manifestations of frustration in low-dimensional spin systems. We present numerical studies of magnetization plateaus in the fascinating spin-1/2 skewed ladder system obtained by alternately fusing five- and seven-membered rings. This system exhibits three significant plateaus at $m = 1/4$, $1/2$ and $3/4$, consistent with the Oshikawa-Yamanaka-Affleck condition. Our numerical as well as perturbative analysis shows that the ground state can be approximated by three weakly coupled singlet dimers and two free spins, in the absence of a magnetic field. With increasing applied magnetic field, the dimers progressively become triplets with large energy gaps to excited states, giving rise to stable magnetization plateaus. Finite-temperature studies show that $m=1/4$ and $1/2$ plateaus are robust and survive thermal fluctuations while the $m=3/4$ plateau shrinks rapidly due to thermal noise. The cusps at the ends of a plateau follow the algebraic square-root dependence on $B$.

preprint2020arXiv

Modeling the spin-Peierls transition of spin-$1/2$ chains with correlated states: $J_1-J_2$ model, CuGeO$_3$ and TTF-CuS$_4$C$_4$(CF$_3$)$_4$

The spin-Peierls transition at $T_{SP}$ of spin-$1/2$ chains with isotropic exchange interactions has previously been modeled as correlated for $T > T_{SP}$ and mean field for $T < T_{SP}$. We use correlated states throughout in the $J_1-J_2$ model with antiferromagnetic exchange $J_1$ and $J_2 = αJ_1$ between first and second neighbors, respectively, and variable frustration $0 \leq α\leq 0.50$. The thermodynamic limit is reached at high $T$ by exact diagonalization of short chains and at low $T$ by density matrix renormalization group calculations of progressively longer chains. In contrast to mean field results, correlated states of 1D models with linear spin-phonon coupling and a harmonic adiabatic lattice provide an internally consistent description in which the parameter $T_{SP}$ yields both the stiffness and the lattice dimerization $δ(T)$. The relation between $T_{SP}$ and $Δ(δ,α)$, the $T = 0$ gap induced by dimerization, depends strongly on $α$ and deviates from the BCS gap relation that holds in uncorrelated spin chains. Correlated states account quantitatively for the magnetic susceptibility of TTF-CuS$_4$C$_4$(CF$_3$)$_4$ crystals ($J_1 = 79$ K, $α= 0$, $T_{SP} = 12$ K) and CuGeO$_3$ crystals ($J_1 = 160$ K, $α= 0.35$, $T_{SP} = 14$ K). The same parameters describe the specific heat anomaly of CuGeO$_3$ and inelastic neutron scattering. Modeling the spin-Peierls transition with correlated states exploits the fact that $δ(0)$ limits the range of spin correlations at $T = 0$ while $T > 0$ limits the range at $δ= 0$.

preprint2020arXiv

Nonquenched rotators ease flocking and memorise it

We introduce a minimal model for a two-dimensional polar flock with nonquenched rotators, and show that the rotators make the usual macroscopic long-range order of the flock more robust than the clean system. The rotators memorise the flock-information which helps in establishing the robustness. Moreover, the memory of the rotators assists in probing the moving flock. We also formulate a hydrodynamic framework for the microscopic model that makes our study comprehensive. Using linearised hydrodynamics, it is shown that the presence of such nonquenched heterogeneities increases the sound speeds of the flock. The enhanced sound speeds lead to faster convection of information and consequently the robust ordering in the system. We argue that similar nonquenched heterogeneities may be useful in monitoring and controlling large crowds.

preprint2020arXiv

Speed inhomogeneity accelerates the information transfer in polar flock

A collection of self-propelled particles (SPPs) shows coherent motion and exhibits a true long range ordered (LRO) state in two dimensions. Various studies show that the presence of spatial inhomogeneities can destroy the usual long-range ordering in the system. However, the effects of inhomogeneity due to the intrinsic properties of the particles are barely addressed. In this paper we consider a collection of polar SPPs moving with inhomogeneous speed (IS) on a two dimensional substrate, which can arise due to varying physical strength of the individual particle. To our surprise, the IS not only preserves the usual long-range ordering present in the homogeneous speed models but also induces faster ordering in the system. Furthermore, The response of the flock to an external perturbation is also faster, compared to Vicsek like model systems, due to the frequent update of neighbors of each SPP in the presence of the IS. Therefore, our study shows that the IS can help in faster information transfer in the moving flock.

preprint2020arXiv

Tunneling density of states in a Y junction of Tomonaga-Luttinger liquid wires: A density matrix renormalization group study

It is well known that the pristine bulk of an interacting one-dimensional system in Tomonaga-Luttinger liquid (TLL) phase shows power law suppression of quasi-particle tunneling amplitude for all values of TLL parameter $g$, in the zero energy limit. We perform a density matrix renormalization group (DMRG) study of a fully symmetric Y junction of TLL wires and observe an anomalous enhancement of the tunneling density of states (TDOS) in the vicinity of the junction for both (a) interacting bosons case and (b) interacting fermions case, when $g>1$. We also observe suppression of TDOS for $g<1$ for both bosonic and fermionic cases. We find that the TDOS enhancements follow different power laws for bosonic and fermionic cases which suggests that these represent distinct fixed points, owing to statistical correlations which play an important role at the Y junction. Analysis of static conductance for the junction indicates that the fixed point for $1<g<3$ resembles the mysterious $M$ fixed point of Y junction predicted by Oshikawa, Chamon, and Affleck [J. Stat. Mech. P02008 (2006)]. We also show that the TDOS enhancement spans over a length scale of $\propto ω^{-1}$ from the junction, for $1<g<3$.

preprint2019arXiv

Melting of Ferromagnetic Order on a Trellis Ladder

The ground state properties of a frustrated spin-1/2 system is studied on a trellis ladder which is composed of two zigzag ladders interacting through rung interactions. The presence of rung interaction between the zigzag ladders induces a non-magnetic ground state, although, each of zigzag ladders has ferromagnetic order in weak antiferromagnetic leg interaction limit. The rung interaction also generates rung dimers and opens spin gap which increases rapidly with rung interaction strength. The correlation between spins decreases exponentially with the distance between them.

preprint2016arXiv

An Efficient Density Matrix Renormalization Group Algorithm for Chains with Periodic Boundary Condition

The Density Matrix Renormalization Group (DMRG) is a state-of-the-art numerical technique for a one dimensional quantum many-body system; but calculating accurate results for a system with Periodic Boundary Condition (PBC) from the conventional DMRG has been a challenging job from the inception of DMRG. The recent development of the Matrix Product State (MPS) algorithm gives a new approach to find accurate results for the one dimensional PBC system. The most efficient implementation of the MPS algorithm can scale as O($p \times m^3$), where $p$ can vary from 4 to $m^2$. In this paper, we propose a new DMRG algorithm, which is very similar to the conventional DMRG and gives comparable accuracy to that of MPS. The computation effort of the new algorithm goes as O($m^3$) and the conventional DMRG code can be easily modified for the new algorithm.

preprint2016arXiv

Boundary-induced spin density waves in linear Heisenberg antiferromagnetic spin chains with $\mathbf{S \ge 1}$

Linear Heisenberg antiferromagnets (HAFs) are chains of spin-$S$ sites with isotropic exchange $J$ between neighbors. Open and periodic boundary conditions return the same ground state energy in the thermodynamic limit, but not the same spin $S_G$ when $S \ge 1$. The ground state of open chains of N spins has $S_G = 0$ or $S$, respectively, for even or odd N. Density matrix renormalization group (DMRG) calculations with different algorithms for even and odd N are presented up to N = 500 for the energy and spin densities $ρ(r,N)$ of edge states in HAFs with $S = 1$, 3/2 and 2. The edge states are boundary-induced spin density waves (BI-SDWs) with $ρ(r,N)\propto(-1)^{r-1}$ for $r=1,2,\ldots N$. The SDWs are in phase when N is odd, out of phase when N is even, and have finite excitation energy $Γ(N)$ that decreases exponentially with N for integer $S$ and faster than 1/N for half integer $S$. The spin densities and excitation energy are quantitatively modeled for integer $S$ chains longer than $5 ξ$ spins by two parameters, the correlation length $ξ$ and the SDW amplitude, with $ξ= 6.048$ for $S = 1$ and 49.0 for $S = 2$. The BI-SDWs of $S = 3/2$ chains are not localized and are qualitatively different for even and odd N. Exchange between the ends for odd N is mediated by a delocalized effective spin in the middle that increases $|Γ(N)|$ and weakens the size dependence. The nonlinear sigma model (NL$σ$M) has been applied the HAFs, primarily to $S = 1$ with even N, to discuss spin densities and exchange between localized states at the ends as $Γ(N) \propto (-1)^N \exp(-N/ξ)$...

preprint2016arXiv

Numerical study of incommensurate and decoupled phases of spin-1/2 chains with isotropic exchange J1, J2 between first and second neighbors

The spin-1/2 chain with isotropic exchange J1, J2 > 0 between first and second neighbors is frustrated for either sign of J1 and has a singlet ground state (GS) for J1/J2 $\ge - 4$. Its rich quantum phase diagram supports gapless, gapped, commensurate (C), incommensurate (IC) and other phases. Critical points J1/J2 are evaluated using exact diagonalization (ED) and density matrix renormalization group (DMRG) calculations. The wave vector $q_G$ of spin correlations is related to GS degeneracy and obtained as the peak of the spin structure factor $S(q)$. Variable $q_G$ indicates IC phases in two $J1/J2$ intervals, [ -4, -1.24] and [0.44, 2], and a C-IC point at J1/J2 = 2. The decoupled C phase in [-1.24, 0.44] has constant $q_G = π/2$, nondegenerate GS, and a lowest triplet state with broken spin density on sublattices of odd and even numbered sites. The lowest triplet and singlet excitations, $E_m$ and $E_σ$, are degenerate in finite systems at specific frustration $J1/J2$. Level crossing extrapolates in the thermodynamic limit to the same critical points as $q_G$. The $S(q)$ peak diverges at $q_G = π$ in the gapless phase with $J1/J2 > 4.148$ and quasi-long-range order (QLRO({$π$})). $S(q)$ diverges at $\pm π/2$ in the decoupled phase with QLRO({$π$}/2), but is finite in gapped phases with finite range correlations. Numerical results and field theory agree at small $J2/J1$ but disagree for the decoupled phase with weak exchange $J1$ between sublattices. Two related models are summarized: one has an exact gapless decoupled phase with QLRO({$π$}/2) and no IC phases; the other has a single IC phase without a decoupled phase in between.

preprint2015arXiv

Anderson transitions in disordered two-dimensional lattices

We numerically analyze the energy level statistics of the Anderson model with Gaussian site disorder and constant hopping. The model is realized on different two-dimensional lattices, namely, the honeycomb, the kagomé, the square, and the triangular lattice. By calculating the well-known statistical measures viz., nearest neighbor spacing distribution, number variance, the partition number and the dc electrical conductivity from Kubo-Greenwood formula, we show that there is clearly a delocalization to localization transition with increasing disorder. Though the statistics in different lattice systems differs when compared with respect to the change in the disorder strength only, we find there exists a single complexity parameter, a function of the disorder strength, coordination number, localization length, and the local mean level spacing, in terms of which the statistics of the fluctuations matches for all lattice systems at least when the Fermi energy is selected from the bulk of the energy levels.

preprint2015arXiv

Density Induced Phases in Active Nematic

We introduce a minimal model for a collection of self-propelled apolar active particles, also called as `active nematic', on a two-dimensional substrate and study the order-disorder transition with the variation of density. The particles interact with their neighbours within the framework of the Lebwohl-Lasher model and move asymmetrically, along their orientation, to unoccupied nearest neighbour lattice sites. At a density lower than the equilibrium isotropic-nematic transition density, the active nematic shows a first order transition from the isotropic state to a banded state. The banded state extends over a range of density, and the scalar order parameter of the system shows a plateau like behaviour, similar to that of the magnetic systems. In the large density limit the active nematic shows a bistable behaviour between a homogeneous ordered state with global ordering and an inhomogeneous mixed state with local ordering. The study of the above phases with density variation is scant and gives significant insight of complex behaviours of many biological systems.

preprint2015arXiv

Efficient Density Matrix Renormalization Group algorithm to study Y-Junctions with integer and half-integer spin

An efficient density matrix renormalization group (DMRG) algorithm is presented and applied to Y-junctions, systems with three arms of $n$ sites that meet at a central site. The accuracy is comparable to DMRG of chains. As in chains, new sites are always bonded to the most recently added sites and the superblock Hamiltonian contains only new or once renormalized operators. Junctions of up to $N = 3n + 1 \approx 500$ sites are studied with antiferromagnetic (AF) Heisenberg exchange $J$ between nearest-neighbor spins $S$ or electron transfer $t$ between nearest neighbors in half-filled Hubbard models. Exchange or electron transfer is exclusively between sites in two sublattices with $N_A \ne N_B$. The ground state (GS) and spin densities $ ρ_r = <S^z_r >$ at site $r$ are quite different for junctions with $S$ = 1/2, 1, 3/2 and 2. The GS has finite total spin $S_G = 2S (S)$ for even (odd) $N$ and for $M_G =S_G$ in the $S_G$ spin manifold, $ρ_r > 0 (< 0)$ at sites of the larger (smaller) sublattice. $S$ = 1/2 junctions have delocalized states and decreasing spin densities with increasing $N$. $S$ = 1 junctions have four localized $S_z = 1/2$ states at the end of each arm and centered on the junction, consistent with localized states in $S$ = 1 chains with finite Haldane gap. The GS of $S$ = 3/2 or 2 junctions of up to 500 spins is a spin density wave (SDW) with increased amplitude at the ends of arms or near the junction. Quantum fluctuations completely suppress AF order in $S$ = 1/2 or 1 junctions, as well as in half-filled Hubbard junctions, but reduce rather than suppress AF order in $S$ = 3/2 or 2 junctions.

preprint2015arXiv

Level crossing, spin structure factor and quantum phases of the frustrated spin-1/2 chain with first and second neighbor exchange

The spin-1/2 chain with isotropic Heisenberg exchange $J_1$, $J_2 > 0$ between first and second neighbors is frustrated for either sign of J1. Its quantum phase diagram has critical points at fixed $J_1/J_2$ between gapless phases with nondegenerate ground state (GS) and quasi-long-range order (QLRO) and gapped phases with doubly degenerate GS and spin correlation functions of finite range. In finite chains, exact diagonalization (ED) estimates critical points as level crossing of excited states. GS spin correlations enter in the spin structure factor $S(q)$ that diverges at wave vector $q_m$ in QLRO($q_m$) phases with periodicity $2π/q_m$ but remains finite in gapped phases. $S(q_m)$ is evaluated using ED and density matrix renormalization group (DMRG) calculations. Level crossing and the magnitude of $S(q_m)$ are independent and complementary probes of quantum phases, based respectively on excited and ground states. Both indicate a gapless QLRO($π/2$) phase between $-1.2 < J_1/|J_2| < 0.45$. Numerical results and field theory agree well for quantum critical points at small frustration $J_2$ but disagree in the sector of weak exchange $J_1$ between Heisenberg antiferromagnetic chains on sublattices of odd and even-numbered sites.

preprint2014arXiv

Magnetic Field Induced Exotic Phases in Isotropic Frustrated Spin-1/2 chain

The frustrated isotropic $J_1-J_2$ model with ferromagnetic $J_1$ and anti-ferromagnetic $J_2$ interactions in presence of an axial magnetic field shows many exotic phases, such as vector chiral and multipolar phases. The existing studies of the phase boundaries of these systems are based on the indirect evidences such as correlation functions {\it etc}. In this paper, the phase boundaries of these exotic phases are calculated based on order parameters and jumps in the magnetization. In the strong magnetic field, $Z_2$ symmetry is broken, therefore, order parameter of the vector chiral phase is calculated using the broken symmetry states. Our results obtained using the modified density matrix renormalization group and exact diagonalization methods, suggest that the vector chiral phase exist only in narrow range of parameter space $J_2/J_1$.

preprint2014arXiv

Spin structure factor and quantum phases of frustrated spin-1/2 chains

The static structure factor S(q) of frustrated spin-1/2 chains with isotropic exchange and a singlet ground state (GS) diverges at wave vector q_m when the GS has quasi-long-range order (QLRO) with periodicity 2π/q_m but S(q_m) is finite in bond-order-wave (BOW) phases with finite-range spin correlations. Exact diagonalization and density matrix renormalization group (DMRG) calculations of S(q) indicate a decoupled phase with QLRO and q_m = π/2 in chains with large antiferromagnetic exchange between second neighbors. S(q_m) identifies quantum phase transitions based on GS spin correlations.

preprint2013arXiv

Decoupled phase of frustrated spin-1/2 antiferromagnetic chains with and without long range order in the ground state

The quantum phases of one-dimensional spin $s= 1/2$ chains are discussed for models with two parameters, frustrating exchange $g = J_2 > 0$ between second neighbors and normalized nonfrustrating power-law exchange with exponent $α$ and distance dependence $r^{-α}$. The ground state (GS) at $g = 0$ has long-range order (LRO) for $α< 2$, long-range spin fluctuations for $α> 2$. The models conserve total spin $S = S_A + S_B$, have singlet GS for any $g$, $α\ge 0$ and decouple at $1/g = 0$ to linear Heisenberg antiferromagnets on sublattices $A$ and $B$ of odd and even-numbered sites. Exact diagonalization of finite chains gives the sublattice spin $ \ < S^2_A \ >$, the magnetic gap $E_m$ to the lowest triplet state and the excitation $E_σ$ to the lowest singlet with opposite inversion symmetry to the GS. An analytical model that conserves sublattice spin has a first order quantum transition at $g_c = 1/4{\rm ln2}$ from a GS with perfect LRO to a decoupled phase with $S_A = S_B = 0$ for $g \ge 4/π^2$ and no correlation between spins in different sublattices. The model with $α= 1$ has a first order transition to a decoupled phase that closely resembles the analytical model. The bond order wave (BOW) phase and continuous quantum phase transitions of finite models with $α\ge 2$ are discussed in terms of GS degeneracy where $E_σ(g) = 0$, excited state degeneracy where $E_σ(g) = E_m(g)$, and $\ < S^2_A \ >$. The decoupled phase at large frustration has nondegenerate GS for any exponent $α$ and excited states related to sublattice excitations.

preprint2012arXiv

Spin parity and broken symmetry in finite spin-1/2 chains with frustrated exchange: quantum transition from high to low spin

Exact diagonalization of finite spin-1/2 chains with periodic boundary conditions is applied to the ground state (gs) of chains with ferromagnetic (F) exchange $J_1 < 0$between first neighbors, antiferromagnetic (AF) exchange $J_2 = αJ_1 > 0$between second neighbors, and axial anisotropy $0 \le Δ\le 1$. In zero field, the gs is in the $S_z = 0$ sector for the relevant parameters and is doubly degenerate at multiple points $γ_m = (α_m, Δ_m)$ in the $α$, $Δ$ plane. Degeneracy under inversion at sites or spin parity or both leads, respectively, to a bond order wave (BOW), to staggered magnetization or to vector chiral (VC) order. Exact results up to $N = 28$ spins directly yield order parameters and spin correlation functions whose weak N dependencies allow inferences about infinite chains. The high-spin gs at $J_2 = 0$ changes discontinuously at $γ_1 = (-1/4, 1)$ to a singlet in the isotropic ($Δ= 1$) chain. The transition from high to low spin $S(α, Δ)$ is continuous for $ Δ< Δ_B = 0.95 \pm 0.01$ on the degeneracy line $α_1(Δ)$. The gs has staggered magnetization between $Δ_A = 0.72$ and $Δ_B$, and a BOW for $Δ< Δ_A$. When both inversion and spin parity are reversed at $γ_m$, the correlation functions $C(p)$ for spins separated by $p$ sites are identical. $C(p)$ minima are shifted by $π/2$ from the minima of VC order parameters at separation $p$, consistent with right and left-handed helices along the z axis and spins in the xy plane. Degenerate gs of finite chains are related to quantum phase diagrams of extended $α$, $Δ$ chains, with good agreement for order parameters along the line $α_1(Δ)$.

preprint2011arXiv

A Density Matrix Renormalization Group Method Study of Optical Properties of Porphines and Metalloporphines

The symmetrized Density-Matrix-Renormalization-Group (DMRG) method is used to study linear and nonlinear optical properties of Free base porphine and metallo-porphine. Long-range interacting model, namely, Pariser-Parr-Pople (PPP) model is employed to capture the quantum many body effect in these systems. The non-linear optical coefficients are computed within correction vector method. The computed singlet and triplet low-lying excited state energies and their charge densities are in excellent agreement with experimental as well as many other theoretical results. The rearrangement of the charge density at carbon and nitrogen sites, on excitation, is discussed. From our bond order calculation, we conclude that porphine is well described by the 18-annulenic structure in the ground state and the molecule expands upon excitation. We have modelled the regular metalloporphine by taking an effective electric field due to the metal ion and computed the excitation spectrum. Metalloporphines have $D_{4h}$ symmetry and hence have more degenerate excited states. The ground state of Metalloporphines show 20-annulenic structure, as the charge on the metal ion increases. The linear polarizability seems to increase with the charge initially and then saturates. The same trend is observed in third order polarizability coefficients.

preprint2011arXiv

Bond order wave (BOW) phase of the extended Hubbard model: Electronic solitons, paramagnetism, coupling to Peierls and Holstein phonons

The bond order wave (BOW) phase of the extended Hubbard model (EHM) in one dimension (1D) is characterized at intermediate correlation $U = 4t$ by exact treatment of $N$-site systems. Linear coupling to lattice (Peierls) phonons and molecular (Holstein) vibrations are treated in the adiabatic approximation. The molar magnetic susceptibility $χ_M(T)$ is obtained directly up to $N = 10$. The goal is to find the consequences of a doubly degenerate ground state (gs) and finite magnetic gap $E_m$ in a regular array. Degenerate gs with broken inversion symmetry are constructed for finite $N$ for a range of $V$ near the charge density wave (CDW) boundary at $V \approx 2.18t$ where $E_m \approx 0.5t$ is large. The electronic amplitude $B(V)$ of the BOW in the regular array is shown to mimic a tight-binding band with small effective dimerization $δ_{eff}$. Electronic spin and charge solitons are elementary excitations of the BOW phase and also resemble topological solitons with small $δ_{eff}$. Strong infrared intensity of coupled molecular vibrations in dimerized 1D systems is shown to extend to the regular BOW phase, while its temperature dependence is related to spin solitons. The Peierls instability to dimerization has novel aspects for degenerate gs and substantial $E_m$ that suppresses thermal excitations. Finite $E_m$ implies exponentially small $χ_M(T)$ at low temperature followed by an almost linear increase with $T$. The EHM with $U = 4t$ is representative of intermediate correlations in quasi-1D systems such as conjugated polymers or organic ion-radical and charge-transfer salts. The vibronic and thermal properties of correlated models with BOW phases are needed to identify possible physical realizations.

preprint2011arXiv

Density matrix renormalization group algorithm for Bethe lattices of spin 1/2 or 1 sites with Heisenberg antiferromagnetic exchange

An efficient density matrix renormalization group (DMRG) algorithm is presented for the Bethe lattice with connectivity $Z = 3$ and antiferromagnetic exchange between nearest neighbor spins $s= 1/2$ or 1 sites in successive generations $g$. The algorithm is accurate for $s = 1$ sites. The ground states are magnetic with spin $S(g) = 2^g s$, staggered magnetization that persists for large $g > 20$ and short-range spin correlation functions that decrease exponentially. A finite energy gap to $S > S(g)$ leads to a magnetization plateau in the extended lattice. Closely similar DMRG results for $s$ = 1/2 and 1 are interpreted in terms of an analytical three-site model.

preprint2011arXiv

Magnetic susceptibility of alkali-TCNQ salts and extended Hubbard models with bond order and charge density wave phases

The molar spin susceptibilities $χ(T)$ of Na-TCNQ, K-TCNQ and Rb-TCNQ(II) are fit quantitatively to 450 K in terms of half-filled bands of three one-dimensional Hubbard models with extended interactions using exact results for finite systems. All three models have bond order wave (BOW) and charge density wave (CDW) phases with boundary $V = V_c(U)$ for nearest-neighbor interaction $V$ and on-site repulsion $U$. At high $T$, all three salts have regular stacks of $\rm TCNQ^-$ anion radicals. The $χ(T)$ fits place Na and K in the CDW phase and Rb(II) in the BOW phase with $V \approx V_c$. The Na and K salts have dimerized stacks at $T < T_d$ while Rb(II) has regular stacks at 100K. The $χ(T)$ analysis extends to dimerized stacks and to dimerization fluctuations in Rb(II). The three models yield consistent values of $U$, $V$ and transfer integrals $t$ for closely related $\rm TCNQ^-$ stacks. Model parameters based on $χ(T)$ are smaller than those from optical data that in turn are considerably reduced by electronic polarization from quantum chemical calculation of $U$, $V$ and $t$ on adjacent $\rm TCNQ^-$ ions. The $χ(T)$ analysis shows that fully relaxed states have reduced model parameters compared to optical or vibration spectra of dimerized or regular $\rm TCNQ^-$ stacks.

preprint2010arXiv

Bond-order wave phase, spin solitons and thermodynamics of a frustrated linear spin-1/2 Heisenberg antiferromagnet

The linear spin-1/2 Heisenberg antiferromagnet with exchanges $J_1$, $J_2$ between first and second neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at $J_2/J_1 = 0.2411$ and is particularly simple at $J_2/J_1 = 1/2$. The BOW phase has a doubly degenerate singlet ground state, broken inversion symmetry and a finite energy gap $E_m$ to the lowest triplet state. The interval $0.4<J_2/J_1<1.0$ has large $E_m$ and small finite size corrections. Exact solutions are presented up to $N=28$ spins with either periodic or open boundary conditions and for thermodynamics up to $N=18$. The elementary excitations of the BOW phase with large $E_m$ are topological spin-1/2 solitons that separate BOWs with opposite phase in a regular array of spins. The molar spin susceptibility $χ_M(T)$ is exponentially small for $T \ll E_m$ and increases nearly linearly with $T$ to a broad maximum. $J_1$, $J_2$ spin chains approximate the magnetic properties of the BOW phase of Hubbard-type models and provide a starting point for modeling alkali-TCNQ salts.

preprint2010arXiv

Ionization potentials of crystalline organic thin films: Position dependence due to molecular shape and charge redistribution

In addition to electronic polarization or charge redistribution, the shape of neutral conjugated molecules yields position-dependent ionization potentials and electron affinities in organic thin films. Self-consistent I(n) and A(n) are computed in each layer n of 10-layer films of prototypical organics on a metal. The depth dependence of I(n) is discussed at surfaces of anthracene, C60 and PTCDA. The shape contribution can be substantial, up to 0.5 eV, and comes primarily from charge-quadrupole interactions.

preprint2010arXiv

Modified DMRG algorithm for the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange: Comparison with field theory at large $J_2/J_1$

A modified density matrix renormalization group (DMRG) algorithm is applied to the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange $J_1$, $J_2$ between first and second neighbors. The modified algorithm yields accurate results up to $J_2/J_1 \approx 4$ for the magnetic gap $Δ$ to the lowest triplet state, the amplitude $B$ of the bond order wave (BOW) phase, the wavelength $λ$ of the spiral phase, and the spin correlation length $ξ$. The $J_2/J_1$ dependences of $Δ$, $B$, $λ$ and $ξ$ provide multiple comparisons to field theories of the zigzag chain. The twist angle of the spiral phase and the spin structure factor yield additional comparisons between DMRG and field theory. Attention is given to the numerical accuracy required to obtain exponentially small gaps or exponentially long correlations near a quantum phase transition.

preprint2007arXiv

DMRG study of scaling exponents in spin-1/2 Heisenberg chains with dimerization and frustration

In conformal field theory, key properties of spin-1/2 chains, such as the ground state energy per site and the excitation gap scale with dimerization delta as delta^alpha with known exponents alpha and logarithmic corrections. The logarithmic corrections vanish in a spin chain with nearest (J=1) and next nearest neighbor interactions (J_2), for J_2c=0.2411. DMRG analysis of a frustrated spin chain with no logarithmic corrections yields the field theoretic values of alpha, and the scaling relation is valid up to the physically realized range, delta ~ 0.1. However, chains with logarithmic corrections (J_2<0.2411 J) are more accurately fit by simple power laws with different exponents for physically realized dimerizations. We show the exponents decreasing from approximately 3/4 to 2/3 for the spin gap and from approximately 3/2 to 4/3 for the energy per site and error bars in the exponent also decrease as J_2 approaches to J_2c.