extfieldsofstudy:
Mathematics
Computer Science
Medicine

text:
A hybrid least squares and principal component analysis algorithm for Raman spectroscopy

The least squares fitting algorithm is the most commonly used algorithm in Raman spectroscopy. In this paper, however, we show that it is sensitive to variations in the background signal when the signal of interest is weak. To address this problem, we propose a novel algorithm to analyze measured spectra in Raman spectroscopy. The method is a hybrid least squares and principal component analysis algorithm. It explicitly accounts for any variations expected in the reference spectra used in the signal decomposition. We compare the novel algorithm to the least squares method with a low-order polynomial residual model, and demonstrate the novel algorithm's superior performance by comparing quantitative error metrics. Our experiments use both simulated data and data acquired from an in vitro solution of Raman-enhanced gold nanoparticles.

Introduction
Background Raman spectroscopy is a powerful technique for analyzing chemical compounds using a laser source. It exploits the Raman effect, which arises from the interaction between laser light and a sample of interest. When incoming photons hit the sample surface, most photons are scattered elastically, after which they continue traveling with the same energy and wavelength. However, a very small fraction of the photons is scattered inelastically, meaning that they lose energy and continue traveling with a longer wavelength. The amount of energy lost by the photons depends on the particular molecules they interact with. In fact, the chemical bonds in the molecules absorb energy in highly specific patterns. As a result, the Raman scattered photons possess highly compoundspecific wavelength spectra. Raman spectroscopy uses these highly specific spectral fingerprints to identify and quantify compound concentrations. The usefulness and power of Raman spectroscopy lie in the fact that it allows rapid sample analysis of single or multiple compounds (known as multiplexed analysis) at high detection sensitivities [1]. Amongst its many promising areas of application, Raman spectroscopy has gained growing interest from the biomedical research community, where it promises to enable sensitive imaging of nanoparticles for both diagnostic and therapeutic applications [1][2][3]. Examples of such applications are Raman colonoscopy for early cancer detection and improved tumor margin detection during surgery.
A critical part of Raman spectroscopy is the use of an appropriate signal analysis algorithm to analyze the measured Raman spectra. This paper focuses on the development of a signal analysis algorithm that is robust to natural variations in both the background and analyte signal. In the following sections, we describe the previous literature on spectral analysis algorithms, highlight their strengths and weaknesses, and explain the need for and novelty of our contribution.

Previous Raman Spectral Analysis Algorithms
Various methods have been used to analyze Raman spectra (i.e. detect and quantify compounds of interest), such as classical least squares [1], least squares with a low-order polynomial background model [4], variable baseline correction [5][6][7], explicit detection and parametric (Gaussian) modeling of Raman peaks [8], principal component regression [9], partial least squares [10,11], and hybrid linear analysis [12]. The classical least squares method can be used when the pure spectra of the compounds of interest (also known as analytes) and an accurate background spectrum are known. While fast and quantitative, this method is sensitive to inaccuracies or variations in the spectra of the compounds of interest and the background. It is also sensitive to noise when the signal of interest is weak. Lutz et. al [4] addressed the problem of background variability by allowing the background spectrum to vary according to a low-order polynomial model. While successful in accounting for slowly varying changes to the background spectrum, it cannot accommodate higher-order variations such as slight peak shifts or changes in the relative amplitudes of peaks, all of which are regularly observed in practice. Many alternative baseline fitting algorithms suffer from the same limitation [5][6][7]. Such methods commonly subtract a smoothly varying baseline estimate from the measured spectra before computing concentration estimates, often using a conventional least squares approach. Again, while suitable to account for smooth background variations such as those caused by autofluorescence, these models are not adequate in the presence of complex sources of variation such as changes in peak position and relative peak amplitudes. Aiming to model such shift and amplitude changes directly, Kode et al. [8] proposed to model spectrum peaks explicitly with 1D Gaussians. By employing a penalized cost function, the method allows the peaks to modestly change position, width and amplitude while computing concentration estimates. However, such peak detection and quantification is sensitive to noise at low signal strengths.
Principal component regression (PCR) and partial least squares (PLS) are implicit methods that require neither prior knowledge of the reference spectra, nor an explicit background model. As such, they can handle background variations that are more complex than the smoothly varying curves discussed above. Known as multivariate calibration methods, they attempt to find a linear model that relates a dependent variable, e.g. analyte concentration, to the measured independent variables, e.g. spectra, for complex mixtures. The model parameters are obtained in the form of a regression or calibration vector b, which is derived from a calibration set of representative mixtures for which the dependent variable is known. Subsequently, the dependent variable (analyte concentration) is predicted by taking the dot product of b with the measured spectrum of an unknown mixture.
While PCR and PLS are broadly applicable, they also ignore valuable knowledge of the analyte spectra when available. This scenario was encountered previously by Berger et al. [12], who possessed a mixture calibration set with variable spectra, knowledge of the analyte concentration, as well as the pure analyte spectrum. In a bid to improve on PCR and PLS, Berger et al. proposed to exploit the additional information using a method they called hybrid linear analysis (HLA). This technique represents the variable background signal as a linear combination of the background signals' principal components. It obtains an accurate calibration vector by estimating the background signals (i.e. the calibration signals without the analyte contribution), and subtracting from the known analyte spectrum its projections onto each of the background signals' principal components (see the section entitled 'Hybrid Linear Analysis (HLA) Method'). This technique can be repeated to derive the calibration vector for any known analyte in the mixture. HLA was shown to significantly outperform PLS (PCR was not tested because its performance is usually similar but slightly inferior to that of PLS). This result makes intuitive sense because HLA uses more physical information to obtain the calibration vector than PLS (or PCR).
In the study reported here, we possess a calibration set of variable background signals excluding the analyte, as well as a calibration set of pure analyte spectra. The analyte consists of Raman-enhanced gold-silica nanoparticles (see the Simulation Results section). Since the background and analyte variations cannot be modeled by smooth curves (see the Results and Discussion section), multivariate calibration techniques are better suited than baseline correction methods. Of these techniques, HLA is more suitable than PCR or PLS due to the availability of pure analyte spectrum information. HLA can be used by skipping the first step of background isolation (since we measure it directly), and working with the mean of the analyte spectrum. However, such an approach ignores the information about the analyte spectrum variation. It also fails to incorporate information about the extent of variation observed in the calibration sets, as captured by the eigenvalues of the principal component analyses.

Hybrid Least Squares and Principal Component Analysis Algorithm
In this study, we propose to model variations in both the background and analyte spectra, in order to increase the robustness of analyte concentration estimates. Improving on HLA, our method incorporates the principal components as well as eigenvalues of the background and analyte calibration sets into a hybrid least squares and principal component analysis (HLP) method. Our method differs from standard multivariate calibration techniques in that it does not derive a calibration vector. Instead, HLP estimates fitting weights for each of the analyte and the background signals. We explain the mathematical details of the method in the section entitled 'Novel Hybrid Algorithm (HLP)'. HLP is tested on both simulated data and experimental data acquired from an in vitro solution of Raman-enhanced gold-silica nanoparticles [1]. The results are presented in the Results and Discussion section, where we demonstrate the improved performance of the novel method compared to that of the least squares method with a low-order polynomial background model [4] and HLA. To reiterate, we chose HLA as a competitive method for comparison, because it is robust to complex variations in the background signal, and incorporates available information about the analyte spectrum. The latter property makes it more suitable for comparison to our HLP algorithm than PCR or PLS. Lastly, the Conclusions and Further Work section draws conclusions and describes avenues for future work.

Classical Least Squares Method
Here we briefly review the classical least squares method with an added low-order polynomial background model, as proposed by Lutz. et al. [4]. The measured spectrum can be modeled as a linear combination of known spectra (a.k.a. reference spectra): where x l is the modeled intensity at wavelength l, K is the number of reference spectra provided, S lk is the value of the reference spectrum of the k th compound at wavelength l, and w k is the weight for the k th compound. In the method proposed by Lutz et al. [4], the spectra S k include the compounds of interest as well as an average background signal and the q+1 components of a q th order polynomial. The concentrations of the various compounds are then estimated by solving for the weights w k that give the closest fit with the measured spectrum. This can be done by writing the problem in the matrix form M = SW, where M is the L61 vector containing the measured spectrum values m l , S is the L6K matrix of reference spectra, and W is the K61 matrix containing the weights w k . The least squares solution is then given bŷ where S {~ST S ð Þ {1 S T is the Moore-Penrose pseudoinverse of the matrix S. Note that this algorithm implicitly assumes that the noise on the signals is Gaussian distributed.

Hybrid Linear Analysis (HLA) Method
As mentioned in the section on previous Raman spectral analysis algorithms, the HLA method proposed by Berger et al. [12] estimates the concentration of an analyte by taking the dot product between a calibration vector b and a measured spectrum m. The calibration vector b is computed from a calibration set of mixture spectra, for each of which the analyte concentration is known, as well as an accurate estimate of the pure analyte spectrum. Let B A denote the calibration set, where the sample spectra are stored in the rows. The corresponding concentrations are stored in a column vector k. The known spectrum of the analyte is represented by a row vector S A , measured at unit concentration. Using B A , k and S A , the calibration vector b is computed as follows: 1. Isolate the background signals by subtracting out the estimated spectral contributions of the analyte: 2. Compute the principal components of the background calibration set B, and store them in the rows of a matrix V. Since all of the spectra in B can be modeled by those in V, the spectra in V act like pure spectra of the background species.
1. Subtract from the reference spectrum S A its projections onto each of the background's principal components, leaving a residual spectrum r: This residual r is the portion of S A that cannot be modeled by the spectra contained in V, i.e. it is orthogonal to V.
2. Normalize r to get the calibration vector b: where '.' denotes the dot product.
This procedure can be repeated for any analyte of interest. As mentioned in the section on previous Raman spectral analysis algorithms, HLA can be applied to our data by skipping step (i), and using the mean spectrum of the analyte calibration set.

Novel Hybrid Algorithm (HLP)
Here we explain our proposed HLP method. It is derived by extending the signal model of Eqn. 1. More specifically, we allow each of the reference spectra S k to vary according to the principal components of variation observed in the background and analyte calibration sets. For each reference spectrum, we penalize deviations from the mean signal in accordance with the eigenvalues obtained from the principal component analyses. In essence, this constrains the variations in the reference spectra to the statistically plausible. Note that the reference spectra can also include the terms of a low-order polynomial background model. However, this is unnecessary since spectrum variations are already modeled by the principal components. In exploratory experiments (not reported here), we verified that the inclusion of a polynomial background model into HLP did not yield further improvements in concentration estimates.
Mathematically, the signal model is extended to where S S k is the mean spectrum observed for compound k, and Z pk is the p th principal component with a non-zero eigenvalue for compound k, observed during prior characterization experiments. Our objective is now to estimate the coefficients W~fw k g and C~fc pk g that give the best signal fit. We do so by using a Bayesian probability framework, where we maximize the posterior probability of W and C, given the measured signal M, the mean reference spectra S S, the principal components Z, and the eigenvalues l (obtained during the same principal component analysis). Using Bayes' theorem and the rules of conditional probability, we can decompose this posterior probability as where N is a normalization constant, and we recognized that P(MD S S,Z,l,W ,C)~P(MD S S,Z,W ,C). We also assumed that S S, Z, W, and C are statistically independent of each other. Note, however, that C and l are not assumed to be statistically independent (see Eqn. 10).
The first term of Eqn. 7 is the data likelihood term. Assuming statistical independence of the samples at each wavelength, it is given by Assuming a zero-mean Gaussian noise model, the probabilities P(m l D S S,Z,W ,C) are given by where s is the standard deviation of the noise in the measured signal. It can be estimated by taking repeated measurements of the same location on a given sample. For the second term in Eqn. 7, we assume that the coefficients c pk are independent, which yields 3.

4.
where P is the number of non-zero eigenvalues. In practice, only the first few eigenvalues and principal components are needed, since they already capture most of the variation seen in the calibration sets. In our experiments, the final concentration estimates were found not to be sensitive to the particular number of principal components used as soon as this number exceeded three or four. By definition, the eigenvalues l pk obtained by the principal component analysis are equal to the variance of the coordinates obtained when projecting all data points on the principal component axis corresponding to l pk , i.e. l pk~s 2 pk . Hence we have Next, the prior probability functions P(l), P( S S), and P(Z) do not contain the variables W and C, and therefore play no part in the optimization problem. Lastly, we assume a uniform distribution for P(W ).
Following standard practice in optimization problems, we optimize the logarithm of the Bayesian cost function given in Eqn. 7. This simplifies the optimization problem by converting multiplications into summations. After dropping constant terms and cancelling common factors, the cost function can be reduced to where b~s 2 . Note that Eqn. 12 takes the familiar form of a penalized maximum likelihood (PML) problem, where b functions as the hyperparameter. The cost function in Eqn. 12 can be efficiently optimized (maximized) by alternatingly solving it as a standard least squares problem in W, and a Tikhonov regularized least squares problem in C. Convergence was observed in all experiments after on the order of 100 iterations. For the experiments in the Results and Discussion section, where L~1015, the time per iteration was 0.02 seconds. All coefficients were initialized with a zero initial guess. The expression for the update of the coefficients w k is, similar to the one in the 'Classical Least Squares Method' section, given byŴ where S {n is the Moore-Penrose pseudoinverse of the n th estimate of the L|K matrix S, which is composed of the signals S lk~ S S lk z P P p~1 c n pk Z lpk , where c n pk are the latest estimates of c pk . To obtain the update steps for the coefficients c pk , it is instructive to substitute Eqn. 6 into Eqn. 12, and to rewrite the latter as where p l~ml { P K k~1 w n k S S lk , and w n k are the latest estimates of w k . To formulate our update step, we store the elements p l into an L|1 vector Q, and the elements w n k Z lpk in an L|KP matrix A~w k Z lpk È É . The matrix C is of size KP|1. Maximizing Eqn. 14 is then equivalent to minimizing the cost function where C~ffi ffiffi b p I s' , and I s' is a KP|KP diagonal matrix that contains the values 1 s pk~1 ffiffiffiffiffiffi ffi l pk p along its diagonal entries. Eqn. 15 is a standard Tikhonov regularized least squares problem, and has the explicit solutionĈ To complete this dicussion, we would like to highlight once again the two key differences between HLA and HLP. First, HLA does not account for variability in the analyte spectrum, whereas HLP does. Second, HLA does not take into account the eigenvalues of the principal component analyses, while they represent important information about the magnitude of the variations observed in the calibration sets. HLP improves on HLA by using the eigenvalues within a Bayesian statistical framework to regularize the weights c pk , and thereby to constrain the allowed variations about the mean spectra to a well-justified statistical range. The resulting improved performance of HLP over HLA is demonstrated in the Results and Discussion section.

Results and Discussion
Here we compare the performance of HLP to that of least squares with a third-order polynomial background model (from hereon referred to as LS-3P), as well as HLA. The third order polynomial was found to give optimal results for Raman signals similar to the ones used here by Lutz et al. [4]. In this paper, we restrict our attention to the case where K = 2. Higher values of K are deferred to future work (see the Further Work section). In other words, the measured spectra are modeled by two reference spectra -one for the analyte (in our case gold-silica nanoparticles), and one for the background. We compare the performance of LS-3P, HLA and HLP for various relative strengths of the analyte and background signal.

Simulation Results
In this section, we simulated the presence of a signal of interest (analyte signal) within a background signal. Our aim was to characterize the accuracy with which the weight (or concentration) of the signal of interest could be recovered, in spite of variability in both the signal of interest and background signal. This analysis was repeated for various dynamic ratios of the two signals. To obtain source signals with a realistic degree of variability to use in the simulation, we collected Raman spectroscopy signals from a real 0.8 nM solution of Raman-enhanced S440 gold-silica nanoparticles produced by Oxonica (now owned by Cabot Security Systems, Boston, MA, USA), as well as signals from a paraffin background material. By performing raster scans across the solution as well as background material, we obtained 106 signals for the S440 nanoparticle solution, and 476 signals for the paraffin background. The collection of signals, the mean signal, and the first two principal components for each are shown in Fig. 1. Note that the bump in the mean paraffin spectrum in the 1600-1800 range appears smaller than could be expected from Fig. 1(b). This is simply due to the fact that the majority of the 476 paraffin spectra contain a smaller bump than the 30 or 40 outlying spectra that visually dominate the plot in Fig. 1(b).
The characterization study was performed by conducting a sequence of experiments as follows. In each experiment, we picked one S440 signal and one paraffin background signal from our database. We then simulated a measured signal by weighting the chosen S440 signal and adding it to the chosen paraffin background signal. This was done for S440 weights of 2 0 , 2 21 , …, 2 212 , and 2 213 . The background signals were not weighted. This was done for all possible combinations of S440 and paraffin signals, at each concentration of S440, yielding a total number of 106|476|14~706,384 simulated signals. For each combination of an S440 and paraffin signal, the remaining signals in the database were used to compute the mean S440 and background signals, as well as their respective principal components. The signal strength (weight w S440 ) of S440 was then recovered by the LS-3P, HLA and HLP algorithms. Finally, the recovered weight w S440 was converted to a concentration estimate c S440 using the formula c S440~c was the concentration of the S440 solution from which the reference spectrum was measured (in this case 0.8 nM).
The performance of each method was evaluated using two metrics: one to measure the closeness of the fit, and one to measure the accuracy of the concentration estimates themselves. The former was the Durbin-Watson statistic of the residual error, as defined by where e l is the residual error at the l th wavenumber. The Durbin-Watson statistic can take values from 0 to 4. A value of 2 signifies no autocorrelation between the successive error values. Values substantially less than 2 indicate positive serial correlation, while those substantially larger than 2 indicate negative serial correlation between the error values. The latter metric was the fractional error, defined as FE~D C est {C true D C true , where C est and C true are the estimated and true concentrations of the S440 signal. An example spectrum, as well as the fitted spectra by LS-3P, HLA and HLP are shown in Fig. 2(a). Note that the fitted LS-3P spectrum shows a wiggle in the 1650-1800 wavenumber range. This is caused by over-fitting of the 3 rd order polynomial background model, which tries to compensate for the variable magnitude of the bump in the 1600-1800 wavenumber range of the paraffin spectrum. As such, the fitted LS-3P spectrum illustrates the limited ability of polynomials to model irregular background variations. While higher order polynomials are able to capture sharper and more irregular background variations, they suffer from increasing degrees of over-fitting and ''whiplash'' effects near the sides of fitted spectra. Lower order polynomials are better behaved, but can account only for slow and smooth background variations. The 3 rd order polynomial model was shown to offer a good compromise between these effects for similar data in [4], but the above mentioned wiggle still reveals its limitations. In essence, the plot illustrates that the variations present in the S440 and paraffin spectra are poorly modeled by a standard polynomial model. From a visual inspection of Fig. 2(a), both HLA and HLP handle the variability more effectively and provide better fits than LS-3P. Fig. 2(b) shows the residual errors of the fits shown in Fig. 2(a). LS-3P clearly yields the greatest residual errors, indicating a failure to model the variations in the component spectra. HLA captures more of the background variation, but not as much as HLP, which shows the smallest residual errors. Fig. 2(c) confirms this trend across all fitted spectra, with superimposed histograms of the Durbin-Watson statistics produced by LS-3P (red), HLA (black), and HLP (green). The histograms show that all three algorithms display positive serial correlation (DW , 2). However, the positive serial correlation of LS-3P tends to be highest (DW values closest to zero), and that of HLP tends to be lowest (DW values closest to two). HLA yields slightly greater degrees of positive serial correlation than HLP, but much less so than LS-3P. Note also that the shape of the DW histograms was similar for all three methods, with each distribution exhibiting two peaks. We found that a particular analyzed spectrum occupied the same relative position in the histograms for each method. In other words, the shifts in the histograms shown in Fig. 2(b) reflect shifts that were true for all of the individual analysed spectra. In short, HLP provided the most right-shifted histogram, and hence yielded the lowest positive serial correlation for all spectra.
Next, Fig. 2(d) shows the estimated concentrations by LS-3P, HLA and HLP, revealing that HLP outperforms both LS-3P and HLA, in the sense that the spread around the true curve is least for HLP. One can also see that HLA produced a large downwards deviation from the true concentration for the second lowest true concentration. This apparent instability yielded 65 values that were lower than the lowest LS-3P value for that concentration. While these 65 cases represent only a small fraction of the total number of signals analyzed, it is worth noting that HLP did not suffer from such instability for any of the analyzed spectra. Next, the curves in Fig. 2(e-f) show the means and standard deviations of the fractional errors as a function of true S440 concentration for both algorithms. The above mentioned instability did not significantly affect HLA's performance curves because the number of signals involved was relatively small. Once again, the HLP algorithm clearly outperforms the LS-3P algorithm at all S440 concentrations, and most markedly so at lower S440 signals/ concentrations. It also outperforms HLA across the range of concentrations tried, doing so most markedly at higher concentrations. Presumably this is because HLP models the variations in the analyte spectrum, whereas HLA does not.

Experimental Results
Variable background phantom. To demonstrate that HLP outperforms LS-3P and HLA on experimental data as well, we designed an experimental phantom where the background signal varied significantly. We placed eight drops of decreasing concentrations of S440 nanoparticles on a thin film of paraffin, which was in turn placed on a background of various colors (see Fig. 3(a)). This background was obtained by printing a color image of a matrix of random numbers between 0 and 1. These printed colors each possessed a distinct Raman spectrum. The first S440 drop had a concentration of 0.8 nM, and subsequent drops were obtained by each time halving the concentration. Raman spectra were acquired on a Renishaw InVia Raman microscope, which was modified by our laboratory for biomedical applications [13]. The integration time of each acquisition was 1 second, and the laser had a wavelength and power of 785 nm and 15 mW, respectively. To characterize the background signal, we first performed a raster scan of the printed color background. The background was covered with a thin paraffin film and had several blank drops of suspension solution (distilled water) placed on top of it. The collection of all acquired signals, as well as the mean background signal and the first two principal components are shown in Fig. 4. The logarithm-transformed (base 2) images of the estimated S440 concentrations are shown in Fig. 3(b-d) for the LS-3P, HLA and HLP algorithms, respectively. The black pixels in the images are points for which the algorithms produced negative weights. It is clear that HLA and HLP succeed in imaging the lowest concentration drop (bottom right), whereas LS-3P breaks down at this low concentration.
To further examine the difference between the three methods, and between HLA and HLP in particular, we provided fitted spectra as well as quantitative results in Figs. 3(e)-(i). Fig. 3(e) shows example fitted spectra computed by LS-3P, HLA and HLP. Fig. 3(f) shows the residuals, obtained by subtracting the fitted spectrum from the measured spectrum. In this example, HLP shows hardly any serial correlation in the residual error (Durbin-Watson = 1.999), unlike LS-3P, which shows clearly structured residuals (Durbin-Watson = 0.127). It also compares favorably to the Durbin-Watson statistic of 1.621 by HLA. Fig. 3(g) shows that the favorable Durbin-Watson statistic of HLP over LS-3P and HLA generally holds for all spectra. HLP achieves values closest to 2, which indicates zero serial correlation. HLA again gives the second best performance. In other words, HLP best captures the variability in the background signals, translating into the closest fits, and the least amount of serial correlation in the residual errors. Note also that the distribution of the Durbin-Watson statistic takes a different shape from the one seen in Fig. 2(c); Fig. 3(c) shows only one peak, whereas Fig. 2(c) shows two peaks. The difference is merely due to the presence of different background materials. The conclusion is however unchanged: HLP right-shifts the entire distribution towards a DW value of 2.
An evaluation of the accuracy of each method in terms of their concentration estimates is provided by Figs. 3(h) and (i). The former shows the estimated concentrations of the droplets; the latter shows the percent error in the estimated concentrations. The concentration estimates were computed by taking the average of the 6 brightest pixels within each droplet. The curves show that LS-3P, HLA and HLP perform very similarly for high nanoparticle concentrations. The performance of LS-3P deteriorates the most rapidly of the three algorithms as the nanoparticle concentration decreases. HLA and HLP yield similar results, though HLP yields better estimates towards the lower nanoparticle concentrations, i.e. when the nanoparticle signal strength is low relative to the background signal. Note that the lowest weight was not shown for LS-3P since it was negative, and hence does not have a real-valued logarithm.
Ex vivo pig colon experiments. In this experiment, we placed S440 drops of decreasing concentrations onto an excised pig colon as shown in Fig. 5(a). All drop locations are identified by a red circle. The pink dots within the red circles reveal the locations of the most concentrated S440 drops. The drops were produced by starting with a stock concentration of 0.8 nM, and serially diluting by a factor of 2 for every subsequent drop. Drops were arranged in 3 columns of 5 drops, starting on the top left of the sample. The drop concentrations decreased down the columns. Next, we raster scanned the sample with a 1 mm step size, a 1 second integration time, and a 785 nm laser with a power  of 15 mW. The total time for this scan was approximately 1 hour. Prior to the placement of the drops, the pig colon was also raster scanned at a step size of 2 mm to obtain a library of background signals. The obtained background signals, as well as their mean and first two principal components are shown in Fig. 6. The objective of the experiment was to compare the S440 concentrations computed by LS-3P, HLA and HLP. The calibration spectra for S440 were obtained at the stock concentration of 0.8 nM. Figs. 5(b-d) show the concentration maps obtained by LS-3P, HLA, and HLP, respectively. HLP and HLA both eliminated the false positive weights generated by LS-3P in the tissue background where no S440 is present. However, HLP captured more of the lowest concentration drop than HLA. Figs. 5(e-f) show example fitted spectra with their corresponding error residuals. It is clear that HLP produced the closest fits, as evidenced by the lowest amount of structure in the error residuals. Fig. 5(g) shows a histogram of the Durbin-Watson statistics for all fits by each algorithm. The HLP histogram shows values that are closest to 2, confirming that HLP's closer fits held widely across all measured spectra.
Lastly, Fig. 5(h) plots the S440 concentrations computed by each algorithm. The concentrations were computed by taking the mean concentration of the 5 brightest pixels in each drop. Note that the S440 particles were not detectable beyond the seventh drop for LS-3P, and not beyond the 8th drop for HLA and HLP. The concentration estimates computed by HLP adhered most closely to the expected linear relationship (as shown by the blue 'Theoretical' line). As before, HLA gave the second best performance, and LS-3P gave the poorest performance.

Conclusions and Further Work
Conclusions. In this work, we showed that the LS-3P algorithm for Raman spectroscopy is sensitive to natural variations in the reference spectra. Our compound of interest was a solution of Raman-enhanced gold-silica nanoparticles. We found that the nanoparticle concentration estimates were sensitive to spectrum variability primarily when the nanoparticle signal (concentration) was weak. We proposed a novel algorithm (HLP) that is more robust to variations in the reference spectra. The HLP method was compared to both LS-3P and Berger et al.'s HLA method [12], where the latter was specifically designed as an improvement over PCR and PLS, and to be robust to variability in the background signal. HLP's superior performance over LS-3P and HLA was shown for both simulated and experimental data. The simulated data was generated by digitally combining weighted instances of experimentally obtained nanoparticle and background spectra. The experimental data was obtained from serially diluted in vitro solutions of Raman-enhanced gold-silica nanoparticles. We confirmed our expectation that Berger et al.'s HLA method significantly outperforms LS-3P. However, our HLP method improved concentration estimates even further for two reasons. First, unlike HLA, it accounts for variability in the reference spectrum of the nanoparticle itself, not just in background signals. This is a useful property since, in practice, a single reference  spectrum will always contain inaccuracies. Second, HLP incorporates the eigenvalues of the principal component analysis. This regularizes how far the reference spectra are allowed to 'stray' from their means, in accordance with experimentally observed variability. Both of these modeling improvements result in the improved performance of HLP compared to HLA and LS-3P, as demonstrated by the experiments in this paper. We also note that our performance metrics (Durbin-Watson statistic and fractional error) were computed over many spectrum instances (N = 1461066476 = 706,384 for the simulated data, N = 180 for the color background experiment, and N = 1,653 for the pig background experiment). In other words, our conclusion that HLP outperforms HLA is strengthened by the fact that the improvement was observed across many different spectra. Furthermore, the improvement was demonstrated across three different experimental set-ups, lending further credibility to our results. Lastly, the improvement makes intuitive sense because of the additional information exploited by HLP compared to HLA, namely the eigenvalues of the background calibration set, and the principal components and eigenvalues of the analyte calibration set.
Further work. In this paper, the HLP algorithm was evaluated for the case where K = 2, i.e. where only a single analyte and mean background were considered. In future work, we will examine the performance of HLP for multiplexed spectroscopy (K.2). Second, the HLP algorithm presented in this paper assumed a Gaussian noise distribution. In our future work, we will evaluate the merits of using a Poisson instead of Gaussian noise model, in an effort to decrease the lowest detectable nanoparticle concentration. Lastly, we recall that the HLP algorithm is equally capable of modeling variations in the background spectrum as in the spectra of compounds of interest such as gold-silica nanoparticles. While the S440 gold-silica nanoparticle spectra were found to be relatively stable in this study, other studies have reported significant variability in the spectra of nanoparticle solutions (eg. [14]). A common cause for such variability is non-uniformity in the particle sizes. In future work, we will examine HLP's ability to improve the robustness of the concentration estimates for such nanoparticles with less stable spectra. Lastly, we aim to investigate the spectral variability of biomolecules as well, due to effects such as conformational changes. As before, we will assess the extent of those variations, as well as the robustness of our novel HLP method to them.


extfieldsofstudy: Medicine Chemistry

text:
Inorganic–organic hybrid materials through post-synthesis modification: Impact of the treatment with azides on the mesopore structure

Hybrid, hierarchically organized, monolithic silica gels, comprising periodically arranged mesopores and a cellular macroscopic network, have been prepared through a co-condensation reaction of tetrakis(2-hydroxyethyl)orthosilicate with chloromethyl-trimethoxysilane or 3-(chloropropyl)-triethoxysilane. Subsequent conversion of the chloro groups into azido groups, by nucleophilic substitution with NaN3 in N,N-dimethylformamide, was conducted upon preservation of the monolithic structure. However, treatment with NaN3 had a strong influence on the structure in the mesoporous regime, with changes such as an increase of mesopore diameter, pore volume and lattice constants, as well as a concomitant decrease of the pore wall thickness, as confirmed by small angle X-ray scattering, transmission electron microscopy, and nitrogen sorption analysis. Similar effects were observed for unmodified silica gels by simple ageing in azide-containing media, whether a relatively small or a sterically demanding counter ion (Na+ or (H3C)4N+) was used. The structural modification did not seem to depend greatly on whether an organic aprotic solvent (N,N-dimethylformamide, 1,1,3,3-tetramethylurea, 1,3-dimethyl-2-imidazolidinone) or a protic solvent that can form hydrogen bonds, such as water, was used.

: Schematic representation of a hexagonally organized pore system and the characteristic sizes From a geometrical point of view, one can consider the lattice constant a to be: a = D + t (1) with D being the mean pore-diameter and t the mean wall thickness between the pores. Variation of pore diameter and wall thickness leads to changes in the electron density distribution and thus to differences in the scattering intensities of the diffraction peaks.
According to Glatter and Kratky [1], the scattering intensity is proportional to the product of structure factor S(q) and the square of the scattering amplitude (form factor) P(q): For the case of a 2-D hexagonal lattice with spherical pores, S(q) is described as the sum of Gaussian functions with maxima located at the reciprocal lattice points. The reciprocal lattice points are equal to the positions of the  [2,3].
The form factor P(q) can be written as [4]: with R being the mean pore radius and J1 the Bessel function of the first kind.
The scattering intensity (2) was then convoluted with the resolution function of the equipment described by a Lorentzian curve and fitted to the experimental values. Fitting parameters were the pore diameter and the wall thickness as described above.

Electron density reconstruction:
A different approach is the reconstruction of the electron densities from a Fourier series and by the appropriate choice of the phases [5,6]. This has been experimentally and theoretically used to model the electron density across the pore for modified and unmodified MCM-41 and SBA-15 materials [7,8]. For a centrosymmetric unit cell, the electron density can be written as a Fourier series of cosines: With Aq being the Fourier coefficients, q  the scattering vector of the reciprocal lattice and r  a vector to the position in the unit cell.
The Fourier coefficients are determined from the peak intensities of the observed (h,k) reflections by , LP(qhk) being the Lorentz polarization-factor [9] and m the multiplicity of the respective reflection [5,6].
In order to obtain the integrated intensities I(q), the background scattering was subtracted from the data using a non-linear baseline. The remaining peaks were then fitted by a Gaussian function, from which the peak area was determined.
The next step in reconstructing the electron density is choosing the right phase. For centro-symmetric structures, this is equivalent to choosing the correct sign (+ or -) of the Fourier coefficients. The method, known as the "swelling method", was originally described by Stamatoff and Krimm [10] for lamellar phases ( α L ) and was extended by Turner and Gruner [5] to hexagonal phases: One compares two phases with slightly differing radii and therefore slightly differing X-ray patterns. Then all possible permutations of signs for the Aq's are calculated and compared, the proper phase set will have the closest correspondence (for more details see [5,6]). In our case, we did not use a swelling of phases, but compared the samples before and after nucleophilic substitution, which led also to slightly differing X-ray patterns. Finally, one has to check if the resulting phase is physically consistent. For example, if the origin was placed in the centre of the pore cross section, the electron density should be a minimum at this point. Towards the pore wall, the electron density should increase and reach a maximum in the silica phase between the pores.    Figure S3: Nitrogen isotherms and SAXS patterns of untreated silica gels, reference silica gels (solvent/60 °C) and azide-treated silica gels (solvent/NaN3/60 °C) in different solvents: TMU (top) and H2O (bottom) Figure S4: Electron-density reconstruction for unmodified silica gels exposed to DMI and DMI/NaN3 (top left), TMU and TMU/NaN3 (bottom left), H2O and H2O/NaN3 (bottom right) and for gels exposed to DMF and DMF/TMAA (top right).


s2fieldsofstudy:
Chemistry
Materials Science

text:
Synthesis of Different Zinc and Zinc Included Nanostructures by High Power Copper Vapor Laser Ablation in Water- Surfactants Solutions

The data of experimental studies of optical characteristics of colloidal solutions, composition and morphology of its dispersed phase, resulting from laser ablation of zinc in aqueous solutions of anionic surfactants --- sodium dodecyl sulfate (SDS), dioctyl sodium sulfosuccinate (AOT) are presented. It is shown that by studying the optical absorption spectra of the colloid, X-ray spectra and AFM-images of extracted from colloid solid phase, it is possible to trace the dynamics of ZnO nanostructures formation from zinc nanoclasters size of several nanometers to ZnO fractal aggregates (FA) size up to hundreds of nanometers. Determinants of this process are the average power and an ablation exposure, the frequency of the laser pulses, the colloid aging time, the type and concentration of surfactant in solution. In the selection of appropriate regimes, along with zinc oxide obtained other nanoproducts --- hydrozincit and organo-inorganic layered composite \ce{[(\beta) - Zn(OH)2 + SDS]}.

Introduction
Over the past decade, among chemical, electro physical and others methods of preparation nanoparticles and nanostructures of metals, oxides, nitrides,the new technique-pulse laser ablation in liquid environment has demonstrated intense interest. 1,2 This method is quite simple, versatile and promotes to have the pure product.A great number of experimental works are devoted for the preparation and investigations of ZnO nanostructures by this method. [3][4][5][6][7][8] Zinc oxide, semiconductor II-VI group, with a wide band gap (Eg ≈ 3.37 eV ), high exciton binding energy (≈ 60 meV ) at room temperature is a very interesting material for optoelectronics (UV lasers, diodes, different models of sensors et al). [9][10][11] The analysis of researches performed on this problem shows that the majority of them were carried out using a low-power Nd:YAG lasers (≈ 1 W ) with low pulse repetition rate (≈ 10 Hz). At the small energy deposition in the ablation process a yield nanoproducts is negligible. The intensification of the process can be achieved by increasing the average laser power output and increase the interaction time with the target material and formed during the ablation colloid, i.e., with increasing of laser pulse repetition rate. In a recent paper, 12 authors using a laser with an average output power ≈ 25 W , picoseconds pulse duration (7 ps) and repetition rate up to 200 kHz, demonstrated the importance exclusion overlapping zones interaction (bubbles) between successive pulses. It is shown that in the case of overlapping zones interactions, efficiency of production of nanoparticles in this experiment decreased by 30-50 percent and more due to defocusing and scattering radiation incident on the target. At an appropriate choice of spatial and temporal distribution of radiation zones on the target, the performance the synthesis of nanoparticles can be increased at hundred-thousand times in comparison with the data of the above works with lowpower Nd:YAG lasers. [3][4][5][6][7][8] In this regard, the study of the nanostructures synthesis by method the ablation different materials in liquids with using high-power lasers with high pulse repetition rate opens up the prospect of effective nanostructures synthesis of new forms and qualities.

Experimental Section
For the experiments was chosen copper vapor laser generated pulses of radiation at wavelengths of 510.6 nm (≈ 60% power output) and 578.2 nm with an average power 15-20 W. The pulse width-20 ns, repetition rate -10 kHz.The radiation was focused on a target located at the bottom of the glass vessel. The focal length of lens -280 mm, a diameter of spot radiation on a surface zinc (≈ 99.5%) target less 100 µm. The vessel was contained deionized water or aqueous solutions of anionic surfactants -sodium dodecyl sulfate (SDS -C 12 H 25 SO 4 Na) and dioctyl sodium sulfosuccinate (AOT -C 20 H 37 NaO 7 S). In the case of AOT was used 0.15 M solution AOT in Nonan (C 9 H 20 -organic compound classes of alkanes) followed by dilution with water up to molar ratios water / AOT -W ≈5, 60, 700. The peculiarity of this solution is the formation of a reverse micelle, the volume of the water pool which depends on W. The volume of the water pool is limited to the size of the nanoparticles. 13 The volume of solutions was about 2 cm 3 . The vessel was rotated during ablation time to change laser beam position on a target surface. The solution temperature was measured at a process of ablation. Laser ablation lasted for 30 min and 3 h.
The optical absorption spectra of obtained colloids were recorded after different aging time by UV-Vis spectrometer SF-46 LOMO. The composition and structure of the solid phase colloids studied by X-ray diffraction (XRD; Dron-2, Cu Kα radiation). Morphology of the precipitate was analyzed by scanning electron microscopy (SEM; Hitachi S405A, 15 kV) and atomic force microscopy (AFM; Solver P47-PRO, method of semicontact topography). The samples for these studies were obtained after a centrifugation at 4000 round/min all of colloids and sediments were washed with deionized water several times to remove main part of surfactant and dried at ≈ 60 0 C temperature.  10-15 nm. In the spectra of other colloids, this band is not observed. In addition, a number of spectra of colloids, the obtained in deionized water, with SDS and AOT, you can see the individual peaks and the characteristic rise of absorption curves in the wavelength range 280-300 nm, which indicates the presence of zinc oxide particles size 1.5-5 nm. 14,15 The peculiarity of the absorption spectra of almost all investigated in this study of colloids is the appearance of a rather narrow ( λ ≈ 20 nm) zone of sharp decrease in the level of absorption (peak of transparency) in the wavelength range 450-480 nm. Peak amplitude increases as with increasing exposure time, and with increasing concentration of SDS and AOT in the initial solution.

Results and Discussion
It is also noteworthy that the peak amplitude decreases significantly with aging time and the peak completely disappears after centrifugation of the colloid fig.2, curve 4, fig.3, curve 2. This means that the appearance of peaks is caused the formation in the solution of fairly larges complexes of particles, which leave it by sedimentation or centrifugation. Interpretation of the absorption spectra of colloids is a difficult task. Indeed, the solid phase colloid consists of products of chemical interaction of zinc with water and surfactant -a nanoparticles of zinc and zinc oxide hydrates, organic-inorganic composites produced in the process of ablation and the subsequent aging of the colloid. All components of the solid phase at different times can change their size and structure due to processes of growth, aggregation and coagulation. Therefore, the overall picture of the optical spectra of a colloids is a superposition of the spectra of plasmon resonances of different nanoparticles, of peaks and bands of exciton absorption of zinc oxide, the absorption bands of other components of colloids. Nature of the changes in the spectra obtained in the present experiment, their features, the authors believe, can be qualitatively explained in terms of physics of fractals, in particular concerning the optical properties of colloidal systems and composite nanostructures, which are typical examples of fractal systems. 21 The theory of optical properties of fractal aggregates (FA) is different from the well-known Mie theory for metal nanoparticles, since taking into account the electrodynamic interaction of neighbors of the particles that make up the fractal aggregate. Corollaries to this theory which confirmed experimentally are a broadening and frequency shift of the spectra of plasmon resonances of a fractal, i.e., the increase of length of the long-wavelength part spectra in comparison with a spectra system of individual, non-interacting particles.Frequency shiftω − ω p, i.e. a mutual electrodynamics influence of neighboring particles of fractal is so great that it leads to a shift in their resonance on amount comparable to the resonant frequency. The length of the long-wavelength range of spectra increases with the size of a fractal. From the theory of optical properties of fractals, in particular, illustrated by the example of silver colloids, it should be possible appearance in the FA spectrum additional peaks and dips in the long-wavelength range of spectra associated with changes in the distribution of particle sizes, including the presence of surfactants in the colloid. There is also a photodynamic effect on colloid particles, the resonances are close to the frequencies of the incident radiation. 22 This effect leads to rather narrow dips -"burn-through" -in the absorption spectra of the colloid.
In these experiments, the above features of the spectra -the rise of the level of absorption especially in the long side of spectra, the emergence of a broad peak in the wavelength range 470-500 nm, the peaks of "transparency" and their disappearance during the sedimentation of large complexes and centrifugation -can be attributed to the provisions of the fractal theory .
Thus, the received optical spectra allow to assume that fractal aggregates of various compositions and shapes are formed in colloids. It proves to be true by X-ray diffraction patterns and an AFM-images of the solid phase extracted from colloids. Fig.4 shows modifications of precipitate X-ray spectra in dependence on a kind and concentration of a surface-active agent in the initial solution. The X-ray spectrum 1 corresponds to the precipitate received after bombarding radiation of zinc in deionized water at τ e = 3.5 h and τ st = 40 h. This spectrum contains accurately expressed peaks belonging to crystalline zinc oxide. This phase has hexagonal crystalline structure with cell parameters a = 0.3294 ± 0.0001 nm, c = 0.5214 ± 0.0001 nm. The average size of crystallites along the crystallographic axis a is L a = 41 nm, along the axis c -L c = 84 nm. In this spectrum zinc peaks are not visible.At adding of 0.001 M SDS in the initial solution two crystal phases are exhibited in spectrum 2. One of them -ZnO (a = 0.3252 ± 0.0001 nm; c = 0.5199 ± 0.0001 nm).
The crystallites sizes along the crystallographic axes a and c are 40 nm. This is less than in experiment with deionized water. The second crystalline phase is presented in the spectrum by a set of the diffraction lines to the greatest degree coinciding with spectral lines of hydrozincite The X-ray diffraction studies fulfilled in the present operation and the data of other exper-iments 11,13,25 show the following dynamics of zinc basis nanoparticles synthesis depending on concentration of anionic surface active agents -SDS and AOT. Only Zn and ZnO are formed at ablation in deionized water. At concentrations being essentially smaller than critical value (0.001 M SDS) hydrozincite is received in the conditions of powerful bombarding radiation by means of copper vapor laser. At achievement of SDS concentration about 0.01 M a layered organic-inorganic composition material is formed in experiments with copper vapor laser and at the durable affecting (∼ 1h) of low-power Nd: YAG-laser radiation. 24 The further rise of SDS and AOT concentration leads to formation of finely dispersed yields of zinc ablation. This fact shows that surface-active agents are able to limit nanoparticles growth testifies to developing process of properties. At the same time increase of surfactant concentration in a solution promotes growth and a union of micelles in fractal frames.

Conclusion
The measured optical spectra, diffraction patterns and AFM-images reflect the dynamics of the formation of nanostructures, resulting from the experiment, a powerful copper vapor laser. The increase in the total time of ablation, the duration of the effects of radiation on the colloid at high repetition rate, rise in temperature of the colloid due to the high average power of irradiation lead to an intensive operating time of nanoparticles of zinc and zinc oxide hydrates, the emergence of clusters and large (up to hundreds of nanometers, and more ) complexes with different structures and forms -fractal aggregates. Based on fractal theory may explain the qualitative features of the optical spectra of the studied colloids. Specifics of the use of surfactants in the process of ablation is manifested in the formation of micelles, consisting of surfactant molecules and nanoparticles.
On the one hand, the surfactant molecules surrounding the nanoparticle, limit their growth and aggregation, on the other hand, there is a possibility of synthesis of various chemical compounds on the basis of surfactant, water and metal that occurs in these experiments. As follows from literature data obtained organic-inorganic composite [(β )−Zn(OH) 2 + DS] is of interest as a promising equality botchy stochastic environment for UV lasers. [26][27][28]


extfieldsofstudy: Chemistry Medicine

s2fieldsofstudy: Chemistry Biology

text:
One-Pot Synthesis of Novel 2,3-Dihydro-1H-indazoles

A copper(I)-mediated one-pot synthesis of 2,3-dihydro-1H-indazole heterocycles has been developed. This synthetic route provides the desired indazoles in moderate to good yields (55%–72%) which are substantially better than those achievable with an alternative two-step reaction sequence. The reaction is tolerant of functionality on the aromatic ring.

Additionally, there have been a few studies in which formation of derivatives of indazole 5 have been observed as a result of trapping experiments on reactive intermediates [9][10][11]. However, other than Zenchoff's limited work, we are unaware of the development of any general routes for the synthesis of this class of heterocycles.
In addition to being interesting compounds in their own right, 2,3-dihydro-1H-indazoles 5 have also served as synthetic intermediates for their more thoroughly-studied counterparts, the 1H-indazoles 4 [8]. Thus, development of synthetic routes towards the synthesis of indazoles 5 similarly creates novel routes for the synthesis of indazoles 4.
The sparse number of investigations into the synthesis of this class of heterocycles is, therefore, quite surprising given that 2,3-dihydro-1H-indazoles have such strong potential for biological activity and as convenient synthetic intermediates. Detailed investigations of their properties, however, demand robust synthetic methods for their preparation. Herein, we describe our initial investigations towards the synthesis of this understudied class of heterocyclic compounds.

Results and Discussion
Ullman-type copper(I)-mediated coupling of bis-BOC protected hydrazine 8 to aryl and vinyl halides in DMF using CuI, 1,10-phenanthroline and Cs 2 CO 3 is a well-established synthetic procedure (see Scheme 2) [14,15]. We initially envisioned the synthesis of bis-BOC protected 2,3-dihydro-1Hindazoles 10 via intramolecular coupling of the free N-H bond of hydrazines 9 at the iodo-substituted carbon atom to form the five-membered indazole nucleus (Scheme 3). Indeed, this proved to be an excellent method for the synthesis of substituted indazoles, providing high yields of 10 (Scheme 3).

Scheme 2.
Ullman-type coupling of hydrazine 8 to halogenated aromatics. Unfortunately, however, synthesis of the required trisubstituted hydrazine starting materials 9 proved to be problematic (Scheme 4). Although we utilized reaction conditions optimized by Rasmussen for selective monoalkylation of 8 [16], only poor yields of monoalkylated product were obtained (i.e., compounds 10a-c were obtained in 25%, 27% and 16% yields, respectively). Surprisingly, bisalkylation of hydrazine 8 was the preferred reaction route. We attempted to increase the yield of monoalkylated product by increasing the amount of starting 8 relative to starting benzyl bromides 11, but this led to difficulties in separation of the product from excess 8. Additionally, other than 11a, which was commercially available, the required substituted ortho-iodobenzyl bromides (compounds 11b-f) needed to be synthesized from the corresponding benzylic alcohols via a two step process of iodination of the aromatic ring (silver CF 3 CO 2 Ag/I 2 ) followed by bromination at the benzylic alcohol position (PBr 3 ). Thus, although the intramolecular coupling reactions afforded high yields (Scheme 3), the effective yields from the synthetically expensive ortho-iodobenzyl bromides were unsatisfactorily low.

Scheme 4. Synthesis of monoalkylated bis-BOC protected hydrazines 9.
To bypass the problem of dialkylation, we considered the possibility of developing a one-pot procedure in which copper-mediated coupling of the initially formed monoalkylated hydrazines 9 to form indazoles 10 might be able to compete with the complicating dialkylation process. A mixture of 8 and ortho-iodobenzyl bromide 11a was added slowly, via syringe pump (over a period of 7 h), to a pre-heated stirring mixture of CuI, 1,10-phenanthroline and Cs 2 CO 3 . After 24 h, we were gratified to note that the only major product detected by TLC analysis of the crude reaction mixture was the desired indazole 10a. Reaction workup, followed by column chromatography afforded 10a in 60% yield, a substantial increase from the effective 22% yield obtained via the two-step process (i.e., 11a→9a→10a) starting from the ortho-iodobenzyl bromide 11a. In addition, this procedure avoided the need for purification of the intermediate hydrazine 9a. Similarly, indazoles 9b and 9c were formed in 60% yields utilizing the one-pot procedure rather than the effective 25% and 15% yields, respectively, from the corresponding two-step procedures.
Given the success of this one-pot procedure, we subjected a number of substituted ortho-iodobenzyl bromides to the same reaction conditions. The yields were consistent at ~60% (Table 1), and the products were obtained in pure form following column chromatography. The reaction was tolerant of electron-donating (i.e., OCH 3 ) and electron-withdrawing (i.e., CO 2 CH 3 ) groups.

General
All chemicals and solvents were used as received (Aldrich, St Louis, MO, USA) Anhydrous DMF was kept under nitrogen and sealed with a septum. Column chromatography was performed using 230-400 mesh silica gel 60. NMR spectra were recorded on a Varian 60 MHz instrument in CDCl 3 as solvent, unless otherwise indicated, and referenced relative to TMS (0.0 PPM). Combustion analysis was performed by Micro Analysis Inc. (Wilmington, DE, USA). Other than unsubstituted 2-iodobenzyl bromide 11a, which was commercially available (Aldrich), benzyl bromides 11b-f were synthesized via standard aromatic iodination (I 2 , CF 3 CO 2 Ag) followed by bromination of the benzylic alcohol (PBr 3 ) following procedures described in the literature [17].

Conclusions
In summary, a convenient synthesis of substituted 2,3-dihydro-1H-indazoles has been developed. This one-pot procedure provides indazoles at considerably higher yields than the corresponding twostep process. This constitutes one of the first general synthetic methods for this class of heterocyclic compounds.