Second-harmonic generation in GFP-like proteins

Theoretically designed M@diaza[2.2.2]cryptand complexes: the role of non-covalent interactions in promoting NLO properties of organic electrides

Organic excess electron compounds with significant nonlinear optical (NLO) properties are widely employed in optoelectronic applications. Herein, single-alkali metals with diaza[2.2.2] cryptand (M@crypt,M=Li, Na, and K) are investigated for optoelectronic and NLO properties by using the density functional theory. Thermodynamic and kinetic stabilities of present complexes are computed through interaction energy (Eint) and ab-initio molecular dynamic (AIMD) simulations. M@crypt complexes carry excess electrons and mimic molecular electrides. Quantum theory of atoms in molecules (QTAIM) analysis and reduced density gradient (RDG) spectra demonstrate the roles of the weak van der Waals (vdW) interactions between metal and complexant. The remarkable hyperpolarizability (βo) value up to 1.41 × 106 au may be credited to the presence of loosely bound excess electrons. The hyper Rayleigh scattering hyperpolarizability (βHRS) is recorded up to 1.31 × 106 au for the K@crypt. Furthermore, frequency-dependent first-order and second-order hyperpolarizability is more prominent at the applied frequency of ω = 0.042823 au. The electron localizing function (ELF) and localized orbital locator (LOL) analysis further disclose the nature of interaction between alkali metal and complexant. The TD-DFT method is adopted to get excited state parameters and absorbance properties. An electron density difference map (EDDM) is exploited to evaluate the orbital contributions in excited states. Hence, the studied electride may become a promising candidate for NLO materials. We anticipate that the present work will provide insight into further development of molecular electride for optoelectronic applications.

TDDFT Investigation of the Raman and Resonant Raman Spectra of Fluorescent Protein Chromophore Models

The off-resonance and resonant Raman spectra have been simulated for models of fluorescent protein chromophores, those of the green fluorescent protein (GFP, called FP1) and of DsRed (called FP2), which presents a longer π-conjugated path, with the aim of providing a systematic investigation of structural but also computational aspects. These were performed at the (time-dependent) density functional theory [(TD)DFT] level. The off-resonance intensities have been calculated from the derivatives of the frequency-dependent polarizability with respect to the normal coordinates while the resonant ones have been evaluated using Huang-Rhys factors determined from the gradients of the excitation energies with respect to the normal coordinates. When applied with the M05 meta-GGA exchange-correlation functional, this simple computational scheme can reproduce most of the experimental Raman signatures of FP1 in its protonated and deprotonated forms, the differences of vibrational signatures of the cis (Z) and trans (E) isomers, as well as their changes as a function of the excitation wavelength. On the other hand, testing the predictions made for FP2 would require new experimental work. It was also observed that simulations with methods that inadequately predict the resonant Raman spectra could nevertheless produce a UV-vis absorption spectrum that is quite similar to the one obtained with better methods, once realistic peak broadening has been applied.

Second-Order Nonlinear Optics Response of the Boron-Dipyrromethenes-Based Mislinked Expanded Porphyrins: Revealing the Role of the -BF2 Group

Here, the mislinked expanded porphyrins singly (labeled A) and doubly (labeled B) neo-confused [22]smaragdyrin, the boron-dipyrromethenes-based mislinked expanded porphyrins singly (labeled C) and doubly (labeled D) neo-confused [22]smaragdyrin, where both C and D include a -BF₂ group, are chosen to serve as the study objects, and theoretical calculations are carried out to study the role of the -BF₂ group in the second-order nonlinear optics (NLO) behaviors. Results highlighted that the -BF₂ group plays an important role for the second-order behaviors in mislinked expanded porphyrins; namely, embedding the -BF₂ group well enhanced the hyper-Rayleigh scattering (HRS) value {β_HRS(0;0,0)}, C{βHRS(0;0,0)}A{βHRS(0;0,0)} = 2.0 and D{βHRS(0;0,0)}B{βHRS(0;0,0)} = 2.9, main owning to the fact that installing -BF₂ increases the electron delocalization degree and decreases the excited energy of the crucial excited state.

All-Atom Quantum Mechanical Calculation of the Second-Harmonic Generation of Fluorescent Proteins

Fluorescent proteins (FPs) are biotags of choice for second-harmonic imaging microscopy (SHIM). Because of their large size, computing their second-harmonic generation (SHG) response represents a great challenge for quantum chemistry. In this contribution, we propose a new all-atom quantum mechanics methodology to compute SHG of large systems. This is now possible because of two recent implementations: the tight-binding GFN2-xTB method to optimize geometries and a related version of the simplified time-dependent density functional theory (sTD-DFT-xTB) to evaluate quadratic response functions. In addition, a new dual-threshold configuration selection scheme is introduced to reduce the computational costs while retaining overall similar accuracy. This methodology was tested to evaluate the SHG of the proteins iLOV and bacteriorhodopsin (bR). In the case of bR, quantitative agreement with respect to experiment was reached for the out-of-resonance low-energy part of the β_HRS frequency dispersion. This work paves the way toward an accurate prediction of the SHG of large structures-a requirement for the design of new and improved SHIM biotags.

Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers

Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.

Harmonic Generation Microscopy 2.0: New Tricks Empowering Intravital Imaging for Neuroscience

Optical harmonic generation, e.g., second- (SHG) and third-harmonic generation (THG), provides intrinsic contrasts for three-dimensional intravital microscopy. Contrary to two-photon excited fluorescence (TPEF), however, they have found relatively specialized applications, such as imaging collagenous and non-specific tissues, respectively. Here we review recent advances that broaden the capacity of SHG and THG for imaging the central nervous system in particular. The fundamental contrast mechanisms are reviewed as they encode novel information including molecular origin, spectroscopy, functional probes, and image analysis, which lay foundations for promising future applications in neuroscience.

Three-stage pH-switchable organic chromophores with large nonlinear optical responses and switching contrasts

First three-stage pH-switchable second-order nonlinear optical chromophores are synthesized and characterized by tunable-wavelength (non)linear spectroscopy, showing remarkably different on–off sequences.

Long- and Short-Range Electrostatic Fields in GFP Mutants: Implications for Spectral Tuning

The majority of protein functions are governed by their internal local electrostatics. Quantitative information about these interactions can shed light on how proteins work and allow for improving/altering their performance. Green fluorescent protein (GFP) and its mutation variants provide unique optical windows for interrogation of internal electric fields, thanks to the intrinsic fluorophore group formed inside them. Here we use an all-optical method, based on the independent measurements of transition frequency and one- and two-photon absorption cross sections in a number of GFP mutants to evaluate these internal electric fields. Two physical models based on the quadratic Stark effect, either with or without taking into account structural (bond-length) changes of the chromophore in varying field, allow us to separately evaluate the long-range and the total effective (short- and long-range) fields. Both types of the field quantitatively agree with the results of independent molecular dynamic simulations, justifying our method of measurement.

Organized chromophoric assemblies for nonlinear optical materials: towards (sub)wavelength scale architectures

Photonic circuits are expected to greatly contribute to the next generation of integrated chips, as electronic integrated circuits become confronted with bottlenecks such as heat generation and bandwidth limitations. One of the main challenges for the state-of-the-art photonic circuits lies in the development of optical materials with high nonlinear optical (NLO) susceptibilities, in particular in the wavelength and subwavelength dimensions which are compatible with on-chip technologies. In this review, the varied approaches to micro-/nanosized NLO materials based on building blocks of bio- and biomimetic molecules, as well as synthetic D-π-A chromophores, have been categorized as supramolecular self-assemblies, molecular scaffolds, and external force directed assemblies. Such molecular and supramolecular NLO materials have intrinsic advantages, such as structural diversities, high NLO susceptibilities, and clear structure-property relationships. These "bottom-up" fabrication approaches are proposed to be combined with the "top-down" techniques such as lithography, etc., to generate multifunctionality by coupling light and matter on the (sub)wavelength scale.

Engineering a genetically-encoded SHG chromophore by electrostatic targeting to the membrane

Although second harmonic generation (SHG) microscopy provides unique imaging advantages for voltage imaging and other biological applications, genetically-encoded SHG chromophores remain relatively unexplored. SHG only arises from non-centrosymmetric media, so an anisotropic arrangement of chromophores is essential to provide strong SHG signals. Here, inspired by the mechanism by which K-Ras4B associates with plasma membranes, we sought to achieve asymmetric arrangements of chromophores at the membrane-cytoplasm interface using the fluorescent protein mVenus. After adding a farnesylation motif to the C-terminus of mVenus, nine amino acids composing its β-barrel surface were replaced by lysine, forming an electrostatic patch. This protein (mVe9Knus-CVIM) was efficiently targeted to the plasma membrane in a geometrically defined manner and exhibited SHG in HEK293 cells. In agreement with its design, mVe9Knus-CVIM hyperpolarizability was oriented at a small angle (~7.3°) from the membrane normal. Genetically-encoded SHG chromophores could serve as a molecular platform for imaging membrane potential.

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.

Impact of proton transfer phenomena on the electronic structure of model Schiff bases: an AIM/NBO/ELF study

Understanding of the electronic structure evolution due to a proton dynamics is a key issue in biochemistry and material science. This paper reports on density functional theory calculations of Schiff bases containing short, strong intramolecular hydrogen bonds where the bridged proton is located: (i) at the donor site, (ii) strongly delocalized, and (iii) at the acceptor site. The mobility of the bridged proton and its influence on the molecular structure and properties of the chosen Schiff base derivatives have been investigated on the basis of Atoms in Molecules, Natural Bond Orbitals, and Electron Localization Function theories. It has been observed that the extent of the bridged proton delocalization is strongly modified by the steric and inductive effects present in the studied compounds introduced by various substituents. It has been shown that: (i) potential energy profiles for the proton motion are extremely dependent on the substitution of the aromatic ring, (ii) the topology of the free electron pairs present at the donor∕acceptor site, as well as their electron populations, are affected qualitatively by the bridged proton position, (iii) the distortion of the molecular structure due to the bridged proton dynamics includes the atomic charge fluctuations, which are in some cases non-monotonic, and (iv) topology of the ELF recognizes events of proton detachment from the donor and attachment to the acceptor. The quantitative and qualitative results shed light onto molecular consequences of the proton transfer phenomena.

Marine natural products from the deep Pacific as potential non-linear optical chromophores

Theoretical analysis using quadratic response theory within the time-dependent density functional theory (TDDFT) formalism shows that the dermacozines, a group of phenazine-based compounds isolated from cultures of Dermacoccus abyssi found in the Mariana Trench, possess large first hyperpolarisability (β) values at common incident laser wavelengths that are highly sensitive to the degree and type of substitution of the core structure. The phenazine moiety is a versatile and tunable chromophore for non-linear optics and this work serves to highlight the potential that (marine) natural products, even those found in the darkest places on the planet, may have for aiding developments in optical materials design.

SCC-DFTB calculation of the static first hyperpolarizability: from gas phase molecules to functionalized surfaces

The performance of the self-consistent charge density functional tight binding (SCC-DFTB) method for calculating the first hyperpolarizability of π-conjugated compounds has been assessed with respect to results obtained with high-level ab initio methods and density functional theory (DFT). The SCC-DFTB method performs similarly or better than DFT with the PBE XC functional. Thus, if for small π-conjugated linkers SCC-DFTB can reproduce trends, for longer chains the first hyperpolarizabilities are overestimated. In the case of push-pull thiophenes, the β values are strongly overestimated, as it is also the case with the B3LYP and PBE XC functionals. On the other hand, the SCC-DFTB method closely reproduces the evolution of β in p-disubstituted benzenes as a function of the donor and acceptor groups, as estimated at the MP2 level. The reliability of SCC-DFTB to determine the bond length alternation and the dihedral angles between the aromatic rings has also been tackled, demonstrating that both are underestimated. Overall, the SCC-DFTB calculations are of the same quality as those performed with the conventional PBE XC functional on which the method was parameterized but the SCC-DFTB calculations are computationally very little demanding, and it can therefore be adopted for very large systems for screening nonlinear optical materials as well as for assessing structure-property relationships. This is illustrated with an application on the first hyperpolarizability of an indolino-oxazolidine molecular switch grafted on a SiO2 surface. This has enabled to pinpoint (i) the effect of the surface on the donor/acceptor character of the linking substituent, (ii) the impact of molecular orientation, (iii) the role of a spacer between the π-conjugated switch and the surface, (iv) the global effect of the surface on the β contrast, and also (v) the fact that the molecular switches can maintain this contrast when adsorbed.

Origin of the Surface-Induced First Hyperpolarizability in the C60/SiO2 System: SCC-DFTB Insight

Using the self-consistent charge density functional tight binding (SCC-DFTB) method, C60 molecules physisorbed on an α-quartz slab are shown to display a first hyperpolarizability, whereas, owing to their symmetry, both the α-quartz slab and C60 molecule have no first hyperpolarizabilities. A larger first hyperpolarizability is achieved when the lowest-lying (five- or six-membered) ring is situated in between two hydroxyl rows, rather than on top, because this situation favors orbital overlaps and charge transfer. Further analysis has demonstrated that (i) the first hyperpolarizability originates from the MO overlap and field-induced charge transfers from the neighboring substrate/adsorbate moieties but not to geometric relaxation of the C60 molecules at the interface and that (ii) larger first hyperpolarizabilities are associated with low surface coverage and with small distances between C60 and the surface. This contribution is a clear illustration of the emergence of second-order nonlinear optical responses (first hyperpolarizability) as a result of breaking the centrosymmetry.

SECOND-HARMONIC GENERATION OF FLUORESCENT PROTEINS

We theoretically study the second-order nonlinear optical properties of six fluorescent proteins (FPs), such as green fluorescent protein (GFP), BFP, enhanced BFP (eBFP), CFP, YFP, and DsRed. To begin with, the geometries of all these FP chromophores are optimized at B3LYP/6-311++G** level in a water medium and the polarized continuum model (PCM in water) method is adopted. Using a time-dependent density functional theory (TDDFT) method, the electronic structures and excited-state properties of chromophores are determined. Here we employ TDDFT combining with the sum-over-states (SOS) method to calculate the first-order hyperpolarizability for second-harmonic generation (SHG) optical process. Moreover, we discuss the origin of the nonlinear optical response and determine what caused the variation of first-order hyperpolarizability. Our calculations show that the charge transfers of π → π* in the central conjugated structure and p → π* charge transfers from the side chain R1 to conjugated structure of chromophores markedly affect the first-order hyperpolarizability.

Design and characterization of molecular nonlinear optical switches

Nanoscale structures, including molecules, supramolecules, polymers, functionalized surfaces, and crystalline/amorphous solids, can commute between two or more forms, displaying contrasts in their nonlinear optical (NLO) properties. Because of this property, they have high potential for applications in data storage, signal processing, and sensing. As potential candidates for integration into responsive materials, scientists have been intensely studying organic and organometallic molecules with switchable first hyperpolarizability over the past two decades. As a result of this, researchers have been able to synthesize and characterize several families of molecular NLO switches that differ by the stimulus used to trigger the commutation. These stimuli can include light irradiation, pH variation, redox reaction, and ion recognition, among others. The design of multistate (including several switchable units) and multifunctional (triggered with different stimuli) systems has also motivated a large amount of work, aiming at the improvement of the storage capacity of optical memories or the diversification of the addressability of the devices. In complement to the synthesis of the compounds and the characterization of their NLO responses by means of hyper-Rayleigh scattering, quantum chemical calculations play a key role in the design of molecular switches with high first hyperpolarizability contrasts. Through the latter, we can gain a fundamental understanding of the various factors governing the efficiency of the switches. These are not easily accessible experimentally, and include donor/acceptor contributions, frequency dispersion, and solvent effects. In this Account, we illustrate the similarities of the experimental and theoretical tools to design and characterize highly efficient NLO switches but also the difficulties in comparing them. After providing a critical overview of the different theoretical approaches used for evaluating the first hyperpolarizabilities, we report two case studies in which theoretical simulations have provided guidelines to design NLO switches with improved efficiencies. The first example presents the joint theoretical/experimental characterization of a new family of multi-addressable NLO switches based on benzazolo-oxazolidine derivatives. The second focuses on the photoinduced commutation in merocyanine-spiropyran systems, where the significant NLO contrast could be exploited for metal cation identification in a new generation of multiusage sensing devices. Finally, we illustrate the impact of environment on the NLO switching properties, with examples based on the keto-enol equilibrium in anil derivatives. Through these representative examples, we demonstrate that the rational design of molecular NLO switches, which combines experimental and theoretical approaches, has reached maturity. Future challenges consist in extending the investigated objects to supramolecular architectures involving several NLO-responsive units, in order to exploit their cooperative effects for enhancing the NLO responses and contrasts.

A THEORETICAL STUDY ON SECOND HARMONIC GENERATION HYPERPOLARIZABILITIES OF PHENYLALANINE POLYPEPTIDES

The second harmonic generation (SHG) hyperpolarizabilities of phenylalanine and homopolypeptides are investigated by configuration interaction among singly excited configurations (CIS) technique combined with the sum-over-states (SOS) method. The geometries of peptides containing phenylalanine ( Phe )n(n = 1–8) are optimized by B3LYP/6-31g(d) method, and they form the special structures like β-sheet (a common protein secondary structure). It is found that the energy gaps of various peptides are reduced and the hyperpolarizabilities are increased with the peptide chains lengthened. We discuss the origin of the second-order nonlinear optical response in phenylalanine homopolypeptides and confirm that the π → π* transitions in the aromatic residue of phenylalanine make the most important contributions to the second-order polarizability. Our results strongly suggest that the hyperpolarizabilities are dominated from the propagation direction of peptide chains.

Improving the second-order nonlinear optical response of fluorescent proteins: the symmetry argument

We have successfully designed and expressed a new fluorescent protein with improved second-order nonlinear optical properties. It is the first time that a fluorescent protein has been rationally altered for this particular characteristic. On the basis of the specific noncentrosymmetry requirements for second-order nonlinear optical effects, we had hypothesized that the surprisingly low first hyperpolarizability (β) of the enhanced yellow fluorescent protein (eYFP) could be explained by centrosymmetric stacking of the chromophoric Tyr66 and the neighboring Tyr203 residue. The inversion center was removed by mutating Tyr203 into Phe203, with minor changes in the linear optical properties and even an improved fluorescence quantum yield. Structure determination by X-ray crystallography as well as linear optical characterization corroborate a correct folding and maturation. Measurement of β by means of hyper-Rayleigh scattering (HRS) as well as their analysis using quantum chemistry calculations validate our hypothesis. This observation can eventually lead to improved red fluorescent proteins for even better performance. On the basis of the specific function (second-harmonic generation), the color of its emission, and in analogy with the "fruit" names, we propose SHardonnay as the name for this Tyr203Phe mutant of eYFP.

Looking at the Green Fluorescent Protein (GFP) chromophore from a different perspective: a computational insight

In the present contribution Density Functional Theory (DFT) has been applied to explore molecular dipole moment, frontier molecular orbital (FMO) features, chemical hardness, and the molecular electrostatic potential surface (MEPS) characteristics for optimized molecular geometry of the Green Fluorescent Protein (GFP) chromophore p-hydroxybenzylideneimidazolinone (HBDI) both in its protonated (neutral) and deprotonated (anion) forms. The distribution of atomic charges over the entire molecular framework as obtained from Natural Bond Orbital (NBO) analysis is found to faithfully replicate the predictions from the MEP map in respect of reactivity map of HBDI (neutral and anion) and possible sites for hydrogen bonding interactions etc. The three dimensional MEP map encompassing the entire molecule yields a reliable reactivity map of HBDI molecule also displaying the most probable regions for non-covalent interactions. The differential distribution of the electrostatic potential over the neutral and anionic species of HBDI is authentically reflected on MEP map and NBO charge distribution analysis. Thermodynamic properties such as heat capacity, thermal energy, enthalpy, entropy have been calculated and the correlation of the various thermodynamic functions with temperature has been established for neutral molecule. More importantly, however, the computational approach has been employed to unveil the nonlinear optical (NLO) properties of protonated (neutral) and deprotonated (anion) HBDI. Also in an endeavor to achieve a fuller understanding on this aspect the effect of basis set on the NLO properties of the title molecule has been investigated. Our computations delineate the discernible differences in NLO properties between the neutral and anionic species of HBDI whereby indicating the possibility of development of photoswitchable NLO device.

Seeing through turbidity with harmonic holography [Invited]

The ability to see inside the body noninvasively is indispensable in modern biology and medicine. Optical approaches to such abilities are of rapidly growing interest because of their nonionizing nature and low cost. However, the problem of opacity due to the optical turbidity of tissues must be addressed before optical means become practical. Harmonic holography amalgamates the capability of holographic phase conjugation with the contrast-forming mechanism of second-harmonic generation, which provides a unique opportunity for imaging through a turbid medium. In this review we give accounts of the effort of imaging through turbid media using harmonic holographic phase conjugation.

Effects of the nature and length of the π-conjugated bridge on the second-order nonlinear optical responses of push-pull molecules including 4,5-dicyanoimidazole and their protonated forms

Theoretical and experimental methods of characterization have been employed to investigate a series of six push-pull molecules containing a 4,5-dicyanoimidazole acceptor (A) unit, a N,N-dimethylamino donor (D) group, and systematically enlarged π-conjugated linkers (π) with the view of assessing their first hyperpolarizabilities (β) as well as their variations upon protonation. The results demonstrated that protonation occurred at the N,N-dimethylamino function, which led to disruption of the electronic charge-transfer delocalization: the A-π-D pattern of the neutral species with large β values was transformed into an A-π-A' pattern, which displayed much smaller β values. This feature makes these systems applicable as pH-triggered nonlinear optics (NLO) switches. In particular, protonation led to a decrease of the hyper-Rayleigh scattering first hyperpolarizabilities by at least one order of magnitude as well as to a decrease of the depolarization ratio, from about five for neutral species, which also supported the 1D-like NLO-phore nature, to about two for the protonated ones.