A hybrid quantum-classical theory for predicting terahertz charge-transfer plasmons in metal nanoparticles on graphene

Metal nanoparticle (NP) complexes lying on a single-layer graphene surface are studied with a developed original hybrid quantum-classical theory using the Finite Element Method (FEM) that is computationally cheap. Our theory is based on the motivated assumption that the carrier charge density in the doped graphene does not vary significantly during the plasmon oscillations. Charge transfer plasmon (CTP) frequencies, eigenvectors, quality factors, energy loss in the NPs and in graphene, and the absorption power are aspects that are theoretically studied and numerically calculated. It is shown the CTP frequencies reside in the terahertz range and can be represented as a product of two factors: the Fermi level of graphene and the geometry of the NP complex. The energy losses in the NPs are predicted to be inversely dependent on the radius R of the nanoparticle, while the loss in graphene is proportional to R and the interparticle distance. The CTP quality factors are predicted to be in the range ∼10-100. The absorption power under CTP excitation is proportional to the scalar product of the CTP dipole moment and the external electromagnetic field. The developed theory makes it possible to simulate different properties of CTPs 3-4 orders of magnitude faster compared to the original FEM or the finite-difference time domain method, providing possibilities for predicting the plasmonic properties of very large systems for different applications.

Triplet state harvesting and search for forbidden transition intensity in the nitrogen molecule

Triplet excited states of the N2 molecule play an important role in electric discharges through air or liquid nitrogen accompanied by various afterglows. In the rarefied upper atmosphere, they produce aurora borealis and participate in other energy-transfer processes connected with atmospheric photochemistry and nightglow. In this work, we present spin-orbit coupling calculations of the intensity of various forbidden transitions, including the prediction of the electric dipole transition moment of the new1 3 Σ g - ←   A 3 Σ u + band, which is strongly prohibited by the (+|-) selection rule, the new spin-induced magneticB ' 3 Σ u - ←   A 3 Σ u + transition, magnetic and electric quadrupole transitions for the B3Πg← X 1 Σ g + Wilkinson band, and the Lyman-Birge-Hopfield a1Πg ← X1Σg transition. Also, two other far-UV singlet-singlet quadrupole transitions are calculated for the first time, namely, the Dressler-Lutz a"1Σg +-X1Σg + and the less studied z1Δg-X1Σg + weak transitions.

Formation and relaxation of K−2 and K−2V double-core-hole states in n-butane

Using a magnetic bottle multi-electron time-of-flight spectrometer in combination with synchrotron radiation, double-core-hole pre-edge and continuum states involving the K-shell of the carbon atoms in n-butane ( n-C4H10) have been identified, where the ejected core electron(s) and the emitted Auger electrons from the decay of such states have been detected in coincidence. An assignment of the main observed spectral features is based on the results of multi-configurational self-consistent field (MCSCF) calculations for the excitation energies and static exchange (STEX) calculations for energies and intensities. MCSCF results have been analyzed in terms of static and dynamic electron relaxation as well as electron correlation contributions to double-core-hole state ionization potentials. The analysis of applicability of the STEX method, which implements the one-particle picture toward the complete basis set limit, is motivated by the fact that it scales well toward large species. We find that combining the MCSCF and STEX techniques is a viable approach to analyze double-core-hole spectra.

X-ray absorption of molecular cations-a new challenge for electronic structure theory

In this paper we put forward some historical notes on the development of computational chemistry toward applications of x-ray spectroscopies. We highlight some of the important contributions by Enrico Clementi as method and program developer and as a supporter of this branch of computational research. We bring up a modern example based on the very recent experimental development of x-ray absorption of cationic molecules. As we show this spectroscopy poses new challenges for electronic structure theory and the electron correlation problem.

First-principles calculations of anharmonic and deuteration effects on the photophysical properties of polyacenes and porphyrinoids

A new method for calculating internal conversion rate constants (k[combining low line]IC), including anharmonic effects and using the Lagrangian multiplier technique, is proposed. The deuteration effect on k[combining low line]IC is investigated for naphthalene, anthracene, free-base porphyrin (H2P) and tetraphenylporphyrin (H2TPP). The results show that anharmonic effects are important when calculating k[combining low line]IC for transitions between electronic states that are energetically separated (ΔE) by more than 20 000-25 000 cm-1. Anharmonic effects are also important when ΔE < 20 000-25 000 cm-1 and when the accepting modes are X-H stretching vibrations with a frequency larger than 2000 cm-1. The calculations show that there is mixing between the S1 and S2 states of naphthalene induced by non-adiabatic interactions. The non-adiabatic interaction matrix element between the S1 and S2 states is 250 cm-1 and 50 cm-1 for the normal and fully deuterated naphthalene structure and this difference significantly affects the estimated fluorescence quantum yield. Besides aromatic hydrocarbons H2P and H2TPP, the k[combining low line]IC rate constant is also calculated for pyrometene (PM567) and tetraoxa[8]circulene (4B) with a detailed analysis of the effect of the vibrational anharmonicity.

Interacting Dirac materials

We investigate the extent to which the class of Dirac materials in two-dimensions provides general statements about the behavior of both fermionic and bosonic Dirac quasiparticles in the interacting regime. For both quasiparticle types, we find common features for the interaction induced renormalization of the conical Dirac spectrum. We perform the perturbative renormalization analysis and compute the self-energy for both quasiparticle types with different interactions and collate previous results from the literature whenever necessary. Guided by the systematic presentation of our results in table 1, we conclude that long-range interactions generically lead to an increase of the slope of the single-particle Dirac cone, whereas short-range interactions lead to a decrease. The quasiparticle statistics does not qualitatively impact the self-energy correction for long-range repulsion but does affect the behavior of short-range coupled systems, giving rise to different thermal power-law contributions. The possibility of a universal description of the Dirac materials based on these features is also mentioned.

X-Ray Absorption Spectrum of the N_{2}^{+} Molecular Ion

The x-ray absorption spectrum of N_{2}^{+} in the K-edge region has been measured by irradiation of ions stored in a cryogenic radio frequency ion trap with synchrotron radiation. We interpret the experimental results with the help of restricted active space multiconfiguration theory. Spectroscopic constants of the 1σ_{u}^{-1} ^{2}Σ_{u}^{+} state, and the two 1σ_{u}^{-1}3σ_{g}^{-1}1π_{g} ^{2}Π_{u} states are determined from the measurements. The charge of the ground state together with spin coupling involving several open shells give rise to double excitations and configuration mixing, and a complete breakdown of the orbital picture for higher lying core-excited states.

Engineering novel tunable optical high-Q nanoparticle array filters for a wide range of wavelengths

The interaction of non-monochromatic radiation with arrays comprising plasmonic and dielectric nanoparticles has been studied using the finite-difference time-domain electrodynamics method. It is shown that LiNbO₃, TiO₂, GaAs, Si, and Ge all-dielectric nanoparticle arrays can provide a complete selective reflection of an incident plane wave within a narrow spectral line of collective lattice resonance with a Q-factor of 10³ or larger at various spectral ranges, while plasmonic refractory TiN and chemically stable Au nanoparticle arrays provide high-Q resonances with moderate reflectivity. Arrays with fixed dimensional parameters make it possible to fine-tune the position of a selected resonant spectral line by tilting the array relative to the direction of the incident radiation. These effects provide grounds for engineering novel selective tunable optical high-Q filters in a wide range of wavelengths, from visible to middle-IR.

Charge-transfer plasmons with narrow conductive molecular bridges: A quantum-classical theory

We analyze a new type of plasmon system arising from small metal nanoparticles linked by narrow conductive molecular bridges. In contrast to the well-known charge-transfer plasmons, the bridge in these systems consists only of a narrow conductive molecule or polymer in which the electrons move in a ballistic mode, showing quantum effects. The plasmonic system is studied by an original hybrid quantum-classical model accounting for the quantum effects, with the main parameters obtained from first-principles density functional theory simulations. We have derived a general analytical expression for the modified frequency of the plasmons and have shown that its frequency lies in the near-infrared (IR) region and strongly depends on the conductivity of the molecule, on the nanoparticle-molecule interface, and on the size of the system. As illustrated, we explored the plasmons in a system consisting of two small gold nanoparticles linked by a conjugated polyacetylene molecule terminated by sulfur atoms. It is argued that applications of this novel type of plasmon may have wide ramifications in the areas of chemical sensing and IR deep tissue imaging.

Collective lattice resonances in arrays of dielectric nanoparticles: a matter of size

Collective lattice resonances (CLRs) in finite-sized $ 2D $2D arrays of dielectric nanospheres have been studied via the coupled dipole approximation. We show that even for sufficiently large arrays, up to $ 100 \times 100 $100×100 nanoparticles (NPs), electric or magnetic dipole CLRs may differ significantly from the ones calculated for infinite arrays with the same NP sizes and interparticle distances. The discrepancy is explained by the existence of a sufficiently strong cross-interaction between electric and magnetic dipoles induced at NPs in finite-sized lattices, which is ignored for infinite arrays. We support this claim numerically and propose an analytic model to estimate a spectral width of CLRs for finite-sized arrays. Given that most of the current theoretical and numerical researches on collective effects in arrays of dielectric NPs rely on modeling infinite structures, the reported findings may contribute to thoughtful and optimal design of inherently finite-sized photonic devices.

Aromaticity and photophysics of tetrasila- and tetragerma-annelated tetrathienylenes as new representatives of the hetero[8]circulene family

The electronic structure, absorption and emission spectra, aromaticity and photophysical behavior of the recently synthesized tetrasilatetrathia[8]circulene and tetragermatetrathia[8]circulene compounds have been studied computationally.

Identification of tautomeric intermediates of a novel thiazolylazonaphthol dye - A density functional theory study

The recently synthesized thiazolylazo dye, 1-[5-benzyl-1,3-thiazol-2-yl)diazenyl]naphthalene-2-ol called shortly BnTAN, is studied by density functional theory (DFT) in three tautomeric forms in order to explain the available ¹H NMR, UV-Vis and FTIR spectra. An experimentally observed IR band at 1678 cm^-1, assigned to the CO bond stretching vibration, supports the notion that BnTAN retains in the less stable keto-form even in the solid state due to an ultrafast single-coordinate intramolecular proton transfer. This finding is also in a good agreement with an X-ray crystallography analysis which indicates an intermediate position of the proton between the -OH and -N=N- groups. Calculations also show that some experimentally observed ¹H NMR signals could be considered as being averaged values between theoretically calculated chemical shifts for the corresponding protons in the keto- and enol-tautomers. At the same time the UV-Vis spectra are almost insensitive to the tautomerization processes as the main single band absorption at 500 nm is present in all tautomers according to our TD DFT simulations. The minor differences in spectral features of the long-wavelength visible region are also noted and discussed with respect to the manifestation of the less stable tautomer form.

Quantum-classical calculations of X-ray photoelectron spectra of polymers-Polymethyl methacrylate revisited

In this work, we apply quantum mechanics/molecular mechanics (QM/MM) approach to predict core-electron binding energies and chemical shifts of polymers, obtainable via X-ray photoelectron spectroscopy(XPS), using polymethyl methacrylate as a demonstration example. The results indicate that standard parametrizations of the quantum part (basis sets, level of correlation) and the molecular mechanics parts (decomposed charges, polarizabilities, and capping technique) are sufficient for the QM/MM model to be predictive for XPS of polymers. It is found that the polymer environment produces contributions to the XPS binding energies that are close to monotonous with the number of monomer units, totally amounting to approximately an eV decrease in binding energies. In most of the cases, the order of the shifts is maintained, and even the relative size of the differential shifts is largely preserved. The coupling of the internal core-hole relaxation to the polymer environment is found to be weak in each case, amounting only to one or two tenths of an eV. The main polymeric effect is actually well estimated already at the frozen orbital level of theory, which in turn implies a substantial computational simplification. These conclusions are best represented by the cases where the ionized monomer and its immediate surrounding are treated quantum mechanically. If the QM region includes only a single monomer, a couple of anomalies are spotted, which are referred to the QM/MM interface itself and to the neglect of a possible charge transfer.

Use of a self-rating scale to monitor depression severity in recurrent GP consultations in primary care - does it really make a difference? A randomised controlled study

Background

Little information is available about whether the use of self-assessment instruments in primary care affects depression course and outcome. The purpose was to evaluate whether using a depression self-rating scale in recurrent person-centred GP consultations affected depression severity, quality of life, medication use, and sick leave frequency.

Methods

Patients in the intervention group met their GP regularly at least 4 times during the 3 months intervention. In addition to treatment as usual (TAU), patients completed a self-assessment instrument (Montgomery-Asberg Depression Rating Scale) on each occasion, and then GPs used the completed instrument as the basis for a person-centred discussion of changes in depression symptoms. The control group received TAU. Frequency of visits in the TAU arm was the result of the GPs’ and patients’ joint assessments of care need in each case.

Depression severity was measured with Beck Depression Inventory-II (BDI-II), quality of life with EQ-5D, and psychological well-being with the General Health Questionnaire-12 (GHQ-12). Data on sick leave, antidepressant and sedatives use, and care contacts were collected from electronic patient records. All variables were measured at baseline and 3, 6, and 12 months. Mean intra-individual changes were compared between the intervention and TAU group.

Results

There were no significant differences between the intervention and control group in depression severity reduction or remission rate, change in quality of life, psychological well-being, sedative prescriptions, or sick leave during the whole 12-month follow-up. However, significantly more patients in the intervention group continued antidepressants until the 6 month follow-up (86/125 vs 78/133, p < 0.05).

Conclusions

When GPs used a depression self-rating scale in recurrent consultations, patients more often continued antidepressant medication according to guidelines, compared to TAU patients. However, reduction of depressive symptoms, remission rate, quality of life, psychological well-being, sedative use, sick leave, and health care use 4-12 months was not significantly different from the TAU group. These findings suggest that frequent use of depression rating scales in person-centred primary care consultations has no further additional effect on patients’ depression or well-being, sick leave, or health care use.

Trial registration

ClinicalTrials.gov Identifier: NCT01402206. Registered June 27 2011(retrospectively registered).

Electronic supplementary material

The online version of this article (doi:10.1186/s12875-016-0578-9) contains supplementary material, which is available to authorized users.

Early onset of inflammation during ontogeny of bipolar disorder: the NLRP2 inflammasome gene distinctly differentiates between patients and healthy controls in the transition between iPS cell and neural stem cell stages

Neuro-inflammation and neuronal communication are considered as mis-regulated processes in the aetiology and pathology of bipolar disorder (BD). Which and when specific signal pathways become abnormal during the ontogeny of bipolar disorder patients is unknown. To address this question, we applied induced pluripotent stem cell (iPSC) technology followed by cortical neural differentiation on adipocyte-derived cells from BD type I patients (with psychotic episodes in psychiatric history) and healthy volunteers (controls). RNA sequencing in iPSC and cortical neural stem cell (NSC) lines were used to examine alterations between the transcriptomes from BD I and control samples during transition from the pluripotent stage towards the neural developmental stage. At the iPSC stage, the most highly significant differentially expressed gene (DEG) was the NLRP2 inflammasome (P=2.66 × 10^-10). Also among 42 DEGs at the NSC stage, NLRP2 showed the strongest statistical significance (P=3.07 × 10^-19). In addition, we have also identified several cytoskeleton-associated genes as DEGs from the NSC stage, such as TMP2, TAGLN and ACTA2; the former two genes are recognised for the first time to be associated with BD. Our results also suggest that iPSC-derived BD-cortical NSCs carry several abnormalities in dopamine and GABA receptor canonical pathways, underlining that our in vitro BD model reflects pathology in the central nervous system. This would indicate that mis-regulated gene expression of inflammatory, neurotransmitter and cytoskeletal signalling occurs during early fetal brain development of BD I patients.

Suppression of surface plasmon resonance in Au nanoparticles upon transition to the liquid state

Significant suppression of resonant properties of single gold nanoparticles at the surface plasmon frequency during heating and subsequent transition to the liquid state has been demonstrated experimentally and explained for the first time. The results for plasmonic absorption of the nanoparticles have been analyzed by means of Mie theory using experimental values of the optical constants for the liquid and solid metal. The good qualitative agreement between calculated and experimental spectra support the idea that the process of melting is accompanied by an abrupt increase of the relaxation constants, which depends, beside electron-phonon coupling, on electron scattering at a rising number of lattice defects in a particle upon growth of its temperature, and subsequent melting as a major cause for the observed plasmonic suppression. It is emphasized that observed effect is fully reversible and may underlie nonlinear optical responses of nanocolloids and composite materials containing plasmonic nanoparticles and their aggregates in conditions of local heating and in general, manifest itself in a wide range of plasmonics phenomena associated with strong heating of nanoparticles.

Quantum mechanics capacitance molecular mechanics modeling of core-electron binding energies of methanol and methyl nitrite on Ag(111) surface

We study a newly devised quantum mechanics capacitance molecular mechanics (QMCMM) method for the calculation of core-electron binding energies in the case of molecules adsorbed on metal surfaces. This yet untested methodology is applied to systems with monolayer of methanol/methyl nitrite on an Ag(111) surface at 100 K temperature. It was found out that the studied C, N, and O 1s core-hole energies converge very slowly as a function of the radius of the metallic cluster, which was ascribed to build up of positive charge on the edge of the Ag slab. Further analysis revealed that an extrapolation process can be used to obtain binding energies that deviated less than 0.5 eV against experiments, except in the case of methanol O 1s where the difference was as large as 1.8 eV. Additional QM-cluster calculations suggest that the latter error can be connected to the lack of charge transfer over the QM-CMM boundary. Thus, the results indicate that the QMCMM and QM-cluster methods can complement each other in a holistic picture of molecule-adsorbate core-ionization studies, where all types of intermolecular interactions are considered.

An experimental and theoretical study of core-valence double ionisation of acetaldehyde (ethanal)

Core-valence double ionisation spectra of acetaldehyde (ethanal) are presented at photon energies above the carbon and oxygen 1s ionisation edges, measured by a versatile multi-electron coincidence spectroscopy technique. We use this molecule as a testbed for analyzing core-valence spectra by means of quantum chemical calculations of transition energies. These theoretical approaches range from two simple models, one based on orbital energies corrected by core valence interaction and one based on the equivalent core approximation, to a systematic series of quantum chemical electronic structure methods of increasing sophistication. The two simple models are found to provide a fast orbital interpretation of the spectra, in particular in the low energy parts, while the coverage of the full spectrum is best fulfilled by correlated models. CASPT2 is the most sophisticated model applied, but considering precision as well as computational costs, the single and double excitation configuration interaction model seems to provide the best option to analyze core-valence double hole spectra.

Calcium-dependent intracellular signal pathways in primary cultured adipocytes and ANK3 gene variation in patients with bipolar disorder and healthy controls

Bipolar disorder (BD) is a chronic psychiatric disorder of public health importance affecting >1% of the Swedish population. Despite progress, patients still suffer from chronic mood switches with potential severe consequences. Thus, early detection, diagnosis and initiation of correct treatment are critical. Cultured adipocytes from 35 patients with BD and 38 healthy controls were analysed using signal pathway reporter assays, that is, protein kinase C (PKC), protein kinase A (PKA), mitogen-activated protein kinases (extracellular signal-regulated kinase (ERK) and c-Jun N-terminal kinase (JNK)), Myc, Wnt and p53. The levels of activated target transcriptional factors were measured in adipocytes before and after stimulation with lithium and escitalopram. Variations were analysed in the loci of 25 different single-nucleotide polymorphisms (SNPs). Activation of intracellular signals in several pathways analysed were significantly higher in patients than in healthy controls upon drug stimulation, especially with escitalopram stimulation of PKC, JNK and Myc, as well as lithium-stimulated PKC, whereas no meaningful difference was observed before stimulation. Univariate analyses of contingency tables for 80 categorical SNP results versus diagnoses showed a significant link with the ANK3 gene (rs10761482; likelihood ratio χ(2)=4.63; P=0.031). In a multivariate ordinal logistic fit for diagnosis, a backward stepwise procedure selected ANK3 as the remaining significant predictor. Comparison of the escitalopram-stimulated PKC activity and the ANK3 genotype showed them to add their share of the diagnostic variance, with no interaction (15% of variance explained, P<0.002). The study is cross-sectional with no longitudinal follow-up. Cohorts are relatively small with no medication-free patients, and there are no 'ill patient' controls. It takes 3 to 4 weeks of culture to expand adipocytes that may change epigenetic profiles but remove the possibility of medication effects. Abnormalities in the reactivity of intracellular signal pathways to stimulation and the ANK3 genotype may be associated with pathogenesis of BD. Algorithms using biological patterns such as pathway reactivity together with structural genetic SNP data may provide opportunities for earlier detection and effective treatment of BD.

Studies of pH-Sensitive Optical Properties of the deGFP1 Green Fluorescent Protein Using a Unique Polarizable Force Field

The aim of this study is to identify the responsible molecular forms for the pH dependent optical properties of the deGFP1 green fluorescent protein mutant. We have carried out static and dynamic type calculations for all four protonation states of the chromophore to unravel the contributions due to finite temperature and the flexible protein backbone on the pH dependent optical properties. In particular, we have used a combined molecular dynamics and density functional-molecular mechanics linear response approach by means of which the optical property calculations were carried out for the chromophore in the explicitly treated solvent and bioenvironment. Two different models were used to describe the environment-electronic embedding and polarizable electronic embedding-accounting for the polarization of the chromophore and the mutual polarization between the chromophore and the environment, respectively. For this purpose a polarizable force field was derived quantum mechanically for the protein environment by use of analytical response theory. While the gas-phase calculations for the chromophore predict that the induced red shift going from low to high pH is attributed to the change of molecular forms from neutral to zwitterionic, the two more advanced models that explicitly account for the protein backbone attribute the pH shift to a neutral to anionic conversion. Some ramifications of the results for the use of GFPs as pH sensors are discussed.

Symmetry breaking and hole localization in multiple core electron ionization

Motivated by recent opportunitites to study hollow molecules with multiple core holes offered by X-ray free electron lasers, we revisit the core-hole localization and symmetry breaking problem, now studying ionization of more than one core electron. It is shown, using a N2 molecule with one, two, three, and four core holes, for example, that in a multiconfigurational determination of the core ionization potentials employing a molecular point group with broken inversion symmetry, one particular configuration is sufficient to account for the symmetry breaking relaxation energy in an independent particle approximation in the case of one or three holes, whereas the choice of point group symmetry is unessential for two and four holes. The relaxation energy follows a quadratic dependence on the number of holes in both representations.

Role of zero-point vibrational corrections to carbon hyperfine coupling constants in organic π radicals

By analyzing a set of organic π radicals, we demonstrate that zero-point vibrational corrections give significant contributions to carbon hyperfine coupling constants, in one case even inducing a sign reversal for the coupling constant. We discuss the implications of these findings for the computational analysis of electron paramagnetic spectra based on hyperfine coupling constants evaluated at the equilibrium geometry of radicals. In particular, we note that a dynamical description that involves the nuclear motion is in many cases necessary in order to achieve a semi-quantitatively predictive theory for carbon hyperfine coupling constants. In addition, we discuss the implications of the strong dependence of the carbon hyperfine coupling constants on the zero-point vibrational corrections for the selection of exchange-correlation functionals in density functional theory studies of these constants.

Dissociative x-ray lasing

X-ray lasing is predicted to ensue when molecules are pumped into dissociative core-excited states by a free-electron-laser pulse. The lasing is due to the population inversion created in the neutral dissociation product, and the process features self-trapping of the x-ray pulse at the gain ridge. Simulations performed for the HCl molecule pumped at the 2p(1/2)→6σ resonance demonstrate that the scheme can be used to create ultrashort coherent x-ray pulses.

Spin-spin and spin-orbit interactions in nanographene fragments: a quantum chemistry approach

The relativistic behavior of graphene structures, starting from the fundamental building blocks--the poly-aromatic hydrocarbons (PAHs) along with other PAH nanographenes--is studied to quantify any associated intrinsic magnetism in the triplet (T) state and subsequently in the ground singlet (S) state with account of possible S-T mixture induced by spin-orbit coupling (SOC). We employ a first principle quantum chemical-based approach and density functional theory (DFT) for a systematic treatment of the spin-Hamiltonian by considering both the spin-orbit and spin-spin interactions as dependent on different numbers of benzene rings. We assess these relativistic spin-coupling phenomena in terms of splitting parameters which cause magnetic anisotropy in absence of external perturbations. Possible routes for changes in the couplings in terms of doping and defects are also simulated and discussed. Accounting for the artificial character of the broken-symmetry solutions for strong spin polarization of the so-called "singlet open-shell" ground state in zigzag graphene nanoribbons predicted by spin-unrestricted DFT approaches, we interpolate results from more sophisticated methods for the S-T gaps and spin-orbit coupling (SOC) integrals and find that these spin interactions become weak as function of size and increasing decoupling of electrons at the edges. This leads to reduced electron spin-spin interaction and hence almost negligible intrinsic magnetism in the carbon-based PAHs and carbon nanographene fragments. Our results are in agreement with the fact that direct experimental evidence of edge magnetism in pristine graphene has been reported so far. We support the notion that magnetism in graphene only can be ascribed to structural defects or impurities.