Search

As the existing techniques to directly detect dark matter particles in laboratory experiments continue to produce negative signals, the research community considers novel approaches including those relying on condensed matter systems with non-trivial quantum properties. For example, finding signatures of dark matter in transport characteristics of solids would be an important step on the road to detect this illusive component of the mass of our Universe. As a first step in this direction we have recently derived the modified kinetic equation taking into account two coupled $U(1)$-gauge fields, one being the standard Maxwell electromagnetic field and other corresponding to the dark sector. The resulting Boltzmann kinetic equation is modified by the Berry curvature which couples to both visible and dark sector gauge fields. The linear and nonlinear longitudinal and/or Hall currents in topological matter may signal the existence of dark particle induced electromagnetic fields. The usage of appropriately designed resonant cavity to enhance the electric field may help to reach the goal. The conservative estimates show that the technique should have the sensitivity of the order of $10^{-9}$ or better and may help to partially fill the gap between XENON1T and haloscope experiments., Comment: 25 pages. 2 figures Revtex, to be published in Phys.Rev.D15

Finding signatures of dark matter in transport characteristics of solids would be an important step on the road to detect this illusive component of the mass of our Universe. This is especially important and timely as the experiments designed to directly detect dark matter particles continue to provide negative results. As a first step in this direction we consider topologically non-trivial Weyl or Dirac semimetals and derive the modified kinetic equation taking into account two coupled $U(1)$-gauge fields, one being the standard Maxwell electromagnetic field and other corresponding to the dark sector. The resulting Boltzmann kinetic equation is modified by the Berry curvature which couples to both visible and dark sector gauge fields. It was revealed that the dark sector induces modifications of transport coefficients due to the appearance of coupling constant between gauge fields and dark sector magnetic field., Comment: RevTex, 23 pages, to be published in Physics. Rev. D15

We propose an efficient method of heat rectification in a simple system consisting of a quantum dot asymmetrically coupled to four mutually perpendicular electrodes. In such a device the Hall-like charge and heat currents appear in response to the voltage bias or temperature difference between one pair of electrodes. Even though both longitudinal (along the bias) and Hall-like (perpendicular to the bias) currents are rectified under appropriate conditions, the rectification factor is typically much bigger for the latter currents. This is true for heat and charge flow. The Hall-like currents are predicted to exist in linear as well as non-linear transport regimes and require broken mirror symmetry but not time reversal symmetry. The linear effect exists only in geometry which breaks two inversion symmetries along two pairs of electrically coupled terminals. The proposed system is attainable within current technology and provides novel platform of simultaneous heat and charge management at the nanoscale., Comment: 10 pages, 9 figures

The transport properties of junctions composed of a central region tunnel-coupled to external electrodes are frequently studied within the single-impurity Anderson model with Hubbard on-site interaction. In the present work, we supplement the model with an important ingredient, namely the charge-bond interaction, also known as correlated or assisted hopping. Correlated hopping enters the second-quantised Hamiltonian, written in the Wannier representation, as an off-diagonal many-body term. Using the equation of motion technique, we study the effect of the correlated hopping on the spectral and transport characteristics of a two-terminal quantum dot. Two different Green functions (GFs) appear: one of them describes the spectral properties of the quantum dot, the other the transport properties of the system. The calculation of the transport GF requires the knowledge of the spectral one. We use decoupling procedures similar to those which properly describe the standard Anderson model within the Kondo regime and outside of it. For an arbitrary ratio $x$ between the amplitudes of correlated and single-particle hopping terms, the transport GF fulfils the $x \leftrightarrow 2-x$ symmetry of the model. The average occupation of the dot also obeys this symmetry, albeit the spectral function of the quantum dot, calculated within an analogous decoupling scheme as for the transport GF, does not. We identify the physical reason for this behavior, and propose a way to cure it. Since the correlated-hopping term breaks the particle-hole symmetry of the model and modifies all transport characteristics of the system, the detailed knowledge of its influence on measurable characteristics is a prerequisite for its experimental detection. Simple, experimentally feasible methods are proposed., Comment: 20 pages, 12 figures

The anomalous spin Hall conductivity in the holographic model of Dirac semimetals with two Dirac nodes protected by the crystal symmetry has been elaborated. Such system besides the chiral anomaly possesses another anomaly which is related to the $\mathbb{Z}_2$ topological charge of the system. The holographic model of the system contains matter action with two $U(1)$-gauge fields as well as the appropriate combination of the Chern-Simons gauge terms. We also allow for the coupling of two gauge fields {\it via} the kinetic mixing parametrised by the coupling $\alpha$. The holographic approach in the probe limit enables us to obtain Hall conductivity. The aim of this work is to describe the phase transitions in the $\mathbb{Z}_2$ Dirac semimetals between the topologically trivial and non-trivial phases. Interestingly the anomalous Hall conductivity plays a role of the order parameter of this phase transition. The holographically found prefactor of the Hall conductivity in the topologically non-trivial phase, depends on the coupling $\alpha$ and the Chern - Simons couplings., Comment: 28 pages, 11 figures, RevTex

Statistical properties of the electron transport flowing through nanostructures are strongly influenced by the interactions, geometry of the system and/or by type of the external electrodes. These factors affect not only the average current induced in the system but also contribute to fluctuations in the flux of charges and their correlations. Due to possible applications of the hybrid nano-systems, containing one or more superconducting electrodes, a detailed understanding of the flow of charge and its fluctuations seems to be of primary importance. Coulomb repulsion between electrons usually strongly affect the current-current correlation function. In this work we study the correlations in the charge flow through such interacting quantum dot contacted to one superconducting and two normal electrodes. This set-up allows for analysis of the Andreev scattering events in the correlations of currents flowing between external electrodes and, in particular, gives access to cross-correlations between currents from/to different normal electrodes. Our approach relies on the master equation technique, which properly captures the Coulomb interactions. We study the finite frequency correlations and find the relaxation processes, related to the high frequency charge and low frequency polarisation fluctuations. The multiterminal structure of here studied single electron device allows to analyze a competition between the intra- and inter-channel correlations. In the appropriate limit of the interacting quantum dot embedded between two normal electrodes our calculations quantitatively describe the recent experimental data on the frequency dependent correlations. This shows a promising potential of the method for description of the hybrid systems with superconducting electrode(s).

The carriers in graphene tuned close to the Dirac point envisage signatures of the strongly interacting fluid and are subject to hydrodynamic description. The important question is whether strong disorder induces the metal-insulator transition in this two-dimensional material. The bound on the conductivity tensor found earlier within the single current description, implies that the system does not feature metal-insulator transition. The linear spectrum of the graphene imposes the phase-space constraints and calls for the two-current description of interacting electron and hole liquids. Based on the gauge/gravity correspondence, using the linear response of the black brane with broken translation symmetry in Einstein-Maxwell gravity with the auxiliary U(1)-gauge field, responsible for the second current, we have calculated the lower bound of the DC-conductivity in holographic model of graphene. The calculations show that the bound on the conductivity depends on the coupling between both U(1) fields and for a physically justified range of parameters it departs only weakly from the value found for a model with the single U(1) field., Comment: 13 pages, Revtex, to be published in Phys.Rev.D15

The quest for good thermoelectric materials and/or high-efficiency thermoelectric devices is of primary importance from theoretical and practical points of view. Low-dimensional structures with quantum dots or molecules are promising candidates to achieve the goal. Interactions between electrons, far-from-equilibrium conditions and strongly non-linear transport are important factors affecting the usefulness of the devices. This paper analyses the thermoelectric power of a two-terminal quantum dot under large thermal $\Delta T$ and voltage $V$ biases as well as the performance of the three-terminal system as a heat engine. To properly characterise the non-linear effects under these conditions, two different Seebeck coefficients are introduced, generalizing the linear response expression. The direct calculations of thermally induced electric and heat currents show, in agreement with recent work, that the efficiency of the thermoelectric heat engine as measured by the delivered power is maximal far from equilibrium. Moreover, the strong Coulomb interactions between electrons on the quantum dot are found to diminish the efficiency at maximum power and the maximal value of the delivered power, both in the Kondo regime and outside of it., Comment: 29 pages, 8 figures

The gauge-gravity duality has been applied to examine the properties of holographic superconducting quantum device (SQUID), composed of two S-N-S Josephson junctions, influenced by dark sector modelled by the additional $U(1)$-gauge field coupled to the ordinary Maxwell one. The dark matter sector is known to affect the properties of superconductors and is expected to enter the current-phase relation. The kinetic mixing between two gauge fields provides a mechanism allowing for the conceivable observation of the effect. We find small but visible effect of the dark matter particle traversing the device, which shows up as a change of its maximal current., Comment: JCAP style, 22 pages, 9 figures

We calculate the magnetoconductivity of the Weyl semimetal with $\mathbb{Z}_2$ symmetry and chiral anomaly utilizing the recently developed hydrodynamic theory. The system in question will be influenced by magnetic fields connected with ordinary Maxwell and the second $U(1)$-gauge field, which is responsible for anomalous topological charge. The presence of chiral anomaly and $\mathbb{Z}_2$ anomalous charge endow the system with new transport coefficients. We start with the linear perturbations of the hydrodynamic equations and calculate the magnetoconductivity of this system. The holographic approach in the probe limit is implemented to obtain the explicit dependence of the longitudinal magnetoconductivities on the magnetic fields., Comment: 35 pages, 4 figures, LaTex, the title was changed, the version meets the printed one

We consider the hydrodynamical model of topological Dirac semi-metal possessing two Dirac nodes separated in momentum space along a rotation axis. It has been argued that the system in question, except the chiral anomaly, is endowed with the other one $\mathbb{Z}_2$. In order to model such a system we introduce two $U(1)$-gauge fields. The presence of the additional $\mathbb{Z}_2$ anomaly leads to the non-trivial modifications of hydrodynamical equations and to the appearance of new kinetic coefficients bounded with the vorticity and the magnetic parts of Maxwell and auxiliary $U(1)$-gauge fields., Comment: 17 pages, RevTex, meets the printed version

Based on the gauge/gravity correspondence we have calculated the thermoelectric kinetic and transport characteristics of the strongly interacting materials in the presence of perpendicular magnetic field. The 3+1 dimensional system with Dirac-like spectrum is considered as a strongly interacting one if it is close to the particle-hole symmetry point. Transport in such system has been modeled by the two interacting vector fields. In the holographic theory the momentum relaxation is caused by axion field and leads to finite values of the direct current transport coefficients. We have calculated conductivity tensor in the presence of mutually perpendicular electric and magnetic fields and temperature gradient. The geometry differs from that in which magnetic field lies in the same plane as an electric one and temperature gradient., Comment: 26 pages, 7 figures

The electrons in graphene for energies close to the Dirac point have been found to form strongly interacting fluid. Taking this fact into account we have extended previous work on the transport properties of graphene by taking into account possible interactions between the currents and adding the external magnetic field directed perpendicularly to the graphene sheet. The perpendicular mag- netic field B severely modifies the transport parameters. In the present approach the quantization of the spectrum and formation of Landau levels is ignored. Gauge/gravity duality has been used in the probe limit. The dependence on the charge density of the Seebeck coefficient and thermo-electric parameters {\alpha}ij nicely agree with recent experimental data for graphene. The holographic model allows for the interpretation of one of the fields representing the currents as resulting from the dark matter sector. For the studied geometry with electric field perpendicular to the thermal gradient the effect of dark sector has been found to modify the transport parameters but mostly in a quantitative way only. This makes difficult the detection of this elusive component of the Universe by studying transport properties of graphene., Comment: The title was changed, the actual version meets the published one

The shear viscosity to the entropy density ratio $\eta/s$ of the anisotropic superfluid has been calculated by means of the gauge/gravity duality in the presence of the {\it dark matter} sector. The {\it dark matter} has been described by the Yang-Mills field analogous to the one describing visible matter sector and it is assumed to interact with the visible field with coupling constant $\alpha$. Close to the superfluid transition temperature ($T_c$) the analytical solution has been given up to the leading order in a symmetry breaking parameter and the ratio of the gravitational constant and Yang-Mils coupling. The tensor element of ratio $\eta/s$ remains unaffected by the {\it dark matter} for the viscosity tensor in the plane perpendicular to the symmetry breaking direction (here $yz$). The temperature dependence and the linear correction in $(1-\alpha)$ in the plane containing this direction (here $xy$) was also revealed. The correction linearly vanishes for temperature tending to the critical one $T\rightarrow T_c$., Comment: 14 pages, some comments and references added

Holographic model of a three-dimensional current carrying superconductor or superfluid with {\it dark matter} sector described by the additional $U(1)$-gauge field coupled to the ordinary Maxwell one, has been studied in the probe limit. We investigated analytically by the Sturm-Liouville variational method, the holographic s-wave and p-wave models in the background of the AdS soliton as well as five-dimensional AdS black hole spacetimes. The two models of p-wave superfluids were considered, the so called $SU(2)$ and the Maxwell-vector. Special attention has been paid to the dependence of the critical chemical potential and critical transition temperature on the velocity of the condensate and {\it dark matter} parameters. The current $J$ in holographic three-dimensional superconductor studied here, shows the linear dependence on $T_c-T$ for both s and p-wave symmetry. This is in a significant contrast with the previously obtained results for two-dimensional superconductors, which reveal the $(T-T_c)^{3/2}$ temperature dependence. The coupling constant $\alpha$, as well as, chemical potential $\mu_D$ and the velocity $S_D$ of the {\it dark matter}, affect the critical chemical potential of the p-wave holographic $SU(2)$ system. On the other hand, $\alpha$, {\it dark matter} velocity $S_D$ and density $\rho_D$ determine the actual value of the transition temperature of the same superconductor/superfluid set up. However, the {\it dark matter} does not affect the value of the current., Comment: Completely rewritten version. Added new references and the discussion of experimental relevance. 36 pages. JHEP style - file included

Based on the gauge/gravity correspondence, the hydrodynamic response coefficients, shear and Hall viscosities, have been studied. The holographic model of Einstein-Maxwell- AdS $(3+1)$-dimensional system additionally coupled with the another gauge field mimicking the {\it dark matter} sector, as well as, gravitational Chern-Simons term bounded with a dynamical scalar field, were taken into account. Condensation of the scalar field provides the parity violating term. Both shear and Hall viscosities have been calculated and their dependence on $\alpha$ - the coupling constant between matter and {\it dark matter} sectors has been studied. To the lowest order in the derivative expansion and perturbation in $\alpha$, the shear viscosity is not influenced by the {\it dark matter}, while the Hall component linearly depends on $\alpha$., Comment: JHEP style, this version meets the published one, JHEP08 (2016) 124

The holographic approach to study strongly coupled superconductors in the presence of dark matter is reviewed. We discuss the influence of dark matter on the superconducting transition temperature of both s-wave and p-wave holographic superconductors. The upper critical field, coherence length, penetration depth of holographic superconductors as well as the metal-insulator transitions have also been analysed. Issues related to the validity of AdS/CFT correspondence for the description of superconductors studied in the laboratory and possible experiments directed towards the detection of dark matter are discussed. In doing so we shall compare our assumptions and assertions with those generally accepted in the elementary particle experiments aimed at the detection of dark matter particles., Comment: 5+ pages, 1 figure, National Conference on Superconductivity 2015, Karpacz, Poland

The thermoelectric effects in three-terminal structures with a quantum dot are considered. We propose the experimentally consistent protocol for determination of the transport coefficients in terms of the local and non-local conductances and thermopowers that can be measured in two steps, applying the `four probe technique'. This proposal is compared with other approaches discussed so far in the literature. As an example we study in detail the thermopower induced by the superconducting electrode in a subgap regime which might be useful for analysis of novel hybrid devices.

Charge and heat transport in nano-structures is dominated by non-equilibrium effects which strongly influence their behaviour. These effects are studied in a setup consisting of three external leads, one of which is considered as a heat reservoir and is tunnel-coupled to two cold electrodes via two independently controlled quantum dots. The energy flow from the hot electrode together with energy filtering provided by quantum dots leads to a voltage bias between the cold electrodes. The heat and charge currents in the device effectively flow in mutually perpendicular directions, allowing for their independent control. The non-equilibrium screening changes the values of the system parameters needed for its optimal performance but leaves the maximal output power and efficiency unchanged. Our results are important from the theoretical point of view as well as for the practical implementation and the control of the proposed heat engine., Comment: 15 pages, 6 figures

The {\it dark matter} seem to be an inevitable ingredient of the total matter configuration in the Universe and the knowledge how the {\it dark matter} affects the properties of superconductors is of vital importance for the experiments aimed at its direct detection. The homogeneous magnetic field acting perpendicularly to the surface of (2+1) dimensional s-wave holographic superconductor in the theory with {\it dark matter} sector has been modeled by the additional $U(1)$-gauge field representing dark matter and coupled to the Maxwell one. As expected the free energy for the vortex configuration turns out to be negative. Importantly its value is lower in the presence of {\it dark matter} sector. This feature can explain why in the Early Universe first the web of {\it dark matter} appeared and next on these gratings the ordinary matter forming cluster of galaxies has formed., Comment: 23 pages, JHEP-style

We analytically study the phase transitions between s-wave holographic insulator/superconductor and metal/superconductor. The problem is solved by the variational method for the Sturm-Liouville eigenvalue problem in the theory with dark matter sector of U(1)-gauge field coupled to the Maxwell field. Additionally in the probe limit we investigate the marginally stable modes of scalar perturbations in the AdS solitonic background, connected with magnetic field in the dark matter sector. We have found that even with dark matter sector the superconducting transition temperature $T_c$ is proportional to charge density $\rho$ in power 1/3. This value seem to be strong coupling modification of the exponent 2/3 known from the Bose - Einstein condensation of charged local pair bosons in narrow band superconductors. The holographic droplet solution is affected by the coupling to the dark matter. Interestingly in the probe limit the critical chemical potential increases with the decreasing coupling to dark matter making the condensation transition harder to appear., Comment: 16 pages, REvTEx, 3 Figures, to be published in Phys.Rev.D15

The holographic approach to building the p-wave superconductors results in three different models: the Maxwell-vector, the SU(2) Yang-Mills and the helical. In the probe limit approximation, we analytically examine the properties of the first two models in the theory with {\it dark matter} sector. It turns out that the effect of {\it dark matter} on the Maxwell-vector p-wave model is the same as on the s-wave superconductor studied earlier. For the non-Abelian model we study the phase transitions between p-wave holographic insulator/superconductor and metal/superconductor. Studies of marginally stable modes in the theory under consideration allow us to determine features of p-wave holographic droplet in a constant magnetic field. The dependence of the superconducting transition temperature on the coupling constant $\alpha$ to the {\it dark matter} sector is affected by the {\it dark matter} density $\rho_D$. For $\rho_D>\rho$ the transition temperature is a decreasing function of $\alpha$. The critical chemical potential $\mu_c$ for the quantum phase transition between insulator and metal depends on the chemical potential of dark matter $\mu_D$ and for $\mu_D=0$ is a decreasing function of $\alpha$., Comment: 25 pages, 3 figures, JHEP style (included), version accepted for publication in JHEP

The interplay between various many body effects in a quantum dot attached to two normal and one superconducting lead is considered in the limit of large superconducting gap. By the proximity effect the superconducting lead induces pairing correlations on the quantum dot. In the subgap region one observes the anomalous tunneling via direct and crossed Andreev scattering, whereas the usual single particle electronic transfer is suppressed. The interactions of electrons on the dot leading to such phenomena as the Coulomb blockade and the Kondo effect severely modify the currents flowing in the system. In particular: (i) they prevent the existence of the negative differential conductance observed for non-interacting quantum dot over the whole range of voltages, (ii) affect the distribution of the currents as function of the applied voltage and (iii) lead to the appearance of additional low bias feature due to the formation of the Abrikosov-Suhl resonance. The non-local correlations in the Coulomb blockade regime are most pronounced for the particle-hole symmetric dot and thus can be easily tuned by means of gate voltage. They are observed even in the Kondo regime and dominate the behavior close to the Abrikosov-Suhl resonance showing convincingly that Kondo correlations do not destroy subtle entanglement between electrons., Comment: 14 pages, 10 figures, to appear in Phys. Rev. B

The interpretation of Kerr rotation measurements in the superconducting phase of Sr$_2$RuO$_4$ is a controversial topic. Both intrinsic and extrinsic mechanisms have been proposed, and it has been argued that the intrinsic response vanishes by symmetry. We focus on the intrinsic contribution and clarify several conflicting results in the literature. On the basis of symmetry considerations and detailed calculations we show that the intrinsic Kerr signal is not forbidden in a general multi- band system but has a rich structure in the near infrared regime. We distinguish different optical transitions determined by the superconducting gap (far infrared) and the inter orbital coupling of the normal state (near infrared). We argue that the low frequency transitions do not contribute to the Hall conductivity while only the inter-orbital transitions in the near infrared regime contribute. Finally, we discuss the difficulties to connect the calculations for the optical Hall conductivity to the experimental measurement of the Kerr angle. We will compare different approximations which might lead to conflicting results., Comment: 9 pages, 8 figures, 1 table

On the occasion of centenary of superconductivity discovery I recall some facts from the first period and attempts to understand the phenomenon. It turns out that most famous physicists of the first half of XX century have tried to solve the puzzle. Bardeen, Cooper and Schrieffer succeeded in 1957. The BCS theory successfully described all known facts and offered new predictions, which soon have been confirmed experimentally contributing to the widespread acceptance of the theory. It have found applications in nuclear physics, theory of neutron stars and cold atomic gases. The discoveries of new superconductors in the last thirty years show that simple BCS model is not enough to understand new unconventional superconductors. The studies of superconductors develop vividly and still fascinate new generations of physicists working in such diverse fields as material science and string theory., Comment: 4+ pages, 3 figures

We present an analysis of the Hall conductivity $\sigma_{xy}(\omega, T)$ in time reversal symmetry breaking states of exotic superconductors. We find that the dichroic signal is non-zero in systems with inter-band order parameters. This new intrinsic mechanism may explain the Kerr effect observed in strontium ruthenate and possibly other superconductors. We predict coherence factor effects in the temperature dependence of the imaginary part of the ac Hall conductivity $ Im\sigma_{xy}(\omega, T)$, which can be tested experimentally., Comment: 4+ pages, 4 figures, published version

The calculations of the local and global properties of two band superconductors have been presented with particular attention to the role of the inter-orbital scattering of pairs. The properties of such superconductors are very different from a single band or typical two band systems with dominant intra-band pairing interactions. The role of Van Hove singularity in one of the bands on the properties of intra-band clean superconductor has been discussed. It leads to marked increase of superconducting transition temperature in the weak coupling limit. We study the inhomogeneous systems in which the characteristics change from place to place by solving the Bogolubov-de Gennes equations for small clusters. The suppression of the superconducting order parameter by the single impurity scattering the fermions between bands is contrasted with that due to intra-band impurity scattering. The results obtained for impure systems have been shown as a maps of local density of states, the order parameter and gap function. They can be directly compared with STM spectra of the real material., Comment: 10 pages 14 figures embedded, revtex format

The boson-fermion (BF) model has been proposed to describe superconductivity in short coherence length superconductors. In this work we use it to study impurity induced inhomogeneities in high temperature superconductors (HTS) as found in numerous scanning tunnelling measurements. The model has been formulated in real space and solved with help of Bogolubov-de Gennes approach. Disorder in the boson or fermion subsystem is directly coupled to the superconducting order parameter and leads to severe changes of superconducting properties like local order parameter and density of states. We present the results for many impurities randomly distributed over otherwise clean and periodic two-dimensional square lattice., Comment: 5 pages with 5 Postscript figures, presented at 13 Krajowa Szkola Nadprzewodnictwa (Ladek-Zdroj 2007), will be published in APP-A, uses appolb.cls

The effect of magnetic field on specific heat of a p-wave spin triplet superconductor has been analyzed. To describe superconductivity we used a single band model with various realizations of p-wave order parameter and analyzed the effect of a Doppler shift in different temperatures., Comment: 5 pages, 4 figures

The calculation of the charge present in central region of the double barrier structure at non-equilibrium conditions is discussed. We propose here a simple method to calculate non equilibrium Green's functions which allows consistent calculations of retarded and distribution functions. To illustrate the approach we calculate the charge on the quantum dot coupled {\it via} tunnel barriers to two external leads having different chemical potentials $\mu_L$ and $\mu_R$. The obtained results have been compared with other approaches existing in the literature. They all agree in the equilibrium situation and the departures grow with increasing the difference $\mu_L-\mu_R$., Comment: 9 pages, 2 (.eps) figures, to be published in Solid State Commun

I address here the question of the mutual interplay of strong correlations and disorder in the system. I consider random version of the Hubbard model. Diagonal randomness is introduced {\it via} random on-site energies and treated by the coherent potential approximation. Strong, short ranged, electron - electron interactions are described by the slave boson technique and found to induce additional disorder in the system. As an example I calculate the density of states of the random interacting binary alloy and compare it with that for non-interacting system., Comment: 5 pages, 2 figures, requires jltp.sty (included)