Literature DB >> 35514650

Isoquinoline gas-phase absorption spectrum in the vacuum ultraviolet between 3.7 and 10.7 eV. New valence and Rydberg electronic states.

Sydney Leach1, Nykola C Jones2, Søren V Hoffmann2, Sun Un3.   

Abstract

VUV photons from a synchrotron source were used to record the gas-phase absorption spectrum of isoquinoline over the range 3.5 to 10.7 eV. The rich spectrum exhibits both broad and sharp features, of varying intensities, that are analyzed into eight valence and eight Rydberg transitions. Previous data on the valence transitions of isoquinoline were essentially limited to solution spectra up to 5.4 eV. Our study increases their number considerably. The features in the 3.96 eV region are discussed in terms of vibronic coupling between the nπ* 11A'' and ππ* 21A' valence electronic states. The intensities of some spectral features are augmented by collective π-electron modes considered to be of plasmon-type. Assignments of the valence transitions were facilitated by our DFT calculations and by earlier Pariser-Parr-Pople MO calculations. The calculation results are compared and their relative value is discussed. The DFT calculations reproduce very well a number of experimentally determined properties of the ground state of isoquinoline, in particular its bond distances and angles, rotational constants, vibrational frequencies and dipole moment. No Rydberg series of isoquinoline have previously been observed. Three of the newly observed Rydberg series converge to the D0 electronic ground state of the ion, while two converge to the D1 and three to the D3 excited electronic states of the cation. Astrophysical applications of the VUV absorption spectrum of isoquinoline, in particular the measured absorption cross-sections, are briefly discussed. A comparison between the absorption spectra of isoquinoline and quinoline highlights their similarities and differences, related to their respective molecular orbitals. This journal is © The Royal Society of Chemistry.

Entities:  

Year:  2019        PMID: 35514650      PMCID: PMC9060705          DOI: 10.1039/c8ra09725a

Source DB:  PubMed          Journal:  RSC Adv        ISSN: 2046-2069            Impact factor:   3.361


Introduction

The azaarenes quinoline and isoquinoline, in which a benzene ring is fused to pyridine, are nitrogen analogues of the polycyclic aromatic hydrocarbon (PAH) naphthalene. They can be considered as examples of PANHs, the term now used by astrophysicists for nitrogen analogues of PAHs. Their presence has been suggested, and in some cases observed, in a number of astrophysical sites: planetary atmospheres,[1] the interstellar medium[2-4] and in meteorites.[5,6] The gas-phase absorption spectra of these azaarenes are of use for furthering their astrophysical potential and their observation, in addition to astrochemical modelling studies of their photochemistry in planetary atmospheres, the ISM and circumstellar media. We note that recent laboratory studies[7,8] have shown that quinoline and isoquinoline can be synthesised under high temperature conditions representing circumstellar envelopes of carbon stars. There have been radiofrequency spectroscopic attempts to observe both C9H7N isomers quinoline and isoquinoline in the circumstellar envelopes of carbon-rich stars.[8] Although no positive signals were observed, it has proved possible to derive upper limits on quinoline and isoquinoline column densities. Since the photon absorption of isoquinoline extends in the VUV into energy regions where photoionization processes occur, the measurement of photoionization quantum yields as function of excitation energy will become necessary in order for the absorption cross-sections to be exploited in cosmochemical modeling, although some coarse grain information can be obtained using rule-of-thumb estimations of photon energy dependent ionization quantum yields.[10,11] This information is notably required in studies of photophysical and photochemical processes that involve superexcited electronic states of neutral molecules.[12] The VUV absorption spectrum of gas-phase quinoline has recently been studied.[13] The present work is a study of the VUV absorption spectrum of its isomer, gas-phase isoquinoline. A comparison between the results of the quinoline and isoquinoline spectral studies is given in a later section (Section 4) which includes a discussion on astrophysical applications. Previous studies of the absorption spectra of isoquinoline have been carried out mainly in solution or in matrices in the near-UV region. The fluorescence and phosphorescence spectra of isoquinoline have been studied in detail.[14] Previous gas phase studies of the electronic transitions of isoquinoline[15-23] are limited to the 3.8–4.5 eV region. They are mainly concerned with the interaction of low-lying quasi-degenerate n–π* 11A′′ and π–π* 11A′ states and with related internal conversion and intersystem crossing behavior. Absorption spectra of isoquinoline in solution have an upper energy limit of 5.4 eV.[14,24,25] The synchrotron radiation source provided photons that enabled us to extend absorption spectra of gas phase isoquinoline up to 10.7 eV. Our spectral analysis identified new valence transitions, and assigned eight Rydberg series of bands converging to several electronic states of the cation. This study also provided data on the UV and VUV absorption cross sections of isoquinoline that are potentially of use of in modeling cosmochemical processes.

Experimental

Absorption spectra measurements

Absorption measurements on isoquinoline were performed on the AU-UV beamline stationed in the ASTRID2 synchrotron facility at Aarhus, Denmark. The beamline, its resolution possibilities, and procedure for determining absorption cross sections are described in detail elsewhere.[26,27] Isoquinoline was introduced at room temperature into a gas cell whose path length (l) is 15.5 cm. The entrance and exit windows of the gas cell are in MgF2, VUV transparent to about 10.8 eV (115 nm). VUV radiation was made available at selected energies by the use of a toroidal dispersion grating. Wavelength scanning was achieved with step sizes mainly ranging from 0.05 to 0.2 nm that were chosen in accord with the degree of spectral feature structure in various regions of the spectrum. Some measurements were made over a limited spectral range with 0.02 nm steps. Photon resolution is 0.08 nm, which corresponds to 2 meV at the mid-point of the wavelength region measured. A UV-VUV sensitive photomultiplier (ET Enterprises 9406B) was used to measure the intensity of light transmitted through the gas (It). A small gap between the exit window of the gas cell and the recording photomultiplier is evacuated so as to maintain transparency at photon energies above 5.5 eV by preventing air contamination. For the energy range below 5.5 eV higher order radiation is avoided by letting air into the gap. The MgF2 entrance window of the cell eliminates higher order radiation of energy above 11 eV. Measurement of the gas pressure is made by a heated capacitance manometer (Inficon CDG100D). The pressure used for the spectral measurements ranged from 0.025 to 0.075 mbar. The maximum vapour pressure of the sample at room temperature was indeed 0.075 mbar. The pressure was chosen so as to keep the average attenuation to within a small range around 50%, in order to avoid saturation effects. The intensity of the synchrotron beam current is measured during each wavelength scan, and the measured spectral intensity is normalised to the beam current. A background scan is recorded with an evacuated cell so as to provide the I0 intensity. Wavelength calibration was carried out by comparison with well-known standard absorption spectra, as follows: nitrogen between 120 and 150 nm; two SO2 bands in the mid (185–225 nm) and long wavelength range (275–315 nm). The measured absorption cross sections are accurate to ± 5%. They were obtained by applying the Beer–Lambert law: I = I0 exp(nσl). The highest purity commercially available (Sigma-Aldrich) isoquinoline sample underwent 4 freeze–pump–thaw cycles and received no further treatment before use.

Computations: methods and results

Calculations of the valence transitions of isoquinoline (Cs symmetry group structure, Fig. 1) were carried out to aid in spectral transition assignments, in particular of the π–π* transitions. Previous reported computations of these transitions are predominantly by Pariser–Parr–Pople (P–P–P) quantum mechanical methods. We carried out time-dependent DFT calculations of structural and spectroscopic properties of isoquinoline for comparison with previous theoretical studies. We note that the TD-DFT formalism, as derived from Kohn–Sham theory, is essentially an empirically corrected version of the random-phase approximation used in Hartree–Fock methods.[28,29]
Fig. 1

Isoquinoline ground state bond distances (Å) and angles (degrees) calculated by a time-dependent DFT method. (see text).

We carried quantum chemical calculations with the GAUSSIAN 16 program, Revision A.03.[30] Geometry optimization of isoquinoline was performed using two different basis sets: the hybrid density functional B3LYP and 6-31+G(D,P) basis set (calculation A) and the B971/6-311 +G(D,P) combination (calculation B). Normal mode analysis confirmed that the final optimized geometries were true potential energy minima. The geometric centres of calculation A and calculation B isoquinoline molecules are the same. The average difference between calculations A and B in the coordinates of the atoms with respect to this center is 0.005 Å with a maximum of 0.007 Å. This suggests that there are no significant differences in the accuracies between the two calculation methods. The calculation B structure was used to carry out TD-DFT calculations of electronic state energies and dipole transition oscillator strengths using the CAM-B3LYP hybrid density-functional[31] along with 6-31+G (d,p) as basis-set.

Geometry and other structural properties

The geometry of isolated isoquinoline has not been measured directly. Our calculated bond distances and bond angles (Fig. 1) are in very reasonable agreement with other calculated values.[32,33] These bond distances and angles have, however, been measured for crystalline isoquinoline in an X-ray crystal structure study;[34] we note that the structure is complicated by an inversion symmetry in the crystal. The atom–atom bond distances obtained with the calculation A and calculation B functionals have an average deviation from the X-ray determined structures generally less than 0.01 Å, the calculation B values being slightly better than those of calculation A. The calculated rotational constants, based on our calculated isoquinoline structures, are in quite good agreement with the experimental values (Table 1) determined by high-resolution microwave and/or optical spectroscopy measurements on ground state isoquinoline.[33,35-37] Another criterion of quality of our calculations is the dipole moment of isoquinoline. Our values agree well with the dipole moment determined by gas-phase polarization measurements,[38] although they are about 9% greater than Stark effect values (Table 1). Polarization measurements of the dipole moment in the gas phase necessitate measurement over a temperature range of at least 100 K, as achieved by Buckingham et al.[38] Stark effect measurements of gas-phase dipole moments are subject to uncertainties that are discussed by Kisiel et al.[33]

Isoquinoline ground state rotational constants and dipole moment: calculated and experimental values

A Rotational constants
(a) DFT calculated: (present study)
A = 3.1921294 GHz, B = 1.2330017 GHz, C = 0.8894428 GHz
(b) Microwave spectroscopy[33]
A = 3.19900020 GHz, B = 1.237931586 GHz, C = 0.892753595 GHz
(c) FIR spectroscopy[36]
Analysis based on ν30. A = 3.2003697 GHz, B = 1.2387657 GHz, C = 0.8937165 GHz
(d) FIR spectroscopy[37]
Analysis based on ν37. A = 3.1980719 GHz, B = 1.23780915 GHz, C = 0.8928502 GHz
B Dipole moment
(a) DFT calculated: 2.7536 Debye (present study)
(b) Polarisation measurement[38]: 2.73 ± 0.03 Debye
(c) Stark effect measurement[33]: 2.5277 Debye

Vibrational frequencies of the ground state

Table 2 presents the ground state harmonic and anharmonic vibrational frequencies of the thirty-one a′ and fourteen a′′ vibrational modes of isoquinoline calculated by our DFT method. A comparison is made with experimental values[39] and with those of other calculations. These provide comparative tests of the various calculation methods of vibrational frequencies. The experimental frequencies of Wait and McNerney[39] are from Raman and infra-red spectroscopic studies in either liquid or gas phase. To make their mode assignments, Wait and McNerney carried out a normal coordinate analysis using the Urey-Bradley Force-Field, including Kekulé modification, using force constants that were derived from naphthalene and s-triazine. More recently, other information, in particular on the frequencies of modes 28, 31, 40, 41, 43 and 45 of isoquinoline was obtained by Martin-Drumel et al.[35] in high resolution gas phase far-infrared spectroscopic studies using a Fourier transform instrument.

Isoquinoline vibrational frequencies (cm−1): harmonic and anharmonic frequency calculation results compared with observed values

Vibrational mode[39]Observed[39]DFT [A] harmon present studyDFT [B] harm present studyDFT harmon[35]DFT [A] anharmon present studyDFT [B] anharmon present studyDFT anharmon[35]Hartree–Fock[32]
1 a′30893189318731893066305730513374
2 a′30593179317631773034303130433372
3 a′30553177317531763070305930483361
4 a′30553166316331643046304830493352
5 a′30243161315831593027302230223347
6 a′30073158315631552998299930053340
7 a′29943124312231203003299329963321
8 a′16251664165816621673161716211846
9 a′15871623161716211586158115851802
10 a′15561605160116041574156915661785
11 a′14971533152815301499149314951671
12 a′14611489148314841457145614591625
13 a′14321464145714591433142614291593
14 a′13821410140514061383137813741525
15 a′13771400139613971364135913671495
16 a′13151364136213631333133113331471
17 a′12731290128312841262125512561388
18 a′12531278127012701256125012511378
19 a′11781241123812391219121412171331
20 a′11401205120012011186118111821293
21 a′11191168116311631150114411451255
22 a′10951160115411551144113811391192
23 a′10341062105910601045104110411135
24 a′10141036103210331018101510151081
25 a′9599579509509449379371030
26 a′778812807807802795796866
27 a′765790786787779775776851
28 a′637/640 (ref. 35)656650651649642644698
29 a′522/521 (ref. 35)530527527525521522566
30 a′504/502 (ref. 35)511507507506501501544
31 a′381362358358360355355386
32 a′′98410059989999829769761127
33 a′′9719939879889779749721118
34 a′′9459769709689589549541096
35 a′′9309469409479389369311067
36 a′′862877872874864862862973
37 a′′825844841845830829830935
38 a′′800789782791782779780863
39 a′′742753751752741739739837
40 a′′611/635 (ref. 35)652648652640637641705
41 a′′482/480 (ref. 35)493489491484481483542
42 a′′459/457 (ref. 35)469465468461460461518
43 a′′355/376 (ref. 35)386384384378377378431
44 a′′201/181 (ref. 35)184184185178180182207
45 a′′182/166 (ref. 35)171171171168161168189
Slope1.047 ± 0.0041.048 ± 0.0041.048 ± 0.0040.994 ± 0.0040.994 ± 0.0040.994 ± 0.0041.105 ± 0.005
Intercept−21 ± 6−27 ± 7−26 ± 716 ± 79 ± 611 ± 63 ± 10
Correlation coefficient0.999810.999800.999800.999750.999790.999800.99961
Included in Table 2 are the vibrational frequencies of isoquinoline calculated in a force-field study by Bandyopadhyay and Manogaran[32] with a geometry obtained with Hartree–Fock theory using a 6-31G** basis set. To reconcile the ab initio values of the harmonic vibrational frequencies in the table with experimental values necessitates 0.7–1.0 scaling for the a′ vibrational modes and 0.5–1.0 for the a′′ modes. Scott and Radom[40] have compared a number of ab initio, semiempirical and DFT calculations of molecular harmonic vibrational frequencies so as to propose useful scaling values. The results presented in Table 2 show that our calculation A and calculation B anharmonic frequencies, which were obtained using the second-order vibrational perturbation theory devised by Barone,[41] in Gaussian 16, require little scaling. There is generally less than 2% deviation of the anharmonic frequencies from observed. They are mainly nearer in value than those calculated by Martin-Drumel et al.[35] with the hybrid B97-1 functional and the 6-311G(d,p) basis set. One can rationalize this as essentially resulting from the different characteristics of second-order vibrational perturbation theory used in the various used versions of Gaussian. A graphical comparison of the various calculation results is presented in Fig. 2, with the relevant statistics given in Table 2.
Fig. 2

Comparison between measured and calculated isoquinoline ground state vibrational frequencies.

Examination of Table 2 and Fig. 2 shows that, as compared with the Hartree–Fock calculation of Bandyopadhyay and Manogaran,[32] the values of the DFT calculations necessitate far less scaling or refinement techniques to come into line with observed values. This is consistent with the conclusions of the comparative study of Scott and Radom.[40]

Electronic state energies and transition oscillator strengths

The calculation B structure was used in the time-dependent DFT method[28] calculation of the transition energies and oscillator strengths involving 161 singlet electronic excited states of isoquinoline up to 11.457 eV (sixty-eight A′ states, eighty-one A′′ states). Since the eigenvalues depend critically on M.O. energy differences in time-dependent DFT calculations these may not provide reliable values for transition energies. The upper energy limit of our calculation was chosen to be greater than our observational limit of 10.7 eV in order to accommodate inaccuracies in the calculated energies. We stress, however, that the DFT calculations allow us to identify strong and weak electronic transitions, since oscillator strengths calculated by DFT methods depend little on the precision of calculated transition energies.[42,43] We note also that the true oscillator strengths can be modified if there exist plasmon-type interaction effects (see later). Table 3 lists the valence transition energies and oscillator strengths calculated by the DFT method. In the energy range up to 11.457 eV there are 33 allowed 1A′ ← 1A′ transitions whose oscillator strength is calculated to be greater than f = 0.01, but only six greater than f = 0.1.

Isoquinoline. DFT calculated singlet valence transition energies and oscillator strengths

TransitionTransition energy eVOscillator strength f
2A′ ← 1A′4.54470.0577
1A′′ ← 1A′4.69450.0020
3A′ ← 1A′4.86060.0340
2A′′ ← 1A′5.54630.0012
4A′ ← 1A′6.10180.6111
3A′′ ← 1A′6.20210.0015
5A′ ← 1A′6.20990.5681
6A′ ← 1A′6.39920.1484
4A′′ ← 1A′6.54180.0004
7A′ ← 1A′6.56720.1306
8A′ ← 1A′6.73280.0186
5A′′ ← 1A′6.75120.0083
6A′′ ← 1A′7.06960.0006
7A′′ ← 1A′7.11850.0001
8A′′ ← 1A′7.14860.0184
9A′ ← 1A′7.18680.0143
9A′′ ← 1A′7.22190.0000
10A′′ ← 1A′7.36240.0055
11A′′ ← 1A′7.47530.0079
10A′ ← 1A′7.55490.0028
11A′ ← 1A′7.67160.0063
12A′′ ← 1A′7.68580.0002
12A′ ← 1A′7.69670.0416
13A′′ ← 1A′7.74320.0006
14A′′ ← 1A′7.79560.0002
15A′′ ← 1A′7.86390.0033
13A′ ← 1A′7.93570.0106
14A′ ← 1A′8.03640.0049
16A′′ ← 1A′8.05280.0022
17A′′ ← 1A′8.08960.0004
15A′ ← 1A′8.12070.1941
16A′ ← 1A′8.25070.0407
18A′′ ← 1A′8.26730.0029
19A′′ ← 1A′8.29800.0001
20A′′ ← 1A′8.35000.0001
17A′ ← 1A′8.37100.1121
18A′ ← 1A′8.38330.0758
21A′′ ← 1A′8.39740.0111
22A′′ ← 1A′8.40800.0006
23A′′ ← 1A′8.47960.0080
19A′ ← 1A′8.49750.0037
20A′ ← 1A′8.51120.0175
24A′′ ← 1A′8.53200.0017
21A′ ← 1A′8.59590.0091
25A′′ ← 1A′8.60720.0015
26A′′ ← 1A′8.73680.0176
22A′ ← 1A′8.77020.0021
23A′ ← 1A′8.88610.0006
27A′′ ← 1A′8.95660.0027
28A′′ ← 1A′8.97630.0001
24A′ ← 1A′8.99920.0020
29A′′ ← 1A′9.10040.0012
30A′′ ← 1A′9.11700.0012
31A′′ ← 1A′9.13340.0030
32A′′ ← 1A′9.17170.0001
25A′ ← 1A′9.18070.0032
33A′′ ← 1A′9.21300.0001
26A′ ← 1A′9.22950.0030
27A′ ← 1A′9.24000.0523
34A′′ ← 1A′9.26510.0020
35A′′ ← 1A′9.30870.0006
36A′′ ← 1A′9.35160.0004
28A′ ← 1A′9.36750.0029
37A′′ ← 1A′9.42440.0004
38A′′ ← 1A′9.43120.0012
29A′ ← 1A′9.44330.0031
39A′′ ← 1A′9.47710.0000
30A′ ← 1A′9.50550.0011
40A′′ ← 1A′9.55160.0215
31A′ ← 1A′9.56380.0018
32A′ ← 1A′9.57780.0285
41A′′ ← 1A′9.59830.0026
42A′′ ← 1A′9.63230.0010
43A′′ ← 1A′9.75020.0000
33A′ ← 1A′9.75520.0535
44A′′ ← 1A′9.77870.0007
34A′ ← 1A′9.80380.0050
35A′ ← 1A′9.89220.0818
36A′ ← 1A′9.92240.0129
37A′ ← 1A′10.02290.0233
38A′ ← 1A′10.03040.0097
45A′′ ← 1A′10.04020.0068
46A′′ ← 1A′10.05620.0023
47A′′ ← 1A′10.11340.0062
48A′′ ← 1A′10.12120.0001
39A′ ← 1A′10.14270.0005
49A′′ ← 1A′10.17850.0101
40A′ ← 1A′10.18650.0883
50A′′ ← 1A′10.19570.0050
41A′ ← 1A′10.21370.0029
42A′ ← 1A′10.25490.0010
51A′′ ← 1A′10.27090.0132
52A′′ ← 1A′10.30080.0000
43A′ ← 1A′10.30850.0195
53A′′ ← 1A′10.31570.0000
54A′′ ← 1A′10.35530.0035
44A′ ← 1A′10.36840.0027
45A′ ← 1A′10.40170.0010
55A′′ ← 1A′10.43050.0432
46A′ ← 1A′10.43670.0032
47A′ ← 1A′10.48390.0072
56A′′ ← 1A′10.48750.0017
57A′′ ← 1A′10.49950.0040
48A′ ← 1A′10.52980.0022
49A′ ← 1A′10.55220.0003
58A′′ ← 1A′10.56410.0000
59A′′ ← 1A′10.577770.0010
60A′′ ← 1A′10.62650.0071
50A′ ← 1A′10.65010.0037
61A′′ ← 1A′10.65990.0002
62A′′ ← 1A′10.69750.0664
63A′′ ← 1A′10.71330.0212
51A′ ← 1A′10.71700.0308
64A′′ ← 1A′10.75140.1542
52A′ ← 1A′10.75700.0312
53A′ ← 1A′10.77930.0376
54A′ ← 1A′10.84080.0034
65A′′ ← 1A′10.85930.0076
66A′′ ← 1A′10.89450.0002
67A′′ ← 1A′10.89890.0111
55A′ ← 1A′10.92730.0102
68A′′ ← 1A′10.95680.0013
56A′ ← 1A′10.97180.0018
69A′′ ← 1A′11.01810.0000
70A′′ ← 1A′11.11840.0330
71A′′ ← 1A′11.12200.0016
72A′′ ← 1A′11.17200.0554
57A′ ← 1A′11.19880.0954
58A′ ← 1A′11.22240.0212
73A′′ ← 1A′11.22830.0006
74A′′ ← 1A′11.24800.0051
75A′′ ← 1A′11.26790.0447
59A′ ← 1A′11.27530.1237
76A′′ ← 1A′11.28180.0161
60A′ ← 1A′11.29680.0186
77A′′ ← 1A′11.31490.0070
78A′′ ← 1A′11.33820.0178
61A′ ← 1A′11.38180.0123
79A′′ ← 1A′11.39320.0427
80A′′ ← 1A′11.40470.0027
62A′ ← 1A′11.41780.0667
63A′ ← 1A′11.42440.0267
81A′′ ← 1A′11.45700.0003

Results and discussion

The absorption spectrum: general remarks

The absorption spectrum of gas phase isoquinoline between 3.6 and 10.7 eV, shown in Fig. 3, includes as an inset a 15× intensity magnification of the region between 3.8 and 5.6 eV. This spectrum contains a great number of features, of which many are small peaks that overlie diffuse continua. A few of the measured features are incompletely resolved shoulders in the peaks. The 161 measured features are seen in greater detail in intensity zooms of a number of parts of the absorption spectrum. We give typical examples that respectively show the spectral regions 3.7–4.5 eV (Fig. 4), 4.5–5.8 eV (Fig. 5), 5.5–6.6 eV (Fig. 6), 6.3–8.5 eV (Fig. 7) and 8.5–10.7 eV (Fig. 8). The band energies reported in Table 4 are recorded for clear peaks forming maxima in the recorded absorption cross-sections and, in the case of less well-resolved features, as gradient particularities in the absorption profile. The precision of peak energy pointing, between ±0.0006 eV and ±0.004 eV, is a function of the spectral energy range and the step size used in wavelength scanning. The measured features are listed in Table 4 as band numbers, band energies and band assignments. The assignments of the valence transitions have been made with the help of quantum mechanical calculations of electronic state energies and transition strengths, as discussed below.
Fig. 3

Absorption spectrum of isoquinoline between 3.6 and 10.8 eV. Rydberg series bands are marked. The origin bands O00 of the set of π–π* valence transitions 21A′ ←11A′ up to 91A′ ←11A′ (Table 4) are indicated as 21A′…. 91A′. Vibrational components of vibronic transitions listed in Table 4 are not denoted, for figure clarity.

Fig. 4

Absorption spectrum of isoquinoline between 3.7 and 4.5 eV. Band numbers are pinpointed.

Fig. 5

Absorption spectrum of isoquinoline between 4.5 and 5.8 eV. Band numbers are pinpointed.

Fig. 6

Absorption spectrum of isoquinoline between 5.5 and 6.6 eV. Band numbers are indicated.

Fig. 7

Absorption spectrum of isoquinoline between 6.3 and 8.5 eV. Band numbers are pinpointed.

Fig. 8

Absorption spectrum of isoquinoline between 8.5 and 10.7 eV. Band numbers are pinpointed.

Isoquinoline gas-phase absorption spectrum: band energies and assignments

Band no.Energy/eVEnergy/cm−1Assignment
13.89331399Hot band (see text)
23.89831439Hot band 500 cm−1 (see text)
33.90531496Hot band (see text)
43.94431810s33 (sequence, see text)
53.94931851s22 (sequence, see text)
63.95431891s11 (sequence, see text)
73.96031939π–π* 21A′ ← 11A′O00: 0+(ππ*) (principal component)
83.96331964
93.96731996a10s11 (see text)
103.970320200: a10(nπ*) (principal component)
113.980321010=: b10(nπ*) (principal component)
123.98932173
133.9943221721A′ ← 11A′φ20
144.01432375
154.01832407
164.0233244821A′ ← 11A′α10
174.03232520
184.03632552
194.03932577
204.025 (sh)32464
214.04532625
224.0503266521A′ ← 11A′β10
234.07232843
244.08032907
254.0843294021A′ ← 11A′α20
264.09533028
274.1133317321A′ ← 11A′α10β10
284.1343334321A′ ← 11A′γ10
294.1433341521A′ ← 11A′α30
304.16133561
314.16833617
324.1753367321A′ ← 11A′α20β10
334.19433827
344.1973385121A′ ← 11A′ γ10α10
354.20433907
364.2093394821A′ ← 11A′γ10β10
374.21834020
384.22434069
394.2583434321A′ ← 11A′γ10α20
404.2673441521A′ ← 11A′α20β20
414.3063473021A′ ← 11A′γ20
424.33634972
434.34035004
444.34535045
454.3963545621A′ ← 11A′γ20β10
464.479 (sh)36125
474.56436814π–π* 31A′ ← 11A′O00
484.6273731631A′ ← 11A′α10
494.7363820431A′ ← 11A′γ10
504.78138561
514.91339626
524.95639973
535.13341400
545.31142836
555.36843296
565.48944272
575.55944836
585.68145822π–π* 41A′ ← 11A′O00
595.71846120π–π* 51A′ ← 11A′O00
605.7714654641A′ ← 11A′β10
615.86447296R1 (n = 3) O00, δ = 0.7409
625.902 (sh)47603R2 (n = 3) O00, δ = 0.724
636.00348417
646.03848699π–π* 61A′ ← 11A′O00
656.09849183
666.1394951461A′ ← 11A′β10
676.21350111
686.25150417
696.32250990π–π* 71A′ ← 11A′O00; R4 (n = 3) O00, δ = 0.8142
706.69353982R3 (n = 3) O00, δ = 0.28
716.75554482R5 (n = 3) O00, δ = 0.621
726.86755386π–π* 81A′ ← 11A′O00
736.9635616081A′ ← 11A′β10
747.08157112π–π* 91A′ ← 11A′O00
757.11757402R1 (n = 4) O00, δ = 0.8969
767.1755787091A′ ← 11A′β10
77 7.261 58564
787.28058717R1 (n = 4) z10; R2 (n = 4) O00, δ = 0.70
797.36059362
807.45360112R6 (n = 3) O00, δ = 0.851
817.52660701R3 (n = 4) O00, δ = 0.32
827.62261475
837.68761999
847.71562225R1 (n = 5) O00, δ = 0.9141
857.78562790R2 (n = 5) O00, δ = 0.726
867.84363258R4 (n = 4) O00, δ = 0.8183; R7 (n = 3) O00, δ = 0.693
877.88063556R1 (n = 5) z10
887.90463750
897.92063879R3 (n = 5) O00, δ = 0.28
907.93063959
917.98664411R1 (n = 6) O00, δ = 0.9989
927.99964516R5 (n = 4) O00, δ = 0.577
938.03864830R2 (n = 6) O00, δ = 0.741
948.12065492R3 (n = 6) O00, δ = 0.24
958.15465766R1 (n = 7) O00, δ = 0.9845
968.19766113R2 (n = 7) O00, δ = 0.608
978.25266557R1 (n = 8) O00, δ = 1.004
988.28166790R2 (n = 8) O00, δ = 0.608
998.31867089R1 (n = 9) O00, δ = 0.9888
1008.35267364R4 (n = 5) O00, δ = 0.8965
1018.38367613
1028.39767726R1 (n = 11) O00, δ = 0.8856
1038.42067912R1 (n = 12) O00, δ = 0.8784
1048.44068073R5 (n = 5) O00, δ = 0.653
1058.58969275R8 (n = 3) O00, δ = 0.259
1068.64069686R4 (n = 6) O00, δ = 0.8848
1078.69070089R5 (n = 6) O00, δ = 0.620
1088.77170743R4 (n = 7) O00, δ = 1.0859
1098.95872251
1108.96572307
1118.97572388
1128.99172517
1139.00772646
1149.03172840R6 (n = 4) O00, δ = 0.847
1159.17473993R7 (n = 4) O00, δ = 0.669
1169.23974517
1179.24674574
1189.31275103
1199.42376004R8 (n = 4) O00, δ = 0.268
1209.49776598
1219.52376808R6 (n = 5) O00, δ = 1.061
1229.53476896
1239.55277042
1249.60777485
1259.66077913
1269.66777969R7 (n = 5) O00, δ = 0.692
1279.68478106
1289.70178243
1299.70978308
1309.72078397
1319.72878461
1329.73678526
1339.77478832R8 (n = 5) O00, δ = 0.338
1349.86079526R6 (n = 6) O00, δ = 0.9804
1359.88779744
1369.91980002R7 (n = 6) O00, δ = 0.68
1379.95980324
1389.96780389
1399.99980647R8 (n = 6) O00, δ = 0.175
14010.03180905R6 (n = 7) O00, δ = 0.9277
14110.07281236R7 (n = 7) O00, δ = 0.559
14210.09681429R8 (n = 7) O00, δ = 0.310
14310.10581502
14410.14081784R6 (n = 8) O00, δ = 0.7937; R7 (n = 8) O00, δ = 0.766
14510.20482300R6 unresolved n = 9, n = 10; R7 (n = 9) O00, δ = 0.668
14610.26482784R6 unresolved n > 10 values
14710.27282849
14810.39783857
14910.41984034
15010.47284462
15110.49884672
15210.53884994
15310.54785067
15410.55685139
15510.56585212
15610.57485285
15710.58385357
15810.59785470
15910.61185583
16010.62085656
16110.65685946
A review by Innes et al.[14] lists the triplet and singlet state energies and oscillator strengths of n–π* and π–π* transitions of isoquinoline, culled from experimental and theoretical studies available in 1988. Data on the π–π* transitions were published earlier by Baba and Yamzaki[44] in their Pariser–Parr–Pople SCF LCAO MO method calculation study of the absorption spectrum of isoquinoline. Their calculation yielded eight π–π* transitions up to 7.13 eV, as did the SCF MO calculations of Nishimoto,[45] who used a fixed core approximation. Eight calculated π–π* transitions are reported in this energy region by Tinland,[46] who used the P–P–P SCF LCAO MO method, and by Ridley and Zerner[47] whose SCF LCAO MO calculations employed a modified INDO method (Table 5).

Isoquinoline π–π* mA′ ← 1A′ electronic singlet state transition energies E (eV) and oscillator strengths (f). Observed and calculated SCF-LCAO-MO and DFT values

Excited state mA′a[24]b[49]c[46]d[45]e[47]f[48]g[44] (Baba)h DFT present studyi Obs. present studyh–i unsignedg–i unsigned
24.5003.8764.0043.974.034.114.064.54473.9600.58470.100
0.026 1Lb 0.040 0.051 0.09 0.01 0.118 0.0577
35.6715.2874.4974.524.614.834.424.86064.5640.29660.144
0.189 1La 0.194 0.180 0.16 0.24 0.147 0.0340
46.2505.5275.315.555.346.10185.6810.42080.331
1Bb 0.207 0.059 0.12 0.020 0.6111
55.6185.505.795.775.646.20995.7180.49190.078
0.831 1.132 0.89 0.02 0.489 0.5681
65.7725.665.895.796.39926.0380.36120.248
0.836 0.257 0.81 0.382 0.1484
76.0685.806.355.986.56726.3220.24520.342
0.070 0.475 0.52 1.151 0.1306
86.2586.186.486.306.73266.8670.13440.567
0.519 0.568 0.08 0.540 0.0186
96.2296.966.847.137.18687.0810.10580.049
0.149 0.001 0.01 0.001 0.0143
107.5549
0.0028
117.6716
0.0063
127.6967
0.0416
137.9357
0.0106
148.0364
0.0049
158.1207
0.1941
168.2507
0.0407
178.3710
0.1121
188.3833
0.0758
198.4975
0.0037
208.5112
0.0175
Calculations that included configuration interactions were carried out in some studies.[44,47,48] On the other hand, as listed in Table 5, more limited π–π* transition calculations were done by Favini et al.,[24] by Mataga[49] and by Wagner et al.[48] Goodman and Harrell[50] made analogous calculations on the n–π* transitions of isoquinoline. Their results and those of our own DFT calculations (Table 5), have helped in assigning the valence transitions of isoquinoline in our gas phase spectrum. Rydberg transition bands of isoquinoline were observed for the first time in our study. There are no reported calculations on Rydberg transitions of isoquinoline. Our analysis and discussion of the π–π* transitions of isoquinoline analysis are based substantially on the calculation results of Baba and Yamakazi.[44] We note that their transition energies and oscillator strengths agree well with those of Nishimoto[45] and are in satisfactory agreement with the values reported by Tinland[46] and by Ridley and Zerner.[47] We remark that the valence transition energies calculated by the DFT method are greater than those of the older, mainly P–P–P, calculations. Our DFT calculations indicate that only a few of the n–π* transitions (Table 3) have oscillator strengths (f ∼ 0.01) possibly of sufficient value to be observable in our spectra (see below).

Absorption 3.8–4.5 eV: 21A′ ← 11A′ π–π* transition origin band region. Vibronic coupling with the 11A′′ nπ* state

The spectral features in this energy region (Fig. 3 and 4) are sharp up to about 4.2 eV. They become more diffuse at higher energies. The maximum absorption cross-section between 3.7 and 5.2 eV is 15 Mb at about 4.8 eV. This region has previously been well studied in solution and in the vapor phase. It contains both n–π* and π–π* transitions. The gas-phase studies of Hiraya et al.,[15,16] assign a very weak n–π* transition, 11A′′ ← 11A′, that appears close to the initial π–π* spectral region. They report the origin band of the 11A′′ ← 11A′ transition to be at 30 821 cm−1 (3.821 eV). This closeness in energy is consistent with our DFT calculations, whose calculated energy separation is 1208 cm−1 (Table 3) as compared with the experimental value 1118 cm−1 (see below), although the calculated order of the first n–π* and π–π* transition energies is reversed with respect to the experimental values (see below). Our calculations give oscillator strengths of f = 0.0020 and 0.0577 respectively for the 11A′′ ← 11A′ n–π* and the 21A′ ← 11A′ π–π* transitions, consistent with their observed relative strengths.[15,16] The weakness of the n–π* transition origin band makes it unobservable in our absorption spectrum of isoquinoline. We note, however, that in the TD-DFT calculations (Table 3), of the 64 n–π* transitions calculated to occur up to 10.8 eV, only ten of them are calculated to have an oscillator strength greater than f = 0.01, (and only one greater than f = 0.1, the 64A′′ ← 1A′ transition at 10.7514 eV, f = 0.1542), are possibly adequate for observation in our spectra. However, these would occur at energies above 7 eV (Table 3) in spectral regions of intense valence π–π* transitions and Rydberg bands. It is known that the order of the first n–π* and π–π* transition energies is reversed in mixed crystals with naphthalene or durene hosts (cited by Hassan and Hollas,[51]) and in solvents where H-bonded complexes can form.[52] This reversed order is given in our DFT calculations on isoquinoline (Table 3). However, in the gas phase, the 11A′′ ← 11A′ n–π* transition origin is reported to lie at 30 821 cm−1,[15,16] which is about 1118 cm−1 below the strongest of a trio of close-lying features seen in the high resolution gas phase spectrum of isoquinoline.[15,16] These three spectral features, traditionally denoted as 0+, 0− and 0= respectively, are in the π–π* 21A′ ← 11A′ origin band region, and they correspond to features in the 3.960 eV band region in our study (Tables 4 and 6). They have been shown to be due to near-resonance vibronic coupling between the n–π* 11A′′ and π–π* 21A′ vibronic levels,[15,16,18] the coupling being induced by out-of-plane ν(a) and ν(b) vibrations of a′′ symmetry in the 11A′′ state.[51]

Isoquinoline absorption origin region

Band no.Absorption cross-section Mb × 10 [relative int.]Energy/eVEnergy/cm−1 [ν – 31939]Energy/cm−1 [relative int.] Fischer and Naaman[18]
112.70 [0.06]3.89331399 [541]31384 (−14)
215.57 [0.08]3.89831439 [−500]31426 (−13)
317.02 [0.08]3.90531496 [−443]31481 (−15)
480.30 [0.40]3.94431810 [−129]31797 (−13) [0.14]
5106.71 [0.53]3.94931851 [−88]31841 (−10) [0.27]
6151.87 [0.75]3.95431891 [−48]31885 (−6) [0.55]
7201.41 [1.00]3.96031939 [0]31925 (−14) [1.00]
870.48 [0.35]3.96331964 [+25]31960 (−4) [0.27]
9102.12 [0.51]3.96731996 [+57]31983 (−13) [0.33]
1038.52 [0.19]3.97032020 [+81]32008 (−12) [0.05]
1122.01 [0.11]3.98032101 [+162]*
1224.44 [0.12]3.98932173 [+234]*
1327.72 [0.14]3.99432214 [+275]32203 (−11)
The relative proportion of the 21A′ state varies across the envelope of the origin band region measured by Hiraya et al,[15] as can be seen in a more refined way in the rotational contour analysis of a high resolution gas phase spectrum of isoquinoline absorption in the 31 925 cm−1 region.[51] The assignments 0+(ππ*), 0− a10(nπ*) and 0= b10(nπ*) in Tables 4 and 6 refer to the principal zero-order vibronic components in the observed band (see below). The spectral effects of vibronic coupling between these non-degenerate but close-lying two states have been studied in detail by Fischer and co-workers.[18,19,53] We recall that the isoquinoline n–π* 11A′′ level is shifted to higher energies, above the π–π* 21A′ level, in an isoquinoline hydrogen-bonded complex with methyl alcohol in a supersonic jet.[52] The π–π* transition origin band has a single feature, at 31 880 cm−1 (3.953 eV) in this case, indicating that the π–π* 21A′ level is no longer perturbed by the n–π* 11A′′ state. The features between 3.893 and 3.994 eV, severely marked by vibronic interactions, are also presented in Table 6, where we compare related gas-phase spectral observations of this complex spectral region. The data of Fischer and Naaman[18] are from absorption spectra recorded on photographic plates. Their reported band energies, used without confirmation in other isoquinoline studies[19,54] are consistently about 10 cm−1 smaller than our measurements. The reported relative intensities of the bands follow the relative trends of our absorption cross section measurements but there are some quantitative differences, due no doubt to the well-known difficulty of measuring band intensities on photographic plates. We remark that the relative intensities of the three main bands, denoted traditionally as 0+, 0− and 0=, reported as 1.00 : 0.63 : 0.20 by Hiraya et al.,[15] based on a jet-cooled multiphoton ionization photoelectron study, are in excellent agreement with those obtained by our absorption cross-section measurements, 1.00 : 0.51 : 0.19 (Table 6). These features have also been shown not to be hot bands.[18] Following the assignments of Fischer and Naaman,[18] bands numbered 1, 2 and 3 are true hot bands that are three members of a sequence progression with a sequence interval ν′ − ν′′ ≈ 40 cm−1. Calculations indicate that in the 2 1A′ and 11A′′ states the sequence forming vibrations have frequencies 135.5 cm−1 and ∼175.5 cm−1,[18] (139 cm−1 and 182 cm−1),[19] respectively. The band number 13 at 3.994 eV can be assigned to excitation of two quanta of the sequence-forming mode ν′(φ) ∼139 cm−1, a mode possibly related to ν′′(44) = 181 cm−1 or ν′′(45) = 166 cm−1 that are out-of-plane vibrations in ground state isoquinoline (Table 2). Band number 7 is the most intense of the origin band group and is chiefly 1ππ* in character.[18] The assignments of bands numbered 7, 9 and 10, (respectively 0+, 0− and 0=) are the principal components of the products of the resonance between the zero-order levels 0(1ππ*), a(1nπ*) and b(1nπ*), where 0(1ππ*) is the vibrationless level of the 21A′ 1ππ* state and a and b refer to out-of-plane modes active in the 11A′′ 1nπ* state. The ν(a) and ν(b) modes and their frequencies have not been identified[18,19] but their frequencies could be of the order of 1000 cm−1.[51] The coupling parameter via ν(a), 26 cm−1, is bigger than through ν(b), 10 cm−1,[51] so that vibronic absorption band a10 is stronger than band b10, as evidenced in our spectra by band number 10 at 3.970 eV being about twice as intense as band number 11 at 3.980 eV (Tables 4 and 6). Some other bands in this region have assignments involving mode s, involved in the sequence structure.[18,51] The relevant ground and excited state zero-order mode frequencies (before resonant level interactions) are = 181 cm−1i.e. mode ν44 (Table 2), = 135.5 cm−1, = 175.5 cm−1. Bands numbered 4–45 (Table 4) are assigned to vibronic components of the π–π* 21A′ ← 11A′ transition. Four vibrational modes are considered to be active in the excited state: ν′(φ) = 139 cm−1, a non-totally symmetric vibration, active in two quanta as discussed above, ν′(α) = 510 cm−1, ν′(β) = 730 cm−1 and ν(γ) ∼1400 cm−1. Related frequencies in the ground state are ν′′(30) = 502 cm−1, ν′′(26) = 778 cm−1, ν′′(13) = 1432 cm−1, ν′′(14) = 1382 cm−1. The strongest of the bands numbered 4–45 are present in the gas phase absorption spectra of Okajima and Lim[55] and of Keith-Jameson et al.[22] but their energies were not reported. Our band frequencies are consistent with those reported up to 4.174 eV in a gas phase absorption spectrum by Fischer and Naaman.[18] The diffuseness of the vibronic bands at energies above 4.05 eV in all reported gas phase absorption spectra of isoquinoline, was attributed by Fischer and Naaman[18] to the fast relaxation processes of internal conversion between the lowest excited singlet ππ and nπ* electronic states, in addition to singlet-triplet intersystem crossing. On the other hand, Fischer and Knight[19] emphasized the effects of spectral congestion in creating the broadness of the bands. It is probable that both internal conversion and spectral congestion are in play.[20,23] We note, for example, that the density of vibronic states in the 21A′ level at about 4.1 eV, ∼1130 cm−1 above the 0+ level, can be estimated to be of the order of 1 per cm−1. The density of quasi-degenerate vibronic levels in the 11A′′ state would be much higher at that energy, thus ensuring the possibility of much vibronic mixing between the n–π* and π–π* levels. The vibronic bands of the 21A′ ← 11A′ transition are also observed in condensed phase absorption spectra[24,25,56-58] but the spectral resolution is insufficient to investigate possible influences of solvent on isoquinoline vibronic level energies as well as to consider further aspects of mixed n–π* and π–π* levels in the absorption bands. The structured features in the 3.8–4.4 eV region is followed by a broad diffuse band whose maximum is at ∼4.76 eV but which contains small features assigned to the 31A′ ← 11A′ π–π* transition, in particular band number 47 at 4.564 eV as the O00 origin band. Our DFT calculation places this transition at 4.86 eV, f = 0.034 (Table 3). Several P–P–P studies calculate this transition to be in the 4.4–4.6 eV region (Table 5). The two lowest singlet π–π* transitions, respectively 21A′ ← 11A′ and 31A′ ← 11A′, were listed by Baba and Yamazaki[44] to have been observed at 3.91 eV, f = 0.023 (calc 0.118) and 4.66 eV, f = 0.085 f (calc) = 0.147 respectively. We remark that Wagner et al.,[48] who calculate the 31A′ ← 11A′ transition at 4.83 eV (f = 0.24), consider that the 31A′ electronic state corresponds to 1La, using the designation of Platt which is based on a free electron model,[60] whereas Zimmerman and Joop[56] made an 1Lb (long axis polarised) assignment as indicated by their measurements of fluorescence polarization. We note also that Favini et al.[24] studied the absorption spectrum of isoquinoline up to about 5.4 eV in cyclohexane and methanol solutions, i.e. in the range 220–400 nm. The 270 nm (4.59 eV) broad band maximum shifts to higher energies in cyclohexane (268 nm, 4.63 eV in methanol and to 265 nm, 4.68 eV in cyclohexane), and the spectrum exhibits more structure in cyclohexane.

Valence transitions: 5.2–7.2 eV

This is also a structured region (Fig. 5–7). Bands numbered 58–60 are assigned to the 41A′ ← 11A′ and 51A′ ← 11A′ π–π* transitions. The oscillator strengths for these transitions differ strongly in their calculated values in Table 5, reflecting differences in the complexity of the electron correlations in the various theoretical methods. The absorption cross-section then rises to a single sharp band at 5.864 eV (band no. 61) and then to a maximum at 6.003 eV (band no. 63). The bands number 61 and 62 are Rydberg features, as described later (Section 3.1.4.1). The absorption cross section continues as a decreasingly intense broad band up to ∼6.8 eV but contains three shoulder features (i) at 6.038 eV (band no. 64) assigned to the π–π* 61A′ ← 11A′ O00 transition, (ii) at 6.139 eV (band no. 66), assigned to a vibronic component of the latter and (iii) at 6.322 eV (band no. 69) assigned to two overlapping transitions, the π–π* 71A′ ← 11A′ O00 band, and the R3 Rydberg n = 3 O00 transition whose Rydberg limit is the first excited ion state at 9.16 eV (see below). Two more π–π* transitions are assigned in the 6.8–7.1 eV region: 81A′ ← 11A′ and 91A′ ← 11A′ and their vibronic components (bands numbered 72–74 and 76). The peak at 6.003 eV has an absorption cross section of 292 Mb. Its high intensity can be considered to be due to plasmon-type effects that result from the existence of collective modes of the π electrons.[61-65] These many-body effects are not considered in Hartree–Fock method single particle calculations. We note that naphthalene, which is closely related to isoquinoline, has a strong absorption peak at 5.89 eV (ref. 66) that is considered as being of plasmon-type,[65] resulting in an increased intensity of the 1Bb π–π* valence transition.[67]

Tests of valence transition energy calculations

The calculated and observed transition energies of isoquinoline are compared in Table 5. The tabulated values of UD, the unsigned difference between the DFT calculated energies and the eight assigned transition energies extend from 0.1058 to 0.5847 eV, with an average value of 0.3305 eV. There have been many time-dependent DFT benchmark test calculations of the electronic state energies of both diatomic and polyatomic molecules.[28,68-78] These calculations are principally of singlet valence excited states although there are a few that also include triplet states and Rydberg electronic states. These benchmark calculations provide comparisons of the efficiencies of various functionals, one of the best being among the hybrid density functional B3LYP, which contains 20% of Hartree–Fock exchange.[78] A typical criterion of the quality of TD-DFT calculation performance has been defined: when, in TD-DFT calculations of excited electronic states, the mean unsigned differences (MUD) between experimental and calculated energies are less than 0.36 eV the TD-DFT calculations are considered to be of value,[78] although it should be stated that in a typical set of benchmark molecules the actual UD values often have a large range, of up to or greater than 1 eV. Our TD-TFT calculation results on isoquinoline are thus consistent with the MUD < 0.36 eV criterion. There are other useful criteria to be considered, dependent on the degree of precision required. Examination of Table 5 shows that our calculated energies are far from so-called spectroscopic accuracy, and from the much less rigorous chemical accuracy of 0.05 eV. The calculated energies can be considered as “coarse-grained values”, yet are useful in spectral assignments. It is of interest to consider what are the possible sources of errors in TD-DFT calculations. These are discussed by Burke et al.,[43] who mention the problem of adequately taking the exchange–correlation potential into account. In addition there is the possibility that the local adiabatic approximation no longer holds. This approximation is of importance in the case of frequency independent exchange–correlation functionals with which the exchange–correlation potential adjusts instantaneously to variations of the electron density. Prilj et al.[79] have discussed the origins of errors in TD-DFT calculations of the excited state energies of heteroaromatic molecules. They stress that Lb states are sensitive to correlation effects, and that La states have considerable double excitation characteristics. Table 5 also lists the UD values cerning the π–π* excitation energies calculated by Baba and Yamazaki[44] by the P–P–P SCF LCAO method including configuration interaction. These are the lowest set of P–P–P calculation UD values of the studies listed in Table 5. The UD values of Baba and Yamazaki range from 0.049 to 0.567 eV, with an average value of 0.232 eV, which is 30% better than the average UD in our DFT calculation results. We note that the UD values decrease with increasing excitation energy in our DFT calculation, but increase in the P–P–P calculations of Baba and Yamazaki.

Rydberg and other transitions: 7.2–10.8 eV

The 7.2–10.8 eV spectral region exhibits strong π–π* transition continua underlying a number of relatively weak features many of which we assigned to Rydberg bands. There are no previous observations of Rydberg bands in isoquinoline spectra. The bands which we assign to Rydberg series are found to converge to the ground and also to excited electronic states of the ion, as discussed below. The broad diffuse band in the 7.5–8.5 eV region has a maximum at 7.88 eV, where the absorption cross section is 75 Mb (Fig. 3 and 7). In the region from 7.7 eV to the maximum at 7.88 eV there is a unique feature, between bands number 84 and 85 (Fig. 3 and 7) exhibiting a sharp dip in intensity to a minimum absorption cross-section of 56.7 Mb at 7.754 eV. The intensity profile in this spectral region strongly resembles some absorption profiles in the 7.5 eV region of the naphthalene spectrum[61,66,80-82] that have been interpreted as anti-resonance features resulting from interaction between Rydberg states and isoenergetic valence states that have a large density of vibronic states inhomogenously broadened by intramolecular vibrational redistribution.[81,82] The situation is that of the statistical limit case of radiationless transitions[83] in which here a Rydberg state interacts with an effective continuum formed by the broadened valence state. We suggest that the 7.7–7.85 eV absorption profile in isoquinoline is another example of such interference effects. We do not attempt to precisely identify the Rydberg and valence states concerned in this sharp antiresonance process but it is interesting to note that three close-lying singlet π–π* valence states, 101A′, 111A′ and 121A′, are calculated by our DFT calculation (Tables 3 and 5) to occur in the 7.5–7.7 eV region where several Rydberg bands are observed (Table 4). Beyond 8.5 eV and up to the 10.8 eV limit, where the absorption cross section is also ∼77 Mb, there is a series of broad diffuse structures with superposed very small features, many of which have been assigned as Rydberg transition features. It is noteworthy that the absorption markedly increases at about 8.53 eV, 9.16 eV and 10.32 eV, corresponding to successive ionization energies of isoquinoline.[84,85] These inflexions of the absorption curve are very prominent and are discussed in more detail below. The origin of the background continua in isoquinoline requires further theoretical and experimental investigation, in particular concerning the possible existence, at the high energy end of the measured absorption spectrum, of transitions other than π–π* and n–π* transitions, since π–σ*, σ–π* and σ–σ* transitions could also exist in this spectral region.[62] The latter are certainly among those transitions given by our DFT calculations of electronic state energies and transition oscillator strengths in the 8–11.4 eV region (Table 3), as can be seen in the various molecular orbital transition components associated with these high energy transitions (not reported).

Rydberg series converging to the ion ground electronic state at 8.53 eV

We identified three Rydberg series that converge on the D0 doublet ground state of the isoquinoline cation. For our analysis we use the equation E(n) = IE − R/(n − δ)2 for the energy E(n) of the nth Rydberg series member, where IE (isoquinoline) = 8.53 eV, the Rydberg constant R = 13.606 eV, and δ is the value of the quantum defect. The n = 3 first member of the R1 Rydberg series (Table 4, Fig. 3) is assigned to band number 61 (R1, n = 3, δ = 0.741). This band is observed in a spectral region notable for π–π* valence transitions. The n = 4, 5, 6, … (Table 4) members have greater values (δ = 1.0 ± 0.1) of their quantum defects. We can rationalize these higher δ values with respect to n = 3 by noting that the radius of the Rydberg orbital in the n = 3 Rydberg state is much the same as that of the ion core. This can result in coupling between Rydberg and π orbitals which would create departure from the hydrogen ion core based model formula used above to determine E(n) values. Naphthalene, a parent molecule of isoquinoline, exhibits an analogous difference between n = 3 and n ≥ 4 Rydberg level quantum defect values in its lowest Rydberg series where the n = 3 members have δ ∼0.68,[86] and n ≥ 4 Rydberg levels δ ∼0.85.[66] Rydberg bands of the R1 series have been observed up to n = 12, the latter being near ionization limit of 8.53 eV. Band resolution limitations are considered to be responsible for the range of quantum defect values around δ ∼0.88 found for the greatest values of n in this Rydberg series. We were unable to assign the n = 10 member, expected to occur at about 8.366 eV for δ = 0.9. This Rydberg level may have been displaced by interaction with a valence state such as the 181A′ valence state calculated to occur at 8.3833 eV (Table 3). The existence of a smaller δ values for the smaller n members of a Rydberg series can be due to a more marked effect of the core potential in comparison to the situation existing in higher n levels. The quantum defect values of the R1 Rydberg series bands lead us to assign R1 as an ns series in which an electron is ejected from the HOMO π−1 orbital (Table 7).

Isoquinoline Rydberg series

Rydberg seriesIon state energy limit of seriesObserved membersAverage quantum defect δRydberg electronic transition
R1D0 8.53 eV n = 3–12 δ ∼ 0.9011A′ → (HOMO)π−1ns
R2D0 8.53 eV n = 3–8 δ ∼ 0.6811A′ → (HOMO)π−1np1
R3D0 8.53 eV n = 3–6 δ ∼ 0.2811A′ → (HOMO)π−1np2
R4D2 9.16 eV n = 3–7 δ ∼ 0.8511A′ → (HOMO−2)π−1ns
R5D2 9.16 eV n = 3–6 δ ∼ 0.6211A′ → (HOMO−2)π−1np1
R6D3 10.40 eV n = 3–8 δ ∼ 0.9011A′ → (HOMO−3)π−1ns
R7D3 10.40 eV n = 3–9 δ ∼ 0.6811A′ → (HOMO−3)π−1np1
R8D3 10.40 eV n = 3–6 δ ∼ 0.2711A′ → (HOMO−3)π−1np2
Some of the R1 series bands correspond to vibrational components of the Rydberg transitions (Table 4). The vibrational frequencies are ν(x) ∼1160 cm −1 and ν(z) ∼1320 cm −1. Galué et al. have measured an IRMPD spectrum of the ground state of the isoquinoline cation.[87] This spectrum exhibits a broad strong band at ∼750 cm−1, two broad bands at about 1000 and 1100 cm−1 respectively, and strong, but broad, bands at ∼1200 and ∼1450 cm−1. The latter vibrational frequencies, close to those of the R1 Rydberg series vibrational components, probably correspond to carbon–carbon stretch modes. A second group of Rydberg bands that converge to the D0 ground state at 8.53 eV was assigned for the n = 3–8 levels, forming the R2 series. Their quantum defect values range from δ = 0.61 to 0.74 (Table 4, Fig. 3). The R2 series bands do not have any recognizable vibrational features. We consider this Rydberg series, assigned as an np1 series (Table 7), to also be associated with the electron ejection from the HOMO π−1 orbital. Another series of Rydberg bands, R3, whose six members, n = 3–8, converge to the ion ground state, with a quantum defects ∼0.28, is assigned as 11A′ → (HOMO)π−1np2. The absorption spectrum in the 8 eV region exhibits a strong background π–π* transition continuum (Fig. 3 and 7). As mentioned above, the absorption intensity in this region reaches a maximum at 7.88 eV (σabs = 75 Mb), and is followed by a minimum at about 8.47 eV, close to IE(isoquinoline) = 8.53 eV. At higher energies there is a broad continuum whose intensity increases to a maximum at 8.64 eV (σabs = 38.8 Mb), drops to a minimum at 8.99 eV, and then rises to a maximum at 9.52 eV (σabs = 55 Mb). The minima correspond in energy to the onset of ionization channels, 8.47 eV reflecting the onset of the first ionization of isoquinoline, and 8.99 eV the onset of higher energy ionization channels mentioned previously. The cross section decreases to a minimum at 10.23 eV, just below the 10.40 eV ionization limit after which it rises to plateau beginning at about 10.4 eV where σabs ∼70 Mb.

Rydberg series converging to the ion excited electronic state at 9.16 eV

The HeI photoelectron spectrum of isoquinoline has been interpreted by Brogli et al.[85] as having two close-lying ion excited states, of π−1 and n−1 character, in the 9.16 eV region. In their diagram correlating molecular orbitals of naphthalene and isoquinoline they effectively assign the n−1 state to be D1 and the π−1 state to be D2, a classification that we follow. We assigned five bands to the R4 Rydberg series converging on the D2 ion excited state, considered as π−1, at 9.16 eV. This limit corresponds to a single clear onset in the absorption cross-section (Fig. 3, 7 and 8) and is also the ionization energy reported by Eland and Danby from their photoelectron spectroscopy study.[84] Brogli et al.[85] gave an uncertain value of 9.30 eV for this ionization energy from their HeI PES spectrum. We tried out both 9.16 and 9.3 eV as possible Rydberg limit values and found that 9.16 eV enabled us to assign five R4 bands, with consistent quantum defects whereas a R4 limit of 9.3 eV was far less satisfactory. The definitive identification of two separate π−1 and n−1 ion states in this energy region requires new high resolution optical and photoelectron spectroscopy studies. The five assigned R4 features are bands number 69 ((n = 3) O00δ = 0.81), 86 ((n = 4) O00, δ = 0.82), 100 ((n = 5) O00δ = 0.90), 106 ((n = 6) O00, δ = 0.89) and 108 ((n = 7) O00, δ = 1.08) (Table 4). The Rydberg transition is 11A′ →(HOMO−2)π−1ns. A second Rydberg series of four members converging to the D2 ion level was assigned: R5, n = 3–6, δ ∼0.62. The corresponding Rydberg transition is 11A′ →(HOMO−2)π−1 np1. Attempts to find corresponding Rydberg series converging to the other electronic excited state of the ion in this energy region, considered as the 11A′ →(HOMO −2)n−1 D1 electronic state, were not successful.

Rydberg series converging to the ion excited electronic state at 10.40 eV

Six bands are assigned to the R6 Rydberg series converging to the ion D3 excited state at 10.40 eV. Their principal quantum numbers are n = 3–8, with quantum defects of the order δ = 0.9 ± 0.1 (Table 4). Band number 145 at 10.204 eV may represent an unresolved overlap of the n = 9 and n = 10 Rydberg bands and band number 146 unresolved further members of this Rydberg series closer to the limit at 10.40 eV. It is interesting that the quantum defects of the n = 3 terms of the various ns Rydberg series increase with the Rydberg limit energy: δ = 0.74 (R1), δ = 0.81 (R4), δ = 0.85 (R6). This indicates that the n = 3 electron penetration to the nucleus is more effectively shielded as the ionization limit increases. Two other Rydberg series converging to the D3 excited state were assigned: R7, with seven components, n = 3–9, δ ∼0.68, considered as a 11A′ →(HOMO−3)π−1np1 series and the more limited R8 series whose four components, n = 3–6, have quantum defects δ ∼0.27 corresponding to a np2 series. We remark that the D3 ← D0 electronic transition spectrum of the cation of isoquinoline, observed by photodissociation spectroscopy,[88] exhibits vibrational components, in particular with frequencies assigned as ν30 = 511 cm−1, ν27 = 795 cm−1 and ν13 = 1412 cm−1. It is often the case that Rydberg vibrational components are similar to those observed in transitions to electronic states of a molecular ion achieved by direct photoionization of a neutral molecule. The vibrational modes excited can, however, be different than those excited in an electronic transition within a molecular ion, since the Franck–Condon factors are not necessarily the same in the two cases. Examination of the Rydberg series converging to the D3 excited state of the isoquinoline cation failed to observed bands that would correspond to vibrational components of the D3 ← D0 cation transition seen in the photodissociation spectra of Dryza.[88] We then compared the cation spectra of isoquinoline[88] with bands occurring in direct photoionization to the D3 excited states of the cation observed in HeI photoelectron spectra.[84,85] This showed that there is apparently a considerable difference in Franck–Condon factors in the two cases. We also examined the analogous spectra in naphthalene and came to the same conclusion as to a difference in Franck–Condon factors in these cases.[66,84,85,88,89]

Comparison between quinoline and isoquinoline absorption spectra

It is of interest to compare the absorption spectra of quinoline[13] and isoquinoline explored over the same energy region 3.5–10.7 eV. There are many similarities but also some interesting differences between the two spectra. Both C9H7N molecules have Cs symmetry and the same number of a′ and a′′ vibrational modes whose corresponding frequencies are of the same order of magnitude but with some differences according to whether the N atom is in the α (quinoline) or the β (isoquinoline) position. Quinoline and isoquinoline are isoelectronic with naphthalene, whose valence transitions are all π–π* but the azamolecules possess non-bonding electrons associated with the nitrogen atom that give rise, in addition, to n–π* transitions. A comparison between the calculated structures of quinoline (Fig. 1 in ref. 13) and isoquinoline (Fig. 1) indicates that changing the nitrogen atom from the α to the β position leads to substantial changes in bond lengths and ring angles in the ring containing the N atom but little modification of the other aromatic ring. Molecular orbital correlations between naphthalene, quinoline and isoquinoline can be used to rationalise a number of spectral properties, here the valence and the Rydberg transitions of the azamolecules. We note that eight valence transitions m1A′ ← 11A′ were observed in quinoline:[13] the transitions m = 2, 3,…8, whose origin bands lie between 3.991 and 6.981 eV, and the m = 12 transition at 7.795 eV. Isoquinoline also has eight assigned m1A′ ← 11A′ valence transitions: m = 2, 3, …9, between 3.960 and 7.081 eV. The energy of the 21A′ ← 11A′ transition in isoquinoline is 31 meV smaller than the corresponding transition in quinoline.[13] On the other hand, the 31A′ ← 11A′ transition in isoquinoline at 4.564 eV is at a higher energy than in quinoline (4.441 eV). These differences can be rationalised on a Koopmans theorem[90] approximation approach, using ionization energies to assess molecular orbital energy differences and considering the optical transitions on a single configuration basis. From the known ionization energies of quinoline and isoquinoline we estimate that the energy of the 21A′ ← 11A′ transition in isoquinoline would lie 98 meV below that in quinoline as compared with the experimental value 31 meV, in agreement with the sign of the difference. A better quantitative agreement is not to be expected since this approach does not take into account effects of interaction between close-lying 21A′ππ* and 11A′′nπ* states, whose zero-order differences in energy differ significantly in quinoline and isoquinoline as discussed below. With this approach we also calculate the energy difference between the 31A′ ← 11A′ and the 21A′ ← 11A′ transitions to be 552 meV in quinoline and 630 meV in isoquinoline; our experimental values, respectively 450 meV and 604 meV, are in good agreement with the calculated values. We remark that in this single configuration Koopmans theorem approach to molecular orbital energies we also neglect the effects of rapid orbital reorganisation that occurs on photoionization.[91] Our DFT and the Baba and Yamazaki[44] P–P–P calculations of π–π* transition energies of both quinoline[13] and isoquinoline gave similar levels of accuracy with respect to observed values, as measured by UD, the unsigned energy difference, the average DFT and P–P–P values being respectively UD = 260 and 340 meV in quinoline,[13] and UD = 331 and 232 meV in isoquinoline. The UD values are of the same order of magnitude but the DFT calculations are overall more accurate than the P–P–P calculations in quinoline, the reverse being found in isoquinoline. The 3.8–4.4 eV spectral region exhibits vibronic interaction effects between the 21A′ππ* and 11A′′nπ* electronic states in both quinoline and isoquinoline. The energy gap between the nπ* and ππ* origin levels is greater in quinoline (1804 cm−1)[16] than in isoquinoline (1118 cm−1). This gives rise to very different vibronic structures in this energy region for the two azamolecules. These structures have been well analysed in isoquinoline (Tables 4 and 6) but more scantily in quinoline where the features are more complex.[16] Both quinoline and isoquinoline have strong broad absorption peaks in the 6 eV region, with maxima at 6.218 eV, σ = 194 Mb in quinoline[7] and 6.003 eV, σ = 292 Mb in isoquinoline. The high intensity is considered to be due in part to plasmon-type effects as discussed above. The larger peak cross-section in isoquinoline indicates a greater influence of the collective modes of the π electrons in this energy region. The difference in the peak energies and intensities marks the difference between the π electron structures in the α-N and β-N azamolecules. Previous to our work on the VUV absorption spectra of quinoline[13] and isoquinoline no Rydberg transitions had been reported for these two molecules. We observed eight Rydberg series in quinoline and eight in isoquinoline. In quinoline there are two series, ns and np1, that converge to the D0 level of the cation (HOMO)π−1, three series, ns, np1 and np2, to the D3 state (HOMO−3)π−1 and three series to D4 (HOMO−4)π−1. Rydberg series were not observed converging to the D1 state (HOMO−1)π−1 or to the D2 level (HOMO−2)n−1. Rydberg series observed in isoquinoline (Table 7) converge to the D0 (HOMO)π−1, D2(HOMO−2)π−1 and D3 (HOMO−3)π−1 electronic states of the ion. A series converging to the D1 level (HOMO−2)n−1 was not observed. That there are differences between the Rydberg series observed for the two molecules can be considered as resulting, in part, from differences in the order of the π and the σ-type n orbitals as expressed in the M.O. correlation diagram.[85] We mention also that an antiresonance feature due to interaction between a Rydberg level and a quasi-continuum of π−π* transition character was observed in isoquinoline at about 7.75 eV; no similar feature was observed in the absorption spectrum of quinoline. Further experimental and theoretical studies on these rich absorption spectra are required to advance our understanding of the valence and Rydberg transitions in these azamolecules.

Conclusion

The gas-phase absorption spectrum of isoquinoline was measured over the photon energy range 3.6–10.7 eV using the ASTRID2 synchrotron radiation facility at Aarhus, Denmark as the photon source, as was previously done over the same energy range for quinoline.[13] Previous absorption spectral studies on isoquinoline were mainly in solution, with an upper energy limit of 5.6 eV. Quantum-mechanical calculations, using time-dependent DFT methods, were carried out to determine properties of the ground state of isoquinoline such as geometry, rotational constants, vibrational frequencies and dipole moment. The results are in very satisfactory agreement with experimental data. Previously published results of Pariser–Parr–Pople (P–P–P) quantum mechanical calculations of electronic state energy levels and dipole transition strengths, as well as the results of our DFT calculations, enabled us to assign eight π–π* singlet m1A′ ← 11A′ valence electronic transitions of isoquinoline. This considerably increases the number of electronic transitions previously known for isoquinoline. We discuss the relative quality and trends of the results of the P–P–P and DFT calculations. Although they are far from spectroscopic or chemical accuracy, these results have been of much use in assigning the valence transitions. No n–π* transitions of isoquinoline were observed, coherent with their low DFT- calculated oscillator strengths. Some of the broad spectral features, such as an intense broad band peaking at 6.003 eV, are suggested to have intensity contributions from plasmon-type collective π electron mode effects. A feature in the 7.7–7.8 eV region has an absorption profile reminiscent of an anti-resonance between Rydberg and valence levels. Rydberg bands of isoquinoline were never previously reported. We observed and assigned eight Rydberg series, three of which, R1, R2 and R3 converge to the D0 ion ground state. Two further Rydberg series, denoted respectively as R4 and R5, converge to the D2 excited state of the cation at 9.16 eV, and three more, the R6, R7 and R8 series, converge to the D3 excited state at 10.40 eV. There appears to be a trend to a lesser penetrability of the 3 s electron as the ionization limit increases with successive ion electronic states, reflected in the increase of the 3 s electron quantum defect. There are 161 measured features in our absorption spectrum, of which we were unable to assign 80, mainly high-energy features. This is a region where we expect to find high energy valence transitions, as well as Rydberg bands converging to highly excited cation states, extending to or above the D3 level energy.[85] In the lower energy region, vibronic interactions between n–π* and π–π* electronic states can cause apparent irregularities in band energies that require very high resolution spectroscopy and associated photophysical measurements to untangle, as has been done for the origin region presented in Table 6. Interesting new studies are thus called for. The isoquinoline absorption cross-sections that we measured over an extended energy region that includes the VUV down to 10.7 eV represent new useful data for astrophysics, in particular for theoretical modeling of the photochemistry of isoquinoline in planetary atmospheres, the ISM and circumstellar media, as mentioned in the introduction where it is noted that quinoline and isoquinoline could be synthesised under high temperature conditions representing circumstellar envelopes of carbon stars.[7,8] Radiofrequency spectroscopic attempts to observe both C9H7N isomers quinoline and isoquinoline, in the circumstellar envelopes of carbon-rich stars have not been successful beyond providing upper limits on quinoline and isoquinoline column densities.[9] Further modelling exploitation of the photon absorption spectra of quinoline and isoquinoline in the VUV will require information on photoionization quantum yields as function of excitation energy, in particular in energy regions involving superexcited electronic states of neutral molecules.[12] In the case of isoquinoline we have used a rule-of-thumb method[10] to evaluate the photoionization quantum yield, γi (13.6 eV) = 0.55, at the 13.6 eV energy limit of the astrophysical HI region. Thus, at 13.6 eV, as in the analogous case of quinoline,[13] about 45% of the absorbed photon energy is available for relaxation processes of the neutral molecule. Such processes include molecular dissociation. At 10.7 eV, the energy upper limit of our absorption spectrum, the rule-of-thumb method indicates that about 75% of the incident photon energy is dissipated in non-ionization photophysical processes. A comparison of the absorption spectra of quinoline and isoquinoline over the same energy range is given in Section 4 in which various spectral aspects of these species are considered as being due, in part, to specific characteristics of the respective molecular orbitals in the two azamolecules. The many similarities, and some interesting differences, in the absorption spectra of these two isomers, can have a bearing on their respective photophysical properties and eventually on their role in astrophysical processes.

Conflicts of interest

There are no conflicts to declare.
  23 in total

1.  Accurate Rydberg excitations from the local density approximation.

Authors:  Adam Wasserman; Neepa T Maitra; Kieron Burke
Journal:  Phys Rev Lett       Date:  2003-12-29       Impact factor: 9.161

2.  Benchmarking the performance of time-dependent density functional methods.

Authors:  Sarom S Leang; Federico Zahariev; Mark S Gordon
Journal:  J Chem Phys       Date:  2012-03-14       Impact factor: 3.488

3.  Electronic Transition Energies: A Study of the Performance of a Large Range of Single Reference Density Functional and Wave Function Methods on Valence and Rydberg States Compared to Experiment.

Authors:  Marco Caricato; Gary W Trucks; Michael J Frisch; Kenneth B Wiberg
Journal:  J Chem Theory Comput       Date:  2010-01-22       Impact factor: 6.006

4.  Interpretation of the vacuum ultraviolet photoabsorption spectrum of iodobenzene by ab initio computations.

Authors:  Michael H Palmer; Trevor Ridley; Søren Vrønning Hoffmann; Nykola C Jones; Marcello Coreno; Monica de Simone; Cesare Grazioli; Malgorzata Biczysko; Alberto Baiardi; Paulo Limão-Vieira
Journal:  J Chem Phys       Date:  2015-04-07       Impact factor: 3.488

5.  Anharmonic vibrational properties by a fully automated second-order perturbative approach.

Authors:  Vincenzo Barone
Journal:  J Chem Phys       Date:  2005-01-01       Impact factor: 3.488

6.  Time-dependent density functional theory: past, present, and future.

Authors:  Kieron Burke; Jan Werschnik; E K U Gross
Journal:  J Chem Phys       Date:  2005-08-08       Impact factor: 3.488

7.  Computation of accurate excitation energies for large organic molecules with double-hybrid density functionals.

Authors:  Lars Goerigk; Jonas Moellmann; Stefan Grimme
Journal:  Phys Chem Chem Phys       Date:  2009-03-23       Impact factor: 3.676

8.  High-Resolution Far-Infrared Spectroscopy of N-Substituted Two-Ring Polycyclic Aromatic Hydrocarbons: An Extended Study.

Authors:  S Gruet; O Pirali; M Goubet; D W Tokaryk; P Brechignac
Journal:  J Phys Chem A       Date:  2015-12-30       Impact factor: 2.781

9.  Electronic spectra of gas-phase polycyclic aromatic nitrogen heterocycle cations: isoquinoline+ and quinoline+.

Authors:  Viktoras Dryza; Julian A Sanelli; Evan G Robertson; Evan J Bieske
Journal:  J Phys Chem A       Date:  2012-04-19       Impact factor: 2.781

10.  Low-Lying ππ* States of Heteroaromatic Molecules: A Challenge for Excited State Methods.

Authors:  Antonio Prlj; María Eugenia Sandoval-Salinas; David Casanova; Denis Jacquemin; Clémence Corminboeuf
Journal:  J Chem Theory Comput       Date:  2016-05-11       Impact factor: 6.006

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.